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which  provide  precise  command  response  over  the  entire  flight 
regime  while  augmenting  the  airframe  stability  characteristics. 
The  design  of  an  air  combat  digital  flight  control  system  for 
the  F-14A  is  described  in  this  report,  and  its  performance  is 
verified  by  the  use  of  a  nonlinear  six-degree-of-freedom  aircraft 
simulation  and  an  air  combat  pilot  model.  Comparison  of  the 
^tracking  accuracy  of  the  DFCS- aircraft  to  that  of  a  conventionall 
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line  and  reduces  the  likelihood  of  control-induced  departure.  In 
"the  realm  of  jsilot  model  validation,  a  method  of  choosing  the 
best  model  among  those  hypothesized  is  described  and  its  operatic 
demonstrated.  Piloted  simulation  data  are  processed  which  indica 
that  the  pilot  does  not  change  his  control  strategy  to  correspond 
to  changes  in  the  aircraft  flight  condition. -Not  only  does  this 
validate  the  pilot-aircraft  stability  boundatV calculation  method 
but  it  has  significant  control  system  design  implications,  since 
design  of  control  system  with  uniform  handling  qualities  is  shown 
to  be  important.  The  methods  and  results  presented  in  this  study 
can  have  substantial  impact  on  the  development,  analysis,  and 
testing  of  high  performance  aircraft,  enhancing  safety ,  reliabili 
and  effectiveness  of  flight  operations.  \ 
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INTRODUCTION 


1 . 


1 . 1  BACKGROUND 

Highly  maneuverable  aircraft  must  achieve  good  han¬ 
dling  qualities  over  a  wide  range  of  flight  conditions.  The 
aircraft  control  system  designer's  task  is  difficult  because 
the  flight  conditions  for  such  aircraft  involve  high  accelera¬ 
tions  and  rapid  rolling.  Additionally,  significant  changes  in 
control  power,  or  even  control  reversals,  can  occur  at  some 
flight  conditions.  The  control  system  developed  for  such  air¬ 
craft  must  be  a  command  augmentation  system  (cas),  which  has 
two  functions  to  perform.  The  first  is  to  act  as  a  stability 
augmentation  system  (sas),  which  modifies  and  damps  the  natural 
fast  modes  of  motion  of  the  aircraft  (short  period,  Dutch  roll 
and  roll  mode)  so  that  the  pilot  has  a  "well-behaved"  platform 
from  which  to  prosecute  the  engagement.  The  second  control 
system  function  consists  of  control  response  augmentation. 

The  control  system  must  modify  the  basic  control  response  of 
the  aircraft  so  that  the  response  appears  uniform  over  the 
range  of  flight  conditions.  Lateral  stick  should  always  pro¬ 
duce  similar  aircraft  roll  response,  longitudinal  stick  should 
always  produce  similar  (although  veloci ty- scaled )  normal  ac¬ 
celeration  response,  and  pedals  should  always  produce  similar 
yawing  motion.  This  uniformity  of  response  is  important,  since 
the  pilot  is  under  a  severe  work  load  in  maneuvering  flight 
and  should  not  be  forced  to  attempt  to  adapt  his  input  strategy 
to  the  specific  aircraft  flight  condition  --  a  difficult  task 
at  best. 


The  complexity  of  the  command-augmentation  task  im¬ 
plies  a  digital  implementation  of  the  aircraft  control  system, 
specifically  a  full-authority  digital  flight  mtrol  system 
(DFCS).  In  principle,  the  pilot's  control  stick  and  pedal 
inpurs  drive  the  DFCS  computers,  which  then  command  the  con¬ 
trol  surface  deflections.  Previous  work  has  been  concerned 
with  stability  and  control  analysis  techniques  (Refs.  1  and  2) 
and  control  system  design  techniques  (Refs.  3  and  4)  that  are 
intended  to  provide  solutions  to  the  control  system  designer's 
difficulty.  This  report  details  work  done  to  extend  and  verify 
these  analysis  and  design  techniques.  Figure  1  outlines  the 
relationships  among  the  flight  vehicle  analysis  tools  and  illus 
trates  the  two  portions  of  the  program  examined  in  this  report: 
control  system  construction  and  validation  using  a  nonlinear 
simulation,  and  pilot  model  hypothesis  testing. 

As  illustrated  in  Fig,  1,  the  performance  of  the  im¬ 
proved  DFCS  designs  developed  in  this  report  are  evaluated 
with  a  nonlinear  6-DOF  aircraft  model.  An  aircraft  response 
comparison  is  made  with  the  same  pilot,  inputs  but  with  differ¬ 
ent  control  systems:  the  DFCS  and  the  conventional  control 
system  (with  and  without  analog  augmentation).  Maneuvers  are 
also  simulated  with  a  conventional  mechanical  control  system 
with  the  DFCS  driving  only  the  limited-authority  sas  actuators. 
The  latter  implementation  of  the  DFCS  has  promise  for  the  retro 
fit.  of  present-day  aircraft,  and  it  has  a  major  impact  on  the 
reliability-redundancy  issue. 

Previous  work  (Ref.  3)  also  demonstrated  the  useful¬ 
ness  of  optimal  control  pilot  models  in  determining  pilot/ 
aircraft  stability  boundaries.  While  Ref.  3  shows  that  these 
stability  boundaries  compare  favorably  with  actual  pilot  re¬ 
sponse,  no  attempt  was  made  there  to  verify  the  details  of  the 
pilot  model  predicted  response.  This  report,  considers  a  method 
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COMPLETE 


Figure  1  Flight  Vehicle  Analysis  Methods 


that  can  be  used  to  verify  pilot  modeling  techniques:  hypothe¬ 
sis  testing.  In  this  method  the  pilot  model  that  best  matches 
the  data  among  a  set  of  hypothesized  pilot  models  is  chosen. 

In  the  application  presented  in  this  report,  data  from  a  full- 
order,  nonlinear  piloted  simulation  (the  Differential  Maneu¬ 
vering  Simulatior  at  NASA  Langley  Research  Center)  is  used  to 
test  the  accuracy  of  the  predictions  of  several  different,  op¬ 
timal  control  pilot  models. 
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1.2 


ORGANIZATION  OF  THIS  REPORT 


The  research  contained  in  this  report  draws  heavily 
on  past  work  TASC  has  performed  for  ONR  (Ref.  3).  Chapter  2  of 
this  report  details  the  four-step  DFCS  design  procedure  that 
has  been  developed  at  TASC;  its  application  to  the  design  of  a 
control  system  for  the  F-14A  is  also  discussed  there  .  The 
DFCS  is  tested  in  Chapter  3  b.  comparing  its  performance  with 
that  of  a  conventional  control  system  in  accurately  controlling 
the  aircraft  along  rapid  ACM  trajectories.  Validating  the 
prediction  of  piloting  control  strategy  using  a  pilot  model  is 
addressed  in  Chapter  4,  where  a  hypothesis  testing  procedure 
is  discussed.  Actual  validation  and  testing  with  piloted  simu¬ 
lation  data  is  included.  The  conclusions  of  this  report  and 
recommendations  concerning  useful  extensions  of  this  work  are 
given  in  Chapter  5. 


^Figure  1  illustrates  the  organization  of  the  report. 
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2. 


AIR  COMBAT  DIGITAL  FLIGHT  CONTROL 
SYSTEM  DESIGN 


Two  digital  flight  control  systems  (DFCS)  for  a  high¬ 
ly  maneuverable  fighter  aircraft  are  designed  and  described  in 
this  chapter.  The  objective  is  to  develop  and  demonstrate 
digi tal  flight  control  design  technology  which  can  improve  the 
performance  and  mission  effectiveness  of  a  fighter  aircraft. 

Based  on  coupled  dynamic  models  of  the  aircraft,  the  two  de¬ 
signs  are  formulated  using  modern  control  techniques  as  a  Type 
0  and  Type  1  DFCS.  The  designs  occur  concurrently,  since  the 
control  laws  are  related  through  an  algebraic  transformation. 
However,  their  implementation  and  response  to  any  disturbances 
and  modeling  errors  are  different.  These  design  differences 
are  also  discussed  here. 

The  two  DFCS  designs  have  a  fly-by-wire  configuration 
and  can  be  easily  implemented  in  a  digital  flight  computer. 

The  three  pilot  commands  used  in  the  two  DFCS  are  the  same  as 
in  Ref.  3  (normal  acceleration,  stability-axis  roll  rate,  and 
sideslip  angle)  and  are  transformed  by  each  control  law  to 
five  separate  control  surface  deflections.  This  unique  feature 
of  the  controllers  whereby  each  control  surface  is  deflected 
independently  in  a  optimal  coordinated  fashion,  will  figure 
prominently  in  many  of  the  discussions  in  the  following  sections. 

2.1  TYPE  0  AND  TYPE  1  DIGITAL  CONTROLLERS 

A  Type  0  DFCS  tracks  constant  commands  without  using 
a  pure  integration  in  the  forward  loop.  The  Type  0  DFCS  relies 
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on  proper  choice  of  the  matrix  which  feeds  forward  the  command 
to  achieve  tracking.  A  Type  1  DFCS  tracks  constant  commands 
with  zero  steady-state  error  by  using  a  pure  integration  in 
the  forward  loop.  The  Type  1  DFCS  asymptotically  estimates 
the  proper  vector  to  feed  forward  to  the  control  to  maintain 
tracking. 


Procedures  for  obtaining  the  Type  0  and  Type  1  DFCS 
structures  using  optimal  control  are  detailed  fully  in  Refs.  3 
to  10.  The  proportional-integral  controller  structure  in  in¬ 
cremental  form  developed  in  Ref.  7  has  been  programmed  into  a 
helicopter  control  minicomputer  and  successfully  flight  tested; 
it  will  be  the  structure  used  in  this  study. 

The  Type  0  and  Type  1  control  laws  are  obtained  by 
forming  a  continuous- time  optimal  control  problem,  transform¬ 
ing  the  problem  to  discrete- time  using  the  sampled-data  regu¬ 
lator,  then  solving  the  discrete- time  problem.  The  continuous¬ 
time  linearized  model  of  the  aircraft  with  adjoined  compensa¬ 
tor  states  is  shown  in  Table  1  along  with  the  optimal  control 
problem  transformation  accomplished  by  the  computer  algorithms 
during  the  design  process.  The  designer  works  in  continuous 
time  by  specifying  the  cost  function  weights  Q  and  R  and^ com¬ 
puter  algorithms  determine  the  discrete  time  weights  Q,  M,  and 

a 

R.  The  optimal  control  problem  is  solved  using  the  Riccati 
equation  in  Table  2  and  the  Type  0  and  Type  1  gains  are  com¬ 
puted.  The  Type  0  and  Type  1  gains  are  found  from  the  same 
optimal  control  solution  and  result  in  the  same  closed-loop 
eigenvalues  and  state  eigenvectors  if  the  number  of  controls 
equals  the  number  of  commands. 

Table  2  shows  the  perturbation  Type  1  control  law 
which  operates  near  the  nominal  controls,  u  ,  and  state,  xQ . 

As  shown  in  Ref.  7,  the  Type  1  control  law  can  be  implemented 
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TABLE  2 

SOLUTION  OF  THE  DISCRETE-TIME  OPTIMAL  CONTROL  PROBLEM 


T-2903 


DISCRETE-TIME  RICCATI  EQUATION 
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TYPE  1  DISCRETE-TIME  OPTIMA!  CONTROL  UW  GAINS 
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TYPE  l  CONTROL  UW  IN  INCREMENTAL  FORM 
ERROR  DYNAMICS l 

*  Aflk.i  -  Cl(ASk  -  Aak.j)  * 

ORIGINAL  VARIABLES  I 

AHk  *  Auk,1  -  CjU£k  -  A5k.l)  *  AtC2  jA2|<-l  *  AXd,k) 


*  lA21  *  C1AU]  lAVk  ‘  AVk-l> 


without,  requiring  the  nominals  if  the  incremental  form  dis¬ 
cussed  in  Ref,  11  is  used.  The  total  value,  zero-order  hold, 
Type  1  control  law  in  incremental  form  is 

yk  =  -k-l  '  C1  ^-k  “  -k-l]  "  AtC2  ^k-1  “  ^d,k]  (1) 

In  F.q.  1,  u^  is  the  (mxl)  control  command  vector  which  is  sent 
to  the  control  actuators  at  time  t^,  x^  is  the  measured  or 


8 


I 


estimated  value  of  the  (nxl)  aircraft  state  vector  in  body 
axes  used  for  feedback  at  time  t^,  ^k-1  t^*e  measured 

or  estimated  value  of  the  output  vector  at  time  t^^  which  is 
to  match  the  pilot  commands,  and  ^  is  the  (£xl)  pilot  com¬ 
mand  vector  at  time  t^.  For  simplicity  all  states  are  assumed 
measured  and  no  estimation  techniques  are  employed.  Detailed 
estimation  techniques  are  discussed  in  Ref.  7. 

Equation  1  is  the  Type  1  control  law  simulated  in 
Chapter  3  with  the  feedback  gain  and  the  integrator  gain  C2 
scheduled  with  flight  condition.  The  prime  on  -*-n  ^9 •  ^ 

indicates  the  control  command  is  reset  always  to  lie  within 
the  control  actuator  range  during  control  saturation  and  to 
prevent  windup  as  discussed  in  Ref.  11.  In  the  nonlinear  sim¬ 
ulations  in  Chapter  3,  maneuvers  are  performed  which  require 
maximum  performance  from  the  aircraft  and  which  briefly  sat¬ 
urate  some  of  the  controls  in  u^. 

The  Type  1  control  law  eliminates  the  need  for  com¬ 
puting  nominal  values  by  taking  advantage  of  the  integral  pro¬ 
perty.  The  Type  0  control  law  does  not  have  this  property  and 
depends  on  knowledge  of  the  nominal.  The  total  value  Type  0 
control  law  is  as  follows  using  Ref.  7, 

Uk  =  Hk-i  ISk.!  -  £]*]  -“K2lHk-l  -  (2) 

*£  =  *o,k  +  n12  ^d.k  -  *o,k>  (3) 

=  Ho.k  +  S122  <*d,k  '  £o,k>  W 

The  states  and  in  Eqs.  3  and  4  are  estimates  of 
the  current  nominal  states  and  controls  of  the  aircraft  based 
on  perturations  from  some  fixed  nominal.  The  states  xQ  ,  uo  , 
and  yo  are  the  fixed  nominal  of  the  aircraft.  Periodically, 
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if  the  pilot  command,  is  significantly  different  from 

the.  fixed  nominal  command  y  .  ,  ,  then  x  .  ,  u  ,,  and  y  . 

O  j  K“  1  "  O  |  K  "“O  j  K  O  j  K. 

must  be  updated  as  the  aircraft  transitions  between  flight 
conditions.  Numerous  procedures  exist  for  determining  the 
fixed  nominal  values  in  flight,  including  trim  maps,  sched¬ 
uling,  or  updating  xQ,  uQ  and  using  low-pass  filtered  ver 
sions  of  x,  u,  and 


In  Eq .  2,  is  the  state  feedback  gain  and  the  gain 
I<2  causes  the  low-pass  filter  action  of  the  Type  0  control 
law.  The  matrices  an^  ^22  are  t^ie  stea^y  state  feedforward 
gains  (Ref.  12)  and  are  obtained  from  the  following  equations 


<*-I> 

H 


Ax* 

0 

Au* 

J 

_A*d_ 
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*1 

R 

r  "1 

<*-n 

r 

nll 

fi12 

H 

r 

D 

r 

J*21 

°22_ 

Ax* 

Au* 


fi12  A^d 
n22  A^d 


(5) 


(6) 

(7) 


The  matrices  <t>  and  I"  are  the  discrete  equivalents  of  the  lin¬ 
earized  aircraft  model  matrices,  F  and  G,  as  shown  in  Table  1. 
The  matrices  H  and  D  transform  the  aircraft  states  and  controls 
to  the  commanded  output,  A^,  as  follows 


A^  =  HAx  +  DAu 

The  vectors  Ax*  and  Au*  are  the  nominal  perturbations  which 
must  be  added  to  current  nominal  states  and  controls  of  the 
aircraft,  xq  and  uQ  when  a  perturbation  in  the  command  vector, 
Ay^ ,  is  requested  by  the  pilot. 


±2 
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In  Eq.  6,  superscript  R  denotes  the  right  inverse. 

The  right  inverse  solves  Eq.  5  when  the  number  of  commands  ,  £, 
is  less  than  the  number  of  controls  m.  The  right  inverse  used 
in  this  report  is  computed  using  the  weighted  pseudoinverse, 


r  — 

R 

r  n 

(«*>-! )  r 

(4>-i)  r 

=  Q"1 

H  D_ 

a 

X 

_ 1 

Singular  value  decomposition  is  used  in  the  computations  in 
Eq.  8  for  improved  numerical  accuracy.  Section  2.2  discusses 
values  chosen  for  the  pseudoinverse  weighting  matrix,  Q  used 
in  Eq.  8. 


The  Type  0  gains  are  related  to  the  Type  1  gains  as 
shown  in  Table  2.  Kj  and  are  determined  from  a  Riccati 
equation  solution  and  C-^  and  C2  are  computed  as  follows: 


[Cl  AtC2  =  AtK| 


<4*-I) 

H 


(9) 


The  values  for  Qg  used  to  obtain  n^2  anc*  fi22  not  ^ave  to 
the  same  as  the  values  in  Qg  used  to  obtain  and  C2  in  Eq . 

9.  The  design  in  this  report  uses  different  values  for  Q  in 

s 

Eq .  6  and  Eq.  9  for  reasons  discussed  in  the  following  sections. 


There  are  usually  an  infinite  number  of  solutions  to 
Eq.  5  when  £<m;  the  pseudoinverse  shown  in  Eq .  9  represents 
the  weighted  least  squares  solution.  If  Eq.  7  is  substituted 
into  Eq.  5  the  equality  is  true  but  if  Eq .  9  is  substituted 
for  and  C2  in  the  equation  in  Table  2,  equality  usually 
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does  not  hold  for  £<m.  If  the  number  of  controls  equals  the 
number  of  commands  then  equality  holds  in  Eq.  5  and  in  Table 
2.  In  the  aircraft  control  design  in  this  report,  the  number 
of  commands,  £ ,  is  less  than  the  number  of  controls,  m.  Ref¬ 
erences  2  and  3  are  earlier  continuous-time  discussions  of  the 
command/con t.rol  requirements  for  the  digital  ACM  controller 
used  here  and  these  references  describe  reasons  why  treating 
all  aircraft  controls  independently  is  a  desirable  feature. 
These  reasons  are  as  follows: 


•  By  varying  the  weights  in  Q  in  Eq.  8 

it  is  possible  to  decrease  control  steady- 
state  positions  while  increasing  aircraft 
final  state  values.  The  result  is  to  in¬ 
crease  the  range  of  the  aircraft  states 
that  can  be  commanded  as  discussed  in 
Section  2.2. 

•  For  maximum  performance  maneuvers  it  is 
possible  to  favorably  saturate  one  con¬ 
trol  while  leaving  the  other  aircraft 
controls  operational.  The  result  is  to 
increase  performance  without  causing 
aircraft  departures  and  is  demonstrated 
in  Chapter  3. 

•  When  the  number  of  commands  and  controls 
are  equal,  the  aircraft  develops  a  finite 
set  of  transmission  zeroes  (Ref.  13) 
which  play  an  important  part  in  tracking 
theory  (Refs.  14  and  12).  These  trans¬ 
mission  zeroes  cannot  be  changed  by  state 
feedback.  Degraded  handling  qualities 
occur  if  any  transmission  zeroes  are  in 
the  right-half  complex  plane  (i.e.,  non¬ 
minimum  phase).  When  the  number  of  con¬ 
trols  is  greater  than  the  number  of  com¬ 
mands,  the  plant  still  has  zeroes  but 
they  can  be  changed  by  adjusting  Q  in 

Eq ,  8.  In  particular,  right-half  plane 
zeroes  can  be  favorably  moved  into  the 
left-half  (minimum  phase)  complex  plane 
as  discussed  in  Section  2.4. 
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•  If  the  number  of  commands  is  less  than 
number  of  the  controls,  the  number  of 
closed-loop  modes  that  must  be  stabil¬ 
ized  is  reduced  when  converting  from  the 
Type  0  to  the  Type  1  control  law.  Extra 
design  freedom  exists  because  some  of 
the  closed-loop  modes  shift  during  the 
transformation  if  Q  in  Eq.  8  is  changed. 
If  Q  is  chosen  appropriately  the  shifts 
are  all  favorable  as  discussed  in  Section 
2.4. 


The  Type  0  controller  Shown  in  Eqs.  2  through  4  is 
considerably  more  complex  to  implement  than  the  Type  1  con¬ 
troller  in  Eq.  1,  Other  disadvantages  of  the  Type  0  controller 
are  discussed  in  the  following  sections. 


2.2  PSEUDOINVERSE  WEIGHTING  MATRIX  FOR  FIGHTER 
AIRCRAFT  APPLICATIONS 


This  section  presents  the  design  of  the  pseudoinverse 
weighting  matrix,  Q_ ,  used  in  Eq.  8  for  specification  of  the 

D 

Type  0  DFCS  feedforward  matrices.  Unity  weighting  is  used  for 
the  continuous-time  Type  0  design  in  Ref.  3  and  it  will  be 
shown  that  this  is  not  necessarily  the  best  choice. 


The  states,  controls,  and  commands  for  the  aircraft 

are 


Ax'*'  =  jjiq  Aw  Av  Ar 
4HT  =  [«s  «sp  Sr] 

A^d  *  [4an  4P  4Pw] 

Three  inputs,  normal  acceleration,  Aa  ,  sideslip,  Ap ,  and  wind 
axis  roll  rate,  Ap^,  are  to  command  five  controls:  stabilator, 


A6g,  main  flaps,  A6m£,  spoilers,  A6gp,  differential  stabila- 
tors,  46^  and  rudder,  A6r  The  DFCS  design  is  nearly  the 
same  as  in  Ref.  3,  except  it  is  to  be  done  in  discrete  time, 
and,  for  reasons  discussed  in  the  next  paragraph,  Au  is  not 
included  in  the  state  vector. 

The  state  Au  is  dropped  from  the  control  design  state 
vector  because  the  discrete  and  continuous  time  Type  0  DFCS 
designs  yield  a  very  poor  Aan  step  response  with  Au  included. 
The  Type  0 . DFCS  at  design  point  2  in  Ref.  3,  Table  14,  has  a 
closed-loop  eigenvalue  with  a  time  constant  of  19.3  sec.  This 
slow  time  constant  disappears  from  the  Type  1  eigenvalues  as 
shown  in  Ref.  3.  The  slow  closed-loop  eigenvalue  is  one  half 
of  a  phugoid  complex  pair  and  the  eigenvector  has  a  large  Au 
contribution  (Ref.  3,  Table  15).  The  slow  eigenvalue  is  re¬ 
moved  from  the  Type  1  eigenvalues  by  using  very  large  gains  on 
Au  (Ref.  3,  Figs.  56  and  58).  When  the  Type  0  Aan  step  re¬ 
sponse  is  simulated,  there  is  a  contribution  to  Aan  from  the 
change  in  Au  in  the  HAx  +  DAu  observation  vector.  Since  the 
Type  0  Au  response  has  a  19.3  sec  time  constant,  the  Type  0 
Aa_  response  also  contains  the  slow  response.  It  is  intriguing 
that  the  Type  0  to  Type  1  transformation  automatically  com¬ 
pensated  for  the  poor  Au  response  and  indicated  superior  Aan 
command  simulations.  The  Type  0  Aa^  response  was  not  attempt¬ 
ed  in  Ref.  3  and  the  Au  problem  was  not  identified  until  the 
digital  Type  0  DFCS  was  simulated.  When  Au  is  removed  from 
the  state  vector,  the  speed  response  of  the  aircraft  becomes 
open  loop  subject  to  any  throttle  activity  of  the  pilot,  (In 
ACM  the  pilot  usually  uses  full  throttle.)  Within  the  time 
frame  of  interest  for  a  Aan  step  response  the  open-loop,  Au 
response  is  close  to  zero  and  does  not  have  a  significant  ef¬ 
fect  on  performance. 
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The  purpose  of  the  pseudoinverse  weighting  matrix, 

Q  ,  in  Eq.  8  is  to  trade  off  steady  state  values  of  the  air- 
craft,  states  against  steady  state  values  of  aircraft  control 
position  needed  to  satisfy  the  command  requirements.  If  more 
weight  is  placed  on  the  controls  in  Q  ,  the  effect  is  to  in- 
crease  the  values  of  the  states  that  are  not  commanded  by 
while  at  the  same  time  decreasing  the  values  of  the  controls 
needed  to  meet  the  command  requirements.  The  net  result  is  to 
increase  the  effective  operating  range  of  the  aircraft  states 
before  control  saturation  is  encountered,  a  very  desirable 
situation . 


A  pseudoinverse  v?eighting  matrix  which  has  the  desired 
weighting  pattern  uses  unity  weighting  for  velocity  states,  v 
and  w,  in  fps,  unity  weighting  for  angular  rates,  p,  q,  and  r 
in  degrees,  and  a  multiple  of  the  inverse  of  the  total  amount 
of  individual  control  travel  for  the  control  position  weights. 


1 


0 


Qc  = 


n 


aSJ 


max 


(10) 


0 


max 


The  control  weighting  pattern  causes  controls  with  the  least 
amount  of  travel  to  have  the  most  weight.  Table  3  shows  the 
changes  in  the  elements  in  anc*  ^12  cause<^  by  changing  the 
value  of  a  in  Eq .  10  between  1  and  2000, 
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TABLE  3 

EFFECTS  OF  CHANGING  THE  PSEUDO INVERSE  WEIGHTING  MATRIX 
ON  QUAD  PARTITION  INVERSE  ELEMENTS 


T-1794 


fl22 

or 

n’2 

ELEMENT 

UNITS 

-  . 

INVERSE* 

PSEUDOINVERSE-CONTROL  WEIGHTS 

UNITY 

2 

ua,; 

UNITY 

20 

mra 

200 

mot 

2000 

Max 

AVA«nc 

...  - *- 

deR/fps 

-0.11B 

-0.378 

-0,321 

-0.263 

-0.119 

-0.0650 

A«mf/AanC 

deg/ fps2 

- 

0,703 

0,557 

0.360 

0.005 

-0,141 

dog/deg 

- 

-1.06 

-0.226 

-0.962 

-0.962 

-0,061 

A*dB/4Bc 

deg/deg 

•  1,22 

-1.04 

-1.18 

-1.05 

-1.05 

-1,05 

AAr/ABc 

deg/deg 

1.11 

1.10 

1.11 

1.10 

1.10 

1.10 

Afisp/APWc 

deg/deg/sec 

- 

-0,117 

-0.0268 

-0.116 

-0. 116 

-0.116 

ASds/APwc 

deg/deg/ sec 

-0. 150 

-0. 120 

-0.145 

-0.120 

-0.129 

-0.120 

AyApwc 

deg/deg/sec 

-0.00174 

-0.00330 

-0.00200 

-0.00338 

-0.00338 

-U. 00330 

Aq/A«nc 

deg/deg/ fpe2 

0.0076 

0.0083 

0 , 0080 

0.0969 

0.0076 

0.0979 

A«/A«t)C 

fps/.fpfl2 

1.40 

0.0135 

0.321 

0.714 

1,48 

1.70 

flv^APwc 

fps/dng/sec 

C.O 

0.0 

0.0 

0.0 

0.0 

0.0 

Ar/APwc 

deg/uee/deg/sec 

0.109 

0,160 

0.160 

0.160 

0.189 

0.169 

Ap/AVwc 

deg/sec/deg/soc 

O.  986 

0.086 

0.986 

0.086 

0.986 

0.986 

Av/ABC 

fpn/deg 

10.4  7 

10.47 

10.47 

10.47 

10.47 

10.47 

Ar/ABC 

deK/seo/di"g 

-0.0684 

-0.0664 

-0.0670 

-0.0666 

-0 . 0660 

-0.0068 

Ap/ABC 

deg/sec/deg 

0.100 

0. 100 

0.  100 

0.100 

0.100 

0.0984 

and  are  removed  from  the  controls  to  construct  the  Inverse, 

mf  sp 


The  flight  condition  used  in  Table  3  is  trimmed  straight 
flight  with  a  forward  velocity,  pitch  rate  and  angle  of  attack 
of  183  m/s  (600  fps),  5  deg/sec,  and  10  deg,  respectively, 

(design  point  1  in  Ref.  2).  Three  important  observations  can 
be  made:  1)  lateral  states  and  control  positions  change  little 
with  increasing  pseudoinverse  control  weights,  2)  the  longitu¬ 
dinal  states  trade  off  Aw  with  control  position,  ..nd  3)  con¬ 
siderable  deviations  in  the  maneuver  flap  steady  state  value 
can  occur  even  through  maneuver  flap  has  the  least  amount  of 
control  travel  before  saturation  (10  deg). 
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The  above  observations  lead  to  the  conclusion  that 
the  weighting  matrix  should  be  chosen  based  on  the  relation¬ 
ship  between  stabilator  and  maneuver  flap.  Figure  2  shows 


A6s/Aanc  and  46mf/Aanc  variations  at  183  m/s  (600  fps)  with 


has  a  small  contribution  in  aiding  A6  , 

s 


A6mf  significantly  increases  the  effec- 


increasing  angle  of  attack  and  pitch  rate  flight  conditions. 
Below  20  deg  0'o,  A6f 
but  beyond  20  deg  a( 

tive  operating  range  of  the  aircraft  to  higher  angles  of  at¬ 
tack  before  A<5g  saturates.  The  weighting  parameter  100/MAX 
(2.5,  10,  1.85,  4.16,  1.66  along  the  control  diagonal  in  Qc) 
produces  the  desired  result.  It  provides  a  ratio  below  \  for 
A6mf/45 s  when  aQ  is  less  than  20  deg.  Above  27  deg  a0  100/MAX 
weighting  causes  A6m£  to  saturate  favorably  in  order  to  satisfy 


an  increase  of  9.8  m/s  (lg)  Aan  command. 


2.3  SYNTHESIS  OF  CONTROLLERS  USING  CONTINUOUS- 
TIME  SPECIFICATIONS 


This  section  presents  the  synthesis  of  the  Type  0  and 
Type  1  digital  control  law  design  at  one  flight  condition. 

One  design  approach  for  digital  flight  control  systems  is  to 
design  an  analog  system  and  then  discretize  it.  The  result 
often  leads  to  a  high  sample  rate  system  which  places  signif¬ 
icant  requirements  on  the  flight  computer  used  to  implement 
the  design.  One  reason  this  practice  continues  is  that  a  large 
body  of  information  and  experience  exists  for  continuous-time 
designs  including  MIL-F-8785B  (Ref.  15).  The  continuous- time  * 
design  in  Ref.  3  could  have  been  emulated  digitally  to  produce 
a  discrete-time  controller. 

Examples  of  analog  designs  implemented  digitally  us¬ 
ing  the  Tustin  approximation  are  shown  in  Ref.  16  for  the  space 
shuttle  (25  Hz  and  12.5  Hz  sampling  time)  and  in  Ref.  17  for  a 
fighter  aircraft  (80  Hz  and  40  Hz  sampling  time).  It  is  clear, 
however,  that  improved  performance  and  reduced  computer  re¬ 
quirements  can  be  obtained  using  a  direct  digital  design  at 
lower  sampling  rates  which  take  better  advantage  of  digital 
system  characteristics. 

The  approach  used  in  this  section  to  obtain  a  direct 
digital  design  which  at  the  same  time  retains  continuous- time 
features  is  the  sampled-data  regulator  discussed  in  Ref.  18. 

The  sampled-data  regulator  transforms  a  continuous-time  system 
and  cost  function  to  their  equivalent  discrete-time  represen¬ 
tation  as  shown  in  Table  1.  The  discrete-time  design  gives  a 
stable  discrete-time  closed-loop  system  with  eigenvalues  inside 
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the  unit  circle,  This  section  shows  that  transforming  the  closed- 
loop  discrete-time  system  back  to  continuous- time  gives  an  equi¬ 
valent  closed-loop  continuous- time  system  that  can  be  used  to 
evaluate  the  digital  design  using  continuous- time  specifica¬ 
tions,  Also  shown  are  the  cost  function  weighting  elements 
and  the  flight  conditions  used  for  point  design.  In  point 
design,  which  is  the  basis  for  all  the  designs  discussed  in 
Ref.  1  to  7 ,  a  number  of  linear  time-invariant  aircraft  models 
are  determined  at  a  set  of  trimmed  flight  conditions  which  span 
the  flight  regime  of  interest.  The  control  law  gains  are  found 
at  each  flight  condition.  The  gains  are  regressed  against  a  set 
of  nonlinear  functions  of  the  major  flight  variables  and  the  re¬ 
sulting  gain  schedules  are  suitable  for  onboard  implementation. 

The  flight  conditions  used  for  point  design  span  a 
large  range  of  trimmed  values  of  velocity,  angle  of  attack, 
and  roll  rate  about  the  velocity  vector.  Thirty  flight  con¬ 
ditions  are  used  with  velocity  either  122,  183,  or  244  m/s, 
wind-axis  roll  rate  either  0,  50,  or  100  deg/sec,  and  angle  of 
attack  varying  between  0  and  33  deg.  Nonzero  sideslip  condi¬ 
tions  are  not  included  since  sideslip  is  usually  commanded  to 
be  zero.  Figure  3  shows  the  flight  conditions  in  a  three  di¬ 
mensional  configuration. 

The  sampling  rate  for  the  DFCS  is  chosen  to  be  20 
samples/sec.  Results  in  Ref.  19  and  20  indicate  that  a  direct 
digital  design  can  go  as  low  as  10  samples/sec  and  still  pro¬ 
vide  adequate  performance  with  unmodeled  body  bending  modes 
and  turbulence  in  the  loop.  The  higher  samplirg  rate  in  the 
DFCS  serves  to  reduce  pilot  awareness  of  the  control  surface 
steps  and  to  keep  the  control  surface  phase  lag  from  exceeding 
30  deg  (Ref.  15). 


4  -U417 


122  183  244 

(400)  (800)  (800) 


VELOCITY,  (m/i) 

Figure  3  Flight  Conditions  for  DFCS  Point  Design 

At  each  flight  condition,  closed-loop  eigenvalues  and 
step-response  simulation  (i.e,,  the  eigenvectors)  of  the  lin¬ 
earized  system  are  used  to  evaluate  the  design,  Closed-loop 
eigenvalues  should  meet  the  requirements  in  MIL-F-8785B  for 
Category  A,  Level  1  flight  of  a  Class  IV  aircraft.  Step  re¬ 
sponse  transients  should  fall  within  envelopes  which  are  known 

i 

to  generate  good  pilot  ratings  as  shown  in  Refs.  10  and  17. 

The  closed-loop  system  for  the  Type  0  DFCS  is 


(ID 
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while  the  closed-loop  system  for  the  Type  1  DFCS  in  position 
form  using  perturbation  variables  and  an  integrator  state,  A£ , 
is 


Ax  4>"rcl 

A£  AtH-AtDC^ 

k+1 


If  the  number  of  commands  and  controls  are  equal,  the  eigen¬ 
values  of  the  two  closed-loop  systems  are  equal  and  the  feed¬ 
forward  gains  are  unique.  If  the  number  of  commands  is  less 
than  the  number  of  controls  <the  situation  for  the  design  here) 
then  the  eigenvalues  between  the  two  differently  dimensioned 
closed-loop-systems  cannot  be  all  equal  but  are  closely  aligned 
depending  on  the  weighting  matrix  in  the  pseudoinverse  identi¬ 
fied  in  Eq.  9.  For  reasons  which  are  discussed  in  the  next 
section,  unity  weighting  (Q  =  I)  is  used  in  Eq.  9. 

The  closed-loop  eigenvalues  and  eigenvectors  are  ad¬ 
justed  by  iterating  on  Q  and  R,  the  quadratic  cost  function 
weights.  Although  emulating  the  continuous- time  design  in 
Ref.  3  is  considered  a  poor  practice,  the  continuous-time  cost 
function  weighting  elements  can  be  used  directly  in  the  sampled- 
data  regulator  digital  design.  Continued  refinement  of  per¬ 
formance,  however,  caused  the  choice  of  Q  and  R  to  differ  from 
those  in  Ref.  3.  The  change  in  the  cost  function  weighting 
elements  is  shown  in  Table  A.  The  primary  changes  are  a  de¬ 
crease  in  Aw,  Aan  and  Av  weighting,  and  an  increase  in  Apw  and 
control  position  weighting.  The  effect  on  the  closed-loop 
system  is  to  decrease  the  short  period  natural  frequency  and 
the  stability-axis  roll  rate  rise  time,  while  increasing  the 
normal  acceleration  and  sideslip  rise  times.  These  effects 
improve  the  command  performance. 
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TABLE  4 

DFCS  WEIGHTS  AT  DESIGN  POINT  I 


J  ■ 


T-3117 


MATRIX 

MATRIX 

TYPE 

MATRIX  ELEMENT 

MAXIMUM  MEAN 
VALUE  IN  REF.  2 

MAXIMUM  MEAN 
VALUE  HERE 

Q 

S  t.  fl  t «! 
Position 

Axial  Velocity,  Au 

12.2  m/s 
(AO  fps) 

Not  used 

Lateral  Velocity,  Av 

3.05  m/s 
(10  fps) 

6.1  in/s 
(20  fps) 

Normal  Velocity,  Aw 

3 . 66  m/s 
(12  fps) 

7.6  m/s 
(25  fps) 

Body  Angular  Kates 

20  deg/sec 

20  deg/sec 

Q 

State 

Rate 

Lateral  Acceleration,  Av 

3 . 66  m/s? 

(12  fps* ) 

3.0  m/s? 

(10  fps2) 

Normal  Acceleration,  aw 

1 . 53  m/s? 

(5  f pR  * 

6.1  m/s? 

(20  fps*) 

Q 

Pilot 

Command- 

Stability-Axis  Normal 

Acceleration  Command,  Aa_ 

n 

0.533  m/s? 
(1.75  fps* 

0.92  m/s? 

(3.0  fps2 

Sideslip 

Sideslip  Command,  ps 

0.9  deg 

1 , 0  deg 

Stability-Axis  Roll 

Rate  Command,  Apw 

2.5  deg/sec 

2.0  deg/sec 

Q 

Control 
Posi tion 

Stabilstor  Deflection,  AAp 

10  deg 

5  deg 

Maneuver  Flap  Deflection, 

Afi  r 

mf 

5  deg 

1  . 25  deg 

Spoiler  Deflection,  A6dR 

27  deg 

3,375  deg 

Differentia]  Stabilator 
Deflection,  Afi^ 

6  deg 

3 , 0  deg 

Rudder  Deflection,  A0r 

15  deg 

7.5  deg 

R 

Control 

Rate 

Stabilator  Hate,  A6g 

3  deg/sec 

3  deg/sec 

Maneuver  Flap  Rate,  A6  p 

A  deg/sec 

A  deg/sec 

Spoiler  R.ite,  % 

HP 

A  deg/sec 

A  deg/sec 

Differential  Stabilator 
Rate,  A*df) 

3  deg/sec 

3  deg/sec 

Rudder  Rate,  A^r 

A  deg/sec 

A  deg/sec 

V, 


The  closed-loop  eigenvalues  for  the  DFCS  using  Eqs. 

11  and  12  all  fall  within  the  unit  circle.  It  is  reasonable 
to  presume  that  if  the  continuous  plant  is  converted  to  an 
equivalent  discrete  system  using  a  matrix  exponential  in  Table 
1,  then  a  discrete  system  can  be  converted  to  an  equivalent 
continuous  system  using  the  natural  logarithm  of  a  matrix.  If 
4>cl  is  either  of  the  closed-loop  system  matrices  in  Eqs.  11 
and  12,  then 


(“cl-1) 


2 


'"“CL-1* 


3 


is  the  equivalent  continuous  system.  The  sequence  for  the 
logarithm  of  a  matrix  converges  if  the  eigenvalues  of  «t>c^  are 
in  the  open  right  half  complex  plane  and  the  magnitude  of  the 
eigenvalues  are  less  than  or  equal  to  two. 

At  the  nominal  design  flight  conditions  (V  »  183  m/s 
(600  fpt,),  a  =  9  deg,  pw  =  0),  Table  5  shows  a  comparison  between 
eigenvalues  of  the  different  closed-loop  systems  in  Eqs.  11 
and  12  for  the  same  Q  and  R.  The  continuous- time  closed- loop 
eigenvalues  determined  by  minimizing  the  continuous-time  cost 
function  in  Table  1  and  the  mapped  discrete-time  eigenvalues 
are  almost  identical.  Furthermore,  changing  the  sampling  time 
to  a  larger  but  not  unreasonable  sampling  time  has  little  ef¬ 
fect  on  the  mapped  eigenvalue  locations,  particularly  in  the 
Type  0  system  which  does  not  involve  using  the  pseudoinverse 
to  compute  the  feedback  gains.  A  comparison  between  Dutch 
roll  mode  eigenvectors  in  Figure  4  shows  the  eigenvectors  are 
also  almost  identical.  It  is  concluded  that  evaluating  the 
digital  design  using  the  mapped  eigenvalues  is  a  useful  pro¬ 
cedure  for  evaluating  a  digital  system's  performance  using 
continuous- time  specifications. 


23 


TABLE  5 

COMPARISON  BETWEEN  CONTINUOUS-TIME  AND  MAPPED 
DISCRETE-TIME  CLOSED-LOOP  SYSTEM  EIGENVALUES 


CLOSED-LOOP  CHARACTERISTICS 
CONTINUOUS  DISCRETE  TIME  TYPE  0 

— UN  ITV '  MEUBBTOgg  TFfCHTlro - 

CLOSED-LOOP  CHARACTERISTICS 
CONTINUOUS  DISCRETE  TIME  TYPE  1 

DYNAMIC  MODE 

Natural 
Frequency , 
Rad/Sec 

Damping 

Ratio 

Time 

Cone  tent. 
See 

Natural 
Frequency, 
Rad/Sec . 

Damping 

Ratio 

Tima 

Conatant 
sac  . . 

Short  Period 

4.19  ' 

0.662 

5.32 

0.735 

Roll  Command 

4,34 

0.686 

— 

4.34 

0.686 

•  ■ 

Dutch  Roll 

3.46 

0.658 

-- 

3.46 

0.658 

•  ** 

Normal  Acceleration 
Command 

-■ 

-- 

0,320 

“  - 

0,363 

Longitudinal  Control 

— 

- 

0,400 

— 

•» 

Lateral  Control 

-- 

0.849 

-- 

-- 

Sideslip  Command 

— 

- 

0.975 

-- 

-- 

0.974 

SAMPLING  TIME  *  0,05  aac 

CLOSED- LOOP  CHARACTERISTICS 
DISCRETE  TIME  TYPE  0 

iMW  psH 

CLOSED-LC 

.DISCR 

OP 1 CHARACTER 1 ST I CS  ° 

,TE  TIME  TYPE  .1 _ 

DYNAMIC  MODE 

Natural 
Frequency , 
Rad/Sec  . 

Damping 
Ratio . 

Time 

Conatant. 

Sec 

Natural 
Frequency , 
Rad/Sec. 

Damping 
Ratio . 

Time 

Conatant 

Sep  .  .. 

Short  Period 

4,21 

0.665 

5.54 

0.715 

Roll  Command 

4.34 

0.682 

*• 

4.34 

0.682 

•  — 

Dutch  Roll 

3.46 

0.656 

-- 

3.46 

0.655 

*  * 

Normal  Acceleration 
Command 

-- 

0.322 

■  • 

0.36B 

Longitudinal  Control 

- 

-- 

0.400 

at  • 

.. 

-- 

Lateral  Control 

- 

-- 

0,849 

-• 

-- 

-- 

Sidealip  Command 

— 

— 

0.978 

-■ 

-- 

0,977 

SAMPLING  TIME  »  0,1  aac. 

CLOSED- LOOP  CHARACTERISTICS 
DISCRETE  TIME  TYPE  C 

UNITY  PSEUDO IN vERsEwErGHTING 
CLOSED-LOOP  CHARACTERISTICS 
_ DISCRETE  TIME  TYPE  1 _ 

DYNAMIC  MODE 

Natural 
Frequency , 
Rad/Sec.  . 

Damping 
.  R*ti.o._  .. 

Time 

Conatant . 
Sec 

Natural 

Frequency, 

-  Raa/Stc  _ 

Damping 

Jilla _ 

Time 

Constant 

Sec 

Short  Period 

0.24 

0.666 

5.75 

0.684 

-- 

Roll  Command 

4.34 

0,671 

-- 

4.34 

0,671 

-- 

Dutch  Roll 

3,46 

0,653 

3.46 

0,653 

-■ 

Normal  Acceleration 
Command 

-- 

0.324 

-  - 

"  " 

0.371 

Longitudinal  Control 

-- 

-• 

0.402 

-- 

-- 

— 

Lateral  Control 

-■ 

-- 

0.848 

-- 

-- 

-- 

Sideslip  Command 

-- 

-- 

0.984 

-- 

0.983 

24 


R-42436 


£11 

n 

Z§0 

o 

□ 

u 

AS 

»P 

A<5r 

CLOSED-LOOP  CONTINUOUS  TIME  TYPE  0 


/-I  22 
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CLOSED-LOOP  DISCRETE-TIME  0,  At  -  0.1  SEC 

Figure  4  Comparison  Between  Continuous-Time  and 
Mapped  Discrete-Time  Closed-Loop  System 
Eigenvectors  for  the  Dutch  Roll  Mode 


Using  MIL-F-8785B  the  desirable  range  of  short  period 
natural  mode  frequency  at  the  nominal  flight  conditions  under 
discussion  is  0.9  rad/sec  to  6.8  rad/sec  with  a  damping  ratio 
between  0.35  and  1.3.  Minimum  Dutch  roll  natural  frequency, 
uj  ,  and  damping  ratio,  requirements  are  0.19  for  0.35 

fo$  ,  and  1.0  rad/sec  for  uj  .  The  roll  mode  time  con- 

Q  Hi  n j 

stant  should  be  less  than  1.0  sec?  All  of  these  requirements 
are  met  by  the  mapped  eigenvalues  of  the  DFCS  at  the  nominal 
flight  condition  as  shown  in  Table  5. 


2.4  COMPARISON  OF  DISCRETE  DFCS  TYPE  0  AND  TYPE  1  RESPONSE 

Type  0  and  Type  1  discrete  DFCS  design  step  responses 
are  presented  in  this  section.  The  flight  condition  for  the 
reference  aircraft  is  the  same  as  the  one  used  for  eigenvalue 
analysis  in  Section  2.3.  The  discrete  state  time  histories 
are  interpolated  to  provide  the  approximate  continuous-time 
simulation.  The  discrete  control  time  histories  are  stair¬ 
cased  in  the  figures  and  indicate  the  step  commands  the  ref¬ 
erence  aircraft  actuators  receive  from  the  control  law.  The 
control  laws  are  simulated  using  the  incremental  form  shown  in 
Eqs.  1  and  2. 

Experiments  with  large  pseudoinverse  weighting  ele¬ 
ments  for  the  states  in  going  from  the  Type  0  to  the  Type  1 
DFCS  have  shown  significant  and  in  some  cases  unstable  eigen¬ 
values  changes.  Care  must  be  taken  in  using  non-unity  pseudo¬ 
inverse  weighting  matrices  when  transforming  from  the  Type  0 
to  the  Type  1  control  design.  Keep  in  mind  that  the  pseudo¬ 
inverse  used  in  Eq.  6  to  find  2  anc*  ^22  anc*  t*)e  Pseu<^°^n“ 
verse  used  in  Eq .  9  to  find  C^  and  C ^  are  similar  but  have 
completely  different  effects  on  the  control  law.  The  former 
affects  steady  state  conditions  while  the  latter  affects 
closed- loop  eigenvalues. 
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To  compare  the  effect  of  eigenvalue  and  step  response 
changes  caused  by  the  pseudoinverse  weighting  matrix,  three 
cases  are  considered.  In  all  three  cases  the  weighting  on  the 
states  is  unity  while  the  control  weightings  are  unity,  100/MAX 
(2.5,  10.0,  1.85,  4.16,  1.66),  and  1000/MAX  (25.0,  100.0,  18.5, 
41.6,  16,6).  The  resulting  steady  state  control  positions  are 


46* 

r -0 . 321  0.0  0.0 

i6£f 

0.557  0.0  0.0 

" Aanc  ' 

4{Sp 

= 

0.0  -0.226  -0.0258 

4(SC 

AS5s 

0.0  -1.18  -0.145 

Apwc 

46* 

0.0  1.11  -0.00209 

for  unity  weighting, 
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46* 

-0.127  0.0  0.0 

A6£f 

0.0272  0.0  0.0 

4i*T,C 

A62p 

= 

0.0  -0.960  -0.116 

4PC 

“Ss 

0.0  -1.052  -0.129 

Ap 

*wc 

1 

<o 

<J 

_ 1 

0.0  1.10  -0.00338 

100/MAX 

control  weighting  and 

'as*  " 

-0,117  0.0  0.0 

A6mf 

0.00236  0.0  0.0 

~Aanc  "" 

= 

0.0  -0.960  -0.116 

APC 

A62s 

0.0  -1.052  -0.129 

_APwc  . 

A6* 

0.0  1.10  -0.00338 

_  *-■ 

(14) 


(15) 


for  1000/WAX  control  weighting.  The  primary  effect  is  to  de¬ 
crease  the  main  flaps  steady  state  position  and  to  increase 
the  spoiler  steady  state  position. 

The  effect  of  increasing  the  control  position  weight¬ 
ing  in  the  pseudoinverse  matrix  on  the  closed-loop  Type  1  mapped 
discrete  eigenvalues  is  shown  in  Table  6.  The  only  signif¬ 
icant.  effect  occurs  in  the  longitudinal  dynamics  where  the 
short  period  natural  frequency  and  damping  ratio  decrease  to 
smaller  but  acceptable  values. 

TABLE  6 

COMPARISON  OF  TYPE  1  DISCRETE  CLOSED-LOOP 
EIGENVALUES  AT  DESIGN  POINT  1  FOR  DIFFERENT 
PSEUDOINVERSE  WEIGHTING  MATRICES 


um 


PSEUDOINVERSE 

WE I OUT 1  NO  MATRIX 

SHORT 

PERIOD 

DUTCH 

ROU. 

ROLL 

COMMAND 

NORMAL  ACCELER¬ 
ATION  COMMAND 

SIDESLIP 

COMMAND 

. . . 

Unity  for  State* 
und  Control* 

Natural  Frequency, 
Had/Sea 

8. 64 

3.46 

4.34 

— 

— 

Damping  Ratio 

0,718 

0.668 

0.682 

— 

— 

Time  Constant,  Sec 

— 

— 

— 

0.368 

0.977 

Unity  for  Staten 
loo/MAX  for  Controls 

Natural  Frnquancy, 
Rud/Sec 

4 .80 

3.46 

4.34 

— 

— 

Damping  Ratio 

0,801 

0,888 

0,683 

— 

— 

Time  Constant,  Sec 

— 

— 

— 

0.384 

0,977 

Unity  for  States 
1000/MAX  for  Controls 

Natural  Frequency, 
Rnd/Sec 

4.63 

3,40 

4.34 

— 

Damping  Ratio 

0.888 

0.688 

0.683 

— 

— 

Tlmo  Constant,  Sec 

— 

— 

0.386 

0.977 
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Since  changing  the  pseudoinverse  weighting  matrix 
primarily  affects  only  the  longitudinal  dynamics,  further  in¬ 
vestigations  will  concentrate  on  Aan  step  responses  shown  in 
Fig.  5.  Figure  6  shows  the  corresponding  control  time  his¬ 
tories.  The  interesting  results  are  the  Type  1  control  steady 
state  positions  in  Fig.  6  and  the  initial  negative  movement  of 
Aan  in  Fig.  5.  Unlike  the  Type  0  control  steady  state  posi¬ 
tion,  the  Type  1  control  steady  state  positions  do  not  agree 
with  Eqs.  14  and  15.  The  Type  0  control  law  compares  the  states 
and  controls  with  Q^Ay^  an<^  ^22^d*  respectively .  Since 
and  ^22  are  exactly  "matched"  to  the  plant,  the  Type  0  control 
laws  drives  the  command  errors,  Ax  -  Ay^  and  Au  -  ^22A^d 
to  zero.  On  the  other  hand,  the  Type, 1  control  law  only  com¬ 
pares  Ay  to  Ay^  and  drives  the  command  error  Ay  -  Ay^  to  zero, 

The  Type  1  control  law  does  not  have  a  reference  that 

would  enable  it  to  choose  among  the  number  of  control  steady 

state  positions  for  A6„  and  A6„r  which  drives  Aart  to  Aa,. 
r  s  mt  n  nc 

The  steady  state  values  for  two  controls  accommodating  only 
one  Aa^  command  error  simultaneously,  depend  on  elements  in 
and  C2  and  on  the  past  trajectory  of  the  plant  states  and  con¬ 
trols.  The  ambiguity  of  the  Type  1  control  steady  state  posi¬ 
tions  for  Aan  is,  of  course,  completely  eliminated  when  the 
additional  control,  A6m£,  is  fixed. 

When  stabilators  are  deflected  to  generate  an  upward 
normal  acceleration,  the  initial  pitching  of  the  aircraft  can 
cause  the  normal  acceleration  response  to  dip  before  arriving 
at  the  proper  value  as  shown  in  Fig,  5.  Dynamically,  the  equa¬ 
tions  of  motion  have  a  zero  in  the  right  half  plane.  By  suit¬ 
ably  adjusting  Qg  in  Eq.  8  and  by  changing  the  control  response 
it  is  possible  to  shift  the  zero  because  of  the  extra  design 
freedom  caused  by  maneuver  flap  (Ref.  13).  Figure  7  shows  the 
more  favorable  normal  acceleration  response  caused  by  using 
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Control  Movements  for  Normal  Acceleration 
Step  Responses 


unity  pseudoinverse  weighting.  The  right  half  plane  zero  has 
been  shifted  to  the  stable  left  half  plane,  and  the  normal 
acceleration  response  has  been  "quickened"  by  the  maneuver 
flap. 


Summarizing  the  results  thus  far,  the  Type  1  discrete 
control  law  changes  closed-loop  eigenvalue  locations  for  changes 
in  the  quad  partition  matrix  pseudoinverse  weighting  matrix, 
while  the  Type  0  control  law  closed-loop  eigenvalues  do  not 
change.  Different  pseudoinverse  weighting  matrices  do  not 
significantly  change  the  effective  operating  range  of  the  Type 
1  control  law  but  can  desirably  change  the  operating  range  of 
the  Type  0  control  law.  Based  on  these  observations,  the  Type 
1  discrete  control  law  is  designed  with  unity  pseudoinverse 
weighting  as  is  done  in  Ref.  3  for  the  continuous  time  case, 
Unity  pseudoinverse  weighting  provides  good  closed-loop  eigen¬ 
value  locations  in  transforming  from  the  Type  0  to  Type  1  DFCS 
as  demonstrated  in  Tables  5  and  6.  The  discrete  time  Type  0 
DFCS  will  use  anc*  ^2  Senerated  using  unity  pseudo  inverse 
state  weighting  and  100/MAX  control  weighting.  This  weighting 
scheme  delay?  the  saturation  of  maneuver  flap  until  very  large 
A&n  commands  are  requested,  increasing  the  normal  operational 
range  of  the  Type  0  discrete  DFCS.  Note  however  that,  as  will 
be  shown  at  the  end  of  this  section,  favorably  saturating  ma¬ 
neuver  flap  does  not  decrease  the  operational  range  of  the 
Type  1  DFCS. 

The  discrete  DFCS  Type  0  and  Type  1  design  step  re¬ 
sponses  is  presented  in  Figs.  7  to  12.  Each  of  the  three  com¬ 
mands  are  individually  stepped  to  unity  and  the  perturbation 
state  and  control  responses  are  plotted.  Overlaid  on  the 
normal  acceleration  step  response  in  Figs.  7  and  10  are  the 
Category  II  boundaries  for  the  C*  response  taken  from  Ref.  17. 

C*  and  an  are  not  equally  comparable  but  the  boundary  docs 
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Figure  9  Type  1  DFCS  A0  Step  Response 
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Figure  10  Type  0  DFCS  Aan  Step  Response 

give  an  indication  of  the  approximate  performance  of  the  DFCS 
transients.  Overlaid  on  the  wind  axis  roll  rate  transient  in 
Figs.  8  and  11  is  the  Category  II  boundary  for  the  roll  rate 
response  also  taken  from  Ref.  17.  Again  p  and  pw  are  not  equal¬ 
ly  comparable  but  the  pw  transient  falls  well  within  the  enve¬ 
lope.  No  step  response  transient  envelope  is  available  for 
sideslip.  Sideslip  response  in  Figs.  8  and  11  during  a  roll 
rate  step  response  is  reduced  until  the  deviation  falls  within 
the  limits  specified  in  Ref.  15. 
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operating  range  of  the  aircraft.  In  the  Type  0  DFCS ,  a  trade¬ 
off  is  being  made  between  small  control  movements  and  speed  of 
response.  The  steady  state  values  for  A6  r  and  A6  for  a  0.305 
m/s z  (1  fps^)  Aan  command  using  the  Type  1  DFCS  are  0.710  deg 
and  -0.377  deg,  respectively,  as  shown  in  Fig.  7.  The  corre¬ 
sponding  Type  0  DFCS  steady  state  values  are  -0.321  for  A6 

s 

and  0.557  for  A6m£  in  Fig,  10.  It  would  appear  that  a  Aan 

command  above  A. 3  m/s^  (-1/2  g)  would  saturate  maneuver  flap 

rendering  the  Type  1  and  Type  0  control  laws  inoperative  for 

2 

Aan  commands  above  4.3  m/s  (1/2  g).  This  interpretation  is 
entirely  incorrect  for  the  discrete  Type  1  DFCS. 

2 

When  a  Aan  command  above  4,3  m/s  (1/2  g)  occurs,  the 

Type  1  control  law  rapidly  saturates  maneuver  flaps  at  10  deg 

of  travel.  The  Type  1  DFCS  reduces  to  the  situation  where 

there  is  one  control,  A 6.  which  is  to  be  commanded  by  the 

s 

control  law  to  meet  the  one  requirement  of  forcing  Aan  to  Aanc. 
The  capability  of  the  discrete  Type  1  DFCS  to  handle  control 
saturation  through  control  reset  is  a  unique  feature  of  the 
incremental  form  and  is  discussed  in  previous  paragraphs.  As 
long  as  one  longitudinal  control  is  not  saturated,  the  Type  1 
control  law  will  continually  try  to  zero  the  command  error 
with  the  unsaturated  control  surface.  The  Type  1  control  law 
fully  exploits  all  the  capabilities  of  the  aircraft  and  will 
accommodate  the  command  until  all  available  controls  saturate . 

A  similar  reset  mechanism  is  also  incorporated  in  the 
Type  0  control  law.  The  Type  0  control  law,  however,  does  not 
have  the  integral  property  and  is  unable  to  automatically  change 
control  strategy.  A  demonstration  of  the  Type  0  and  Type  1 
control  law  Aan  response  for  a  normal  acceleration  command  of 
12.2  m/s  (1.25  g)  is  shown  in  Fig.  13.  The  Type  0  control  law 
has  a  steady  state  error  while  the  Type  1  control  law  does  not . 
The  Type  1  maneuver  flap  response  unsaturates  because  stabilator 
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Effects  of  a  Normal  Acceleration  Command 
Which  Saturates  Maneuver  Flap 


position  compensates  for  the  saturation  and  at  the  same  time 
the  command  error  changes  sign  because  Aan  overshoots  the  com¬ 
mand.  The  resulting  sign  change  in  the  maneuver  flap  incre¬ 
mental  command,  Av^,  causes  the  control  position  for  A6m£  to 
reenter  operational  status.  The  Type  0  control  law  compares 
the  maneuver  flap  position  with  ^22^d'  ^e  TyPe  ®  command 
error  does  not  change  sign  and  maneuver  flaps  remains  saturat¬ 
ed.  Similar  behavior  in  the  nonlinear  simulations  are  shown 
in  Chapter  3. 


2.5  DESIGN  PERFORMANCE  AT  THIRTY  FLIGHT  CONDITIONS 

The  design  results  for  the  entire  range  of  flight 
conditions  are  presented  in  this  section.  The  remaining  twenty- 
nine  point  designs  are  performed  with  the  cost  function  weight¬ 
ing  elements  shown  in  Table  4  for  the  primary  design  held  essen¬ 
tially  constant.  The  closed-loop  mapped  eigenvalues  for  the 
Type  1  and  Type  0  DFCS  are  shown  in  Tables  7  and  8,  respectively. 
The  rise  time,  overshoot,  and  settling  time  for  each  command 
at  each  flight  condition  for  the  Type  1  DFCS  is  shown  in  Table 
9. 


The  rise  time  is  the  time  it  takes  the  response  to 
reach  90  percent  of  the  commanded  value.  The  overshoot  is  the 
maximum  peak  of  the  response  expressed  as  a  percentage  of  the 
command.  The  settling  time  is  the  time  required  for  the  re¬ 
sponse  to  settle  within  5  percent  of  the  commanded  value.  The 
step  response  characteristics  in  Table  9  show  little  overshoot 
and  acceptable  rise  times  with  behavior  much  like  those  in 
Figs.  7  to  12, 

There  is  some  eigenvalue  variation  in  transforming 
from  the  Type  0  (Table  8)  to  the  Type  1  (Table  7)  DFCS  but  the 
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TABLE  9 


100.0  I  0.55  I  -  I  0.70  I  2.85  I  -  I  3 .65  I  0-50  |  2.57. 


general  trend  is  to  improve  the  eigenvalue  location.  The  larg¬ 
est  eigenvalue  change  in  transforming  from  the  Type  0  to  the 
Type  1  DFCS  is  the  normal  acceleration  command  mode  which  be¬ 
comes  faster  (i.e.,  more  stable).  The  short  period  mode  exhi¬ 
bits  the  required  increase  in  natural  frequency  with  increased 
civnamic  pressure.  When  the  aircraft  is  rolling  (pw  /  0.0)  the 
short  period  damping  ratios  all  decrease  in  Table  7  and  8.  On 
the  other  hand  the  Dutch  roll  damping  improves  when  rolling. 
Except,  for  two  flight  conditions  (VQ  =  122  m/s,  oQ  =  17  deg, 
and  aQ  =  21.4  deg,  pWQ  -  100  deg/fv.c)  all  closed-loop  eigenval¬ 
ues  remain  with  the  requirements  specified  by  MIL-F-8785B . 

These  two  flight  conditions  are  considered  to  to  be  extreme 
and  the  transgression  in  minor  (short  period  damping  ratio  = 
0.33  instead  of  0.35). 


2.6  EFFECTS  OF  GAIN  SCHEDULING 

After  the  optimal  point  designs  are  complete,  the 
resulting  thirty  sets  of  gains  are  regressed  against  functions 
cf  the  trimmed  flight  conditions  in  order  to  find  highly  cor¬ 
related  relationships.  Four  function  sets  are  found  to  be 
highly  correlated  and  are  used  for  gain  scheduling.  They  are 


dO  Ofi(]  *) 

Gain  “  a, a  +  - +  aQa  +  a,  (Schedule  1) 

1  a  -30a  +  280  1311  * 


Gain  =  a^  V  +  ^2 +  a3an  +  a4 


(Schedule  2)  (16) 


Gain  =  «jPKV  +  a2pw0  +  a3pwan  + 


„  .  a2pw280  2 

Gam  =  alPua  -  ^  +  a3P„anV  +  a/. 


a  -30o  +  280 
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where  to  are  the  regression  coefficients  and  a^  is  the 

regression  constant.  These  gain  schedules  are  derived  from 

flight  conditions  at  one  altitude  (6,100  m,  20,000  ft)  but  the 

major  effects  of  altitude  can  be  incorporated  by  replacing  the 
2  12 

V  terms  with  j  p#  V  .  The  terra  is  the  local  air  density. 
Schedule  1  and  Schedule  2  are  for  the  primary  longitudinal  and 
lateral-directional  gains  while  the  latter  two  sets  are  for 
the  cross-coupling  gains  that  are  non-zero  when  pw  is  non-zero. 
Each  gain  is  scheduled  using  the  function  set  that  yields  the 
highest  correlation  coefficient. 

The  Type  0  DFCS  described  by  Eqs.  2  to  4  requires  80 
gains  to  be  scheduled.  The  Type  1  DFCS  described  by  Eq .  1 
requires  the  scheduling  of  only  40  gains.  The  scheduled  gains 
generally  have  correlation  coefficients,  p,  greater  than  0.8 
with  the  optimal  gains.  The  effect  of  the  gain  schedule  on 
closed-loop  eigenvalues  is  shown  in  Figs.  14  and  15  for  a  sweep 
in  angle  of  attack.  The  Type  1  DFCS  gain  schedule  design  re¬ 
mains  within  the  requirements  of  MIL-F-8785B  throughout  the 
angle  of  attack  sweep.  The  Type  0  DFCS,  however,  has  a  prob¬ 
lem  with  the  sideslip  command  mode  as  shown  in  Fig.  15.  Near 
16  deg  «o  the  closed-loop  mode  briefly  goes  unstable  using 
schedule  1. 

From  previous  work,  in  Refs.  1  and  3,  it  is  known  that 
the  fighter  aircraft  under  consideration  has  control  reversals 
occuring  in  the  lateral-directional  axis  near  16  deg  aQ.  An 
Important  part  of  this  study  is  to  determine  if  the  gain  sched¬ 
ule  constructed  here  can  accurately  capture  control  gain  sign 
changes.  In  summary  the  Type  1  gain  schedules  performed  very 
veil  while  the  initial  Type  0  gain  schedules  had  difficulties. 

The  instability  of  the  Ap  command  mode  in  Fig.  15  can 
be  traced  to  two  gains  which  are  plotted  in  Fig.  16b;  the  Av 
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Figure  14  Effects  of  the  Gain  Schedule  on  Closed- 
Loop  Type  1  DFCS  Mapped  Eigenvalues 

to  A6r  gain  and  the  Av  to  A6dg  gain.  Schedule  1  has  the  high¬ 
est  correlation  with  these  two  gains  and  was  originally  imple¬ 
mented.  Schedule  1  closely  matches  the  gains  except  near  16 
deg  ac.  Schedule  2  is  a  different  gain  schedule  shown  in  Eq . 
16  and  has  lower  correlation  with  the  optimal  gain.  Schedule 
2,  however,  maintains,  stability  as  shown  in  Fig.  15. 


Examples  of  Type  1  scheduled  gains  are  shown  in  Figs. 
16a  and  Fig.  17.  The  scheduled  gains  closely  match  the  op¬ 
timal  gains  particularly  near  16  deg  aQ.  The  Aw  to  A6g  gain 
in  Fig.  17  and  the  Av  to  Afir  gain  in  Fig.  16a  show  the  "bump" 
near  16  deg  aQ  with  a  roll  off  as  aQ  increases  that  the  func¬ 
tion  280/(a^-30a  +  280)  tries  to  capture.  The  Type  1  gain 
scheduled  design  step  responses  in  Fig.  18  at  the  primary  de¬ 
sign  point  are  essentially  the  same  as  the  optimal  control 
step  responses. 
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Figure  17  Control  Gain  Variation  With  Angle  of  Attack 
(V  =  183  m/s  (600  fps),  pw  =  0  deg/sec) 


2 . 7  CHAPTER  SUMMARY 

This  chapter  presents  the  design  of  Type  0  and  Type  1 
digital  control  laws  for  a  fighter  aircraft.  The  Type  1  design 
is  successful  and  is  suitable  for  full  evaluation  in  a  nonlinear 
simulation  as  is  done  in  Chapter  3.  Numerous  difficulties  are 
encountered  in  simultaneously  designing  the  Type  0  control 
law.  It  is  concluded  that  the  Type  0  control  structure  with 
the  low-pass  filter  effect  and  no  integral  compensation  is  not 
a  suitable  controller  for  the  fighter  aircraft. 
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a)  OPTIMAL  CONTROL  DESIGN 


b)  SCHEDULED  GAIN  DESIGN 


Figure  18  DFGS  Step  Response  Performance  at  the 
Design  Flight  Condition 
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Although  both  the  Type  0  and  Type  1  laws  are  direct 
digital  designs,  the  cost  function  weights  and  desired  closed- 
loop  eigenvalues  and  eigenvectors  are  chosen  using  continuous¬ 
time  specifications.  This  is  accomplished  by  using  the  sampled- 
data  regulator  and  by  mapping  the  resulting  discrete  closed-loop 
system  to  an  equivalent  continuous  system  using  the  log  of  a 
matrix. 


There  are  four  steps  in  the  Type  1  and  Type  0  DFCS 
designs  that  represent  the  "art"  in  the  design  procedure.  De¬ 
tailed  discussions  of  each  of  the  four  steps  are  presented  in 
this  chapter,  The  first  step  is  the  choice  of  the  pseudoin¬ 
verse  weighting  matrix  for  inverting  the  quad  partition  matrix 

~(*-I>  r” 

H  D 

used  in  calculating  the  steady-state  feedforward  matrix.  The 
quad  partition  matrix  is  not  square  because  three  commands  are 
driven  by  five  controls .  The  pseudoinverse  allows  for  the 
optimal  blending  of  steady-state  control  surfaces  to  satisfy 
the  command  requirements.  The  second  step  is  the  choice  of  the 
optimal  control  quadratic  weights  subject  to  what  constitutes 
a  "good"  design  for  a  fighter  aircraft.  The  third  step  is  the 
choice  of  the  pseudoinverse  weighting  matrix  for  inverting  the 
quad  partition  matrix  when  transforming  the  control  gains  from 
the  Type  0  to  Type  1  structure.  The  primary  objective  of  the 
third  step  is  to  reduce  any  adverse  closed- loop  eigenvalue 
variations  caused  by  the  transformation.  A  secondary  objective 
is  to  place  system  zeroes  in  more  desirable  locations.  The 
fourth  step  is  the  choice  of  the  schsduling  functions  used  in 
the  control  law  gain  schedules.  The  objective  is  to  develop 
gain  schedules  that  allow  control  gains  to  be  computed  as  a 
function  of  flight  condition  for  onboard  implementation  and 
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give  nearly  the  same  performance  as  the  optimal  control  point 
design  gains. 

Evaluation  of  the  gain  schedule  is  an  important  part 
of  the  design  because  the  reference  aircraft  is  known  to  have 
an  adverse  yaw  problem  at  high  angles  of  attack.  The  Type  1 
gain  schedule  effectively  handles  the  control  gain  sign  changes, 
while  the  initial  Type  0  gain  schedule  demonstrates  difficul¬ 
ties  which  are  subsequently  surmounted.  A  demonstration  of 
the  successful,  gain  scheduled  Type  1  DFCS  in  demanding  non¬ 
linear  maneuvers  is  presented  in  the  next  chapter. 


3. 


EVALUATION  OF  A  FIGHTER  AIRCRAFT  DIGITAL  FLIGHT 
CONTROL  SYSTEM  BY  NONLINEAR  SIMULATION 


3.1  OVERVIEW 

The  previous  chapter  describes  the  construction  of  a 
fighter  aircraft  digital  flight  control  system  (DFCS).  The 
system  is  designed  to  produce  accurate  and  well-damped  re¬ 
sponse  to  commands,  and  it  should  be  capable  of  operating  over 
a  wide  range  of  velocity,  normal  acceleration,  and  roll  rate. 
The  control  system  is  not  sensitive  to  control  surface  actua¬ 
tor  position  limiting,  as  occurs  in  air  combat  maneuvering 
(ACM)  flight.  To  test  all  of  these  attributes,  this  chapter 
details  the  construction,  execution,  and  results  of  a  series 
of  experiments  which  compare  the  ACM  tracking  accuracy  of  the 
digital  flight  control  system  to  the  conventional  control  sys¬ 
tem  presently  on-board  the  aircraft. 

The  comparative  tests,  structured  as  shown  in  Fig.  19, 
are  based  on  a  simulation  which  is  identical  from  experiment, 
to  experiment  except  for  the  control  system  structure.  The 
same  ACM  tasks  flown  by  the  same  pilot  model  are  used,  and  the 
same  nonlinear  6  degree-of- freedom  aircraft  model  is  used. 

The  full  authority  digital  flight  control  system  is  primarily 
compared  to  the  conventional  mechanical  linkages  with  limited- 
authority  analog  sas.  A  third  control  system  consisting  of 
the  mechanical  linkages  and  digitally  driven,  limited-authority 
sas  actuators  is  given  limited  testing.  Section  3.2  discusses 
the  details  of  the  pilot  and  aircraft  simulation,  while  Section 
3.3  presents  the  control  system  models.  The  test  results  are 
included  in  Sections  3.4  through  3.6. 
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Figure  19  Test  Structure 


3.2  SIMULATION  DESCRIPTION  AND  TEST  PROCEDURE 

The  components  of  the  ACM  simulation  are  the  nonlin¬ 
ear,  6-degree-of- freedom  (DOF)  aircraft  simulation,  the  ACM 
pilot  model,  the  ACM  task  definition,  and  the  aircraft  control 
system.  All  except  the  latter  are  discussed  in  this  section. 

It  is  important  to  emphasize  that  in  each  of  the  experiments, 
only  the  control  system  is  changed.  This  allows  a  direct  com¬ 
parison  of  the  control  system  effectiveness. 


3.2.1  Nonlinear  6-DOF  Aircraft  Simulation 


The  aircraft  dynamics  are  modeled  by  the  nonlinear,  6 
degree-of- freedom  rigid-body  equations  of  motion.  The  six 
kinematics  equations  are  defined  as 
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where  xE  is  the  aircraft  position  vector,  Ug  is  the  aircraft 
Euler  angle  vector,  Vg  is  the  aircraft  body-axis  velocity  vec 
tor,  Wg  is  the  aircraft  body-axis  angular  rate  vector,  Hg  is 
the  body- to-earth  axis  direction  cosine  matrix  and  Lg*  is  the 
body-to-Euler  angular  rate  transformation  matrix.  The  corre¬ 
sponding  six  dynamics  equations  complete  the  12  rigid  body 
state  equations 
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F/m  +  Hg  gE  -  wBvB 

1B1-  "  1Blu)BIB-B 

(18) 


where  F  and  M  are  the  force  and  moment  vectors  due  to  aero¬ 
dynamics  and  thrust  effects,  and  (~)  represents  the  cross- 
product  equivalent  operation.  The  aircraft  mass  is  m,  and  the 
rotational  inertia  matrix  is  Ig.  The  cross  product  terms  model 
the  inertial  couplings  which  are  important  in  maneuvering  flight 


The  aerodynamic  forces  and  moments  are  nonlinear  func¬ 
tions  of  the  state  of  motion  of  the  aircraft  and  the  control 
surface  deflections.  The  thrust  forces  and  moments  (which  are 
entirely  within  the  plane  of  symmetry)  are  functions  of  the 
flight  condition  and  throttle  position.  The  following  equa¬ 
tions  indicate  the  effects  included  in  this  model: 


^  s  cx(a  *P  >Q>6s,(5mf  ,(5sp^ 

T(M,h,6rp)  COS  6rp 

F  = 

da  S  Cy(a,p,p)?,6gp,6ds,6r) 

4* 

0 

S  cz(M,a,p,q,6s,6mf,dsp)^ 

-T(M,h,6T)  sin  0T 

(19) 
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0 
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The  quantities  0p  and  z-p  are  the  thrust  misalignment  and  mo¬ 
ment  arm,  respectively.  The  functional  relationships  implied 
in  Eqs .  19  and  20  are  given  primarily  by  tabular  data,  and  the 
aircraft  simulation  uses  an  efficient  linear  interpolation 
routine  to  calculate  the  aerodynamic  and  thrust  forces  and 
moments  at  any  specific  state. 

Control  surface  position  limits  are  included  in  the 
aircraft  model.  First-order,  rate-limited  actuator  models  are 
used  in  preliminary  testing,  but  are  not  included  in  most  pro¬ 
duction  aircraft  simulations  for  reasons  of  simulation  effi¬ 
ciency.  This  simplification  should  not  have  a  significant 
impact  on  conclusions  based  on  a  comparison  of  different  flight 
control  techniques  of  comparable  bandwidth.  Table  10  lists 
the  control  channels  and  the  displacement  limits  used  in  all 
the  experiments  reported  in  this  chapter,  along  with  the  rate 
limits  and  actuator  time  constants  that  can  be  used. 

3.2.2  ACM  Pilot  Model 

The  ACM  pilot  model  contains  an  acceleration-oriented 
trajectory  generator  that  describes  the  specific  maneuver  the 
pilot  wishes  to  perform,  and  a  set  of  piloting  gains  that  trans 
late  the  ACM  trajectory  into  stick  and  rudder  pedal  displace¬ 
ments.  Hence,  this  pilot  model  is  a  trajectory  tracking  pilot 
model .  Alternatively,  the  ACM  trajectory  internal  to  the  pilot 
model  can  be  viewed  as  the  maneuver  the  pilot  must  perform  to 
accurately  track  a  specific  opponent's  maneuver,  The  extension 
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TABLE  10 


CONTROL  CHANNEL  DETAILS 


T-3115 


CONTROL 

TOSITION  LIMITS 

RATE  LIMIT 

TIME  CONSTANT 

Throttle.  6^ 

0  to  2  units 

0.4  unit/ssc 

0.5  sec 

Stabilator,  fl 

-30  to  10  deg 

30  deg/sec 

0.05  sec 

Maneuver  Flap, 

0  to  10  deg 

30  deg/sec 

0.05  sec 

Spoiler,  6gp 

-55  to  55  deg 

250  deg/sec 

0.05  sec 

Differential  Stabilator,  6^ 

-12  to  12  deg 

30  deg/sec 

0.05  sec 

Rudder,  4 

-30  to  30  deg 

30  deg/sec 

0.05  nee 

• 

of  this  model  to  include  a  target  model,  tracking  kinematics 
and  right  dynamics  is  possible  but  unnecessarily  complex  for 
the  purposes  of  the  tests  reported  in  this  chapter 


The  ACM  trajectory  generator  commands  consist  of  the 
three  components  of  earth  relative  acceleration  expressed  in  wind 
axes  (V  ,  a  ,  a  )  and  the  component  of  body  angular  rate  along 
the  velocity  vector  (p  ).  Wind  axes  consist  of  an  x-axis  along 
the  velocity  vector  and  a  z-axis  in  the  body  x-z  plane.  The  roll 
angle  about  the  velocity  vector  between  the  wind  z-axis  and  the 
vertical  plane  is  denoted  by  <|>  ,  which  forms  one  of  the  ACM  tra¬ 
jectory  generator  states.  The  others  are  velocity  magnitude 
(V),  flight  path  angle  (y),  and  velocity  heading  (£,).  The 
nonlinear  state  equations  that  describe  the  trajectory  gen¬ 
erator  dynamics  are  given  in  Eq  21. 


r.  " 

r  -  I 

V 

v 

c 

K 

p  +  t  sin  y  j 

= 

c 

i 

(1/V  cos  y )  ( anc  sin<frv  +  a  co3*v) 

• 

Y 

(1/V)(anc  cos<^v  -  ayc  sin<t>v) 
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Note  again  that  anc  and  a  are  components  of  earth-relative 
acceleration,  and  differ  from  ideal  accelerometer  outputs  in 
these  axes  by  the  acceleration  of  gravity. 

The  transformation  of  the  trajectory  generator  out¬ 
puts  to  stick  and  pedal  deflections  occurs  by  two  routes,  as 
illustrated  in  Fig.  20.  The  ACM  commands  themselves  (acceler¬ 
ations  and  roll  rate)  provide  pilot  control  deflections  direct¬ 
ly  by  way  of  the  command/control  pilot  gain  matrix.  These 
inputs  are  "open- loop"  in  the  sense  that  they  do  not  depend  on 
the  actual  aircraft  path- following  accuracy.  These  open-loop 
inputs  will  dominate  the  pilot's  input  over  the  short  term  and 
in  very  rapid  maneuvers,  such  as  a  rolling  reversal.  Command/ 
control  gains  for  the  ACM  pilot  model  are  derived  primarily 
from  the  command/control  steady-state  ft  matrices  discussed  in 
the  previous  chapter  of  this  report.  Some  modif ications  to 
these  gains  were  made  based  on  initial  simulation  tests. 


ACTUM  VELOCITY 
ROLL  iNCLE 

Figure  20  ACM  Pilot  Model 


By  differencing  the  ideal  ACM  vel.ocity  and  roll  angle 
(V,  <(>  ,  £ ,  y)c  with  the  actual  aircraft  values,  the  maneuver 
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errors  can  provide  error-correcting  pilot  inputs  (Fig.  20).  In 
principle,  velocity  error  drives  throttle,  roll  error  drives 
lateral  stick,  flight  path  angle  error  (transformed  into  air¬ 
craft  axes)  drives  longitudinal  stick  and  heading  error  (trans¬ 
formed  into  aircraft  axes)  can  drive  rudder  pedal.  In  actual 
implementation  in  this  study,  since  lateral  acceleration  (di¬ 
rect  sideforce)  is  not  within  the  capability  of  this  airframe 
to  any  significant  extent,  ayc  is  always  zero,  and  sideslip 
error  is  used  to  drive  rudder. 

The  error  feedback  gains  are  derived  from  the  command 
gains  by  assuming  a  2-sec  time  constant  for  error  correction. 
Table  11  lists  the  pilot  model  gains  for  the  high  subsonic 
flight  condition  at  which  testing  is  performed.  Note  that  the 
large  throttle  error  gain  is  chosen  so  that  the  pilot  model 
selects  full  afterburner  early  in  an  ACM  maneuver,  which 
minimizes  energy  loss.  More  detail  concerning  the  rudder 
command  gain  is  given  in  the  results  sections  later  in  this 
chapter;  the  rudder  input  is  necessary  to  roll  the  aircraft  at 

TABLE  11 

ACM  PILOT  MODEL  GAINS 


CONTROL 

COMMAND 

KRROR 

t— — - 

POSITION  LIMIT 

Tht  nt  i  Jr 

A.  5  units  A.j.'Vfi 

1 - 

0.07 

en  1 1  h  A^./fps  av 

0  t  o 

2  unit  s'' 

l.onjt i  t  ud  1  tin  1 

St  Irk 

0.7S8  In 

0.  161 

»  1"  *,onAI<iR  A, 

-4 . 0 

to  5.5  In, 

-10,2  cm  li»  10.  cm 

Lai  ci'fll  Stick 

0.04  in  A,at/deR/sec  pw 

c 

0.02 

in  fi ] n 

-4.5 

to  4.5  i  ii . 

•11.4  cm  m  |),4  cm 

Rudder 

0.0 

or 

-0.1 

ir  f'ped/<lrR  4,1 

-:i.o 

to  3.0  in, 

-7.6  cm  to  7.6  cm 

0.012  in  6,,,.,,/deR/sFC  pu 

*  Tli  ro  1 1 1  c  units  0,0 
1.0 
2,0 


idle 

full  mi  1 1 1  ary  power 
full  afterburner 


high  angle  of  attack  with  a  conventional  control  system.  The 
table  also  lists  the  pilot  control  deflection  limits. 

3.2.3  ACM  Tracking  Tasks 

There  are  two  tracking  tasks  used  in  this  study,  a 
wind-up  turn  and  a  rolling  reversal.  They  are  specified  in 
terms  of  the  velocity-axis  accelerations  and  roll  rates  in¬ 
volved.  Figure  21  illustrates  the  time  history  of  the  com¬ 
manded  wind-up  turn.  The  commanded  longitudinal  acceleration 
is  zero.  This  set  of  commands  results  in  a  rapid  roll  into  a 
steady,  high-g  right  turn  followed,  10  sec  later,  by  a  rapid 
roll  into  a  high-g  left  turn.  The  actual  capabilities  of  the 
aircraft  are  such  that  this  maneuver  exceeds  the  steady  turn 
rate  at  the  flight  condition  chosen  (800  fps  TAS  at  20,000  ft, 
244  m/s  at  6100  m).  Hence,  the  aircraft  will  slow  down  in 
this  maneuver  and  increasing  or  will  be  required  to  maintain 
the  high-g  flight  condition.  The  net  effect  of  all  this  is 
that  the  bank-to-bank  maneuver  10  sec  into  the  trajectory  will 


n.4?H 


Figure  21 


Wind-up  Turn  Trajectory  Inputs 


occur  at  high  angle  of  attack.  This  maneuver  is  flown  both  by 
an  open-loop  pilot  model  (all  error  feedback  gains  set  to  zero) 
and  a  closed-loop  pilot  model,  which  includes  all  of  the  error 
feedback  gains.  The  latter  case  is  a  form  of  target  tracking. 
In  all  of  these  tests  the  roll  command- to-rudder  gain  is  zero. 

The  second  maneuver  used  in  these  tests  is  described 
by  the  inputs  illustrated  in  Fig.  22.  This  maneuver  is  called 
a  rolling  reversal,  and  is  designed  to  convert  ("reverse")  a 
defensive  tactical  situation  into  an  offensive  one.  The  man¬ 
euver  is  a  combination  of  high-g  pull-up  and  360  deg  roll. 

The  pull-up  starts  1  sec  before  the  roll.  This  combination 
results  in  a  trajectory  which  rises  sharply  and  then  dives  (by 
means  of  a  sharp  pull-up  while  rolling  inverted)  back  to  rough¬ 
ly  the  original  trajectory  Since  this  is  a  "pre-programmed" 
maneuver,  these  tests  are  performed  without  any  error  feedback 
in  the  pilot  model.  This  maneuver  is  so  rapid  that  velocity 
changes  are  not  as  severe  as  in  the  wind-up  turn.  One  diffi¬ 
culty  in  flying  this  trajectory  involves  the  high  angle-of- 
attack  roll  effectiveness  of  lateral  stick  in  the  conventional 
aircraft.  It  is  low  enough  so  that  lateral  stick  alone  will 
not  lead  to  the  completion  of  the  maneuver.  For  this  reason, 
the  roll  rate  command- to-pedals  gain  is  inserted  in  the  ACM 
pilot,  model.  The  DFCS  does  not  need  the  pedal  input  to  follow 
this  maneuver,  as  will  be  shown  in  Section  3.5. 

3.3  CONTROI,  SYSTEM  CANDIDATES 

Three  control  systems  are  tested  in  the  experiments 
reported  in  the  following  sections.  Figure  19  shows  that  they 
all  interconnect  the  pilot’s  control  stick  and  pedals  with  the 
aircraft  control  surfaces.  Thi*ottle  and  engine  control  are 
not  discussed  further  in  this  section,  since  the  trajectories 
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Figure  22  Rolling  Reversal  Trajectory  Inputs 

that  are  tested  require  full  afterburner  to  minimize  energy 
loss  during  the  high-g  maneuvers. 

The  conventional  mechanical  control  system  is  illus¬ 
trated  in  Fig.  23,  and  it  consists  of  three  separate  channels. 
The  pilot's  lateral  stick  also  drives  the  spoilers,  and  the 
maneuver  flap  is  not  deflected.  The  pitch  and  yaw  sas  actua¬ 
tors  are  limited-authority  channels  driven  by  washed-out  rate 
feedbacks.  Yaw  sas  also  involves  lateral  acceleration  feed¬ 
back  for  turn  coordination,  The  pitch  sas  limit  is  ±3  deg 
stabilator,  while  the  yaw  sas  limit  used  in  this  study  is  £9.5 
deg  rudder.  The  variable  limit  on  pilot's  rudder  input  is 
about  ±15  deg  in  the  ACM  flight  condition. 

The  roll  channel  contains  a  rate  feedback  and  a  com¬ 
mand  augmentation  feedforward.  The  pilot's  differential  sta¬ 
bilator  can  never  be  larger  than  ±7  deg,  but  the  roll  cas  can 
add  up  to  5  deg  of  differential  stabilator  to  this.  Hence, 
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neither  the  pilot  nor  the  cas  have  full  differential  stabila- 
tor  authority,  and  in  a  case  of  roll  cas  failure  or  shutdown, 
the  effective  pilot  roll  control  power  (due  to  both  the  me¬ 
chanical  linkage  and  the  roil  cas)  is  reduced  considerably. 

The  roll  cas  signal  is  composed  of  a  roll  rate  feedback  and  a 
lagged  feedforward  of  pilot  lateral  stick.  This  roll  command 
augmentation  produces  high  effective  pilot  roll  control  power 
while  avoiding  the  large  roll  accelerations  that  would  result 
if  the  mechanical  channel  had  full  differential  stabilator 
authority.  Any  or  all  of  the  three  independent  sas  channels 
can  be  turned  off. 

The  inputs  to  the  analog  sas  computer  consist  of  nor¬ 
mal  and  lateral  acceleration  and  the  roll,  pitch  and  yaw  angular 
rates.  Later  modifications  to  the  conventional  control  system 
(described  in  Ref.  21  but  not  modeled  here)  also  use  angle  of 
attack  and  Mach  number  to  schedule  automatic  maneuver  flap  and 
aileron-rudder  interconnect.  The  nose  probe  that  senses  angle 
of  attack  can  also  sense  sideslip,  although  this  capability  is 
not  presently  utilized.  Another  input  to  the  sas  computer  is 
pilot  lateral  stick  position,  which  drives  the  roll  cas  lag 
compensation.  The  outputs  from  the  analog  sas  computer  are 
commands  to  the  rudder,  stabilator  and  differential  stabilator 
sas  actuators. 

Another  of  the  flight  control  system  candidates  is 
the  full-authority  digital  flight  control  system  (DFCS),  de¬ 
scribed  in  Chapter  2.  As  illustrated  in  Fig.  24,  all  control 
calculations  occur  in  a  digital  flight  computer  (actually  a 
set  of  parallel  computers  for  reliability  reasons)  and  the 
outputs  command  full-authority  actuators  for  stabilator,  man¬ 
euver  flap,  spoiler,  differential  stabilator,  and  rudder.  The 
command  inputs  are  pilot  longitudinal  and  lateral  stick  posi¬ 
tion  and  rudder  pedal  position,  which  are  interpreted  in  the 
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flight  computer  as  commands  of  normal  acceleration,  roll  rate 
about  the  velocity  vector  and  sideslip.  Although  these  pilot- 
control  -to-command  relations  would  be  scheduled  wirh  flight 
condition  in  an  actual  implementation,  they  are  fixed  for  the 
purposes  of  this  study  at  1.32g 's/in,  25  deg/sec/in  and  -10 
deg/in  respectively.  The  inputs  to  the  digital  flight  com¬ 
puter  are  similar  to  those  available  to  the  conventional  ana¬ 
log  sas  computer  in  that  normal  acceleration,  angle  of  attack, 
and  roll,  pitch  and  yaw  rates  are  needed.  Additionally,  side¬ 
slip  is  required  for  the  control  algorithm  and  dynamic  pres¬ 
sure  is  needed  for  gain  scheduling.  In  principle,  sideslip 
command  and  feedback  can  be  eliminated  in  favor  of  lateral 
accelerometer  feedback  and  lateral  specific  contact  force  com¬ 
mand,  resulting  in  inputs  essentially  identical  to  those  used 
by  the  conventional  control  system. 

The  third  control  system  tested  in  the  ACM  experi¬ 
ments  is  a  digitally-augmented  conventional  control  system, 
and  it  is  illustrated  in  Fig.  25.  In  this  control  system,  the 
conventional  mechanical  linkages  between  pilot  and  control  sur¬ 
face  main  actuators  are  unchanged,  and  the  limited-authority 
sas  actuators  are  retained.  The  sas  computer,  however,  is 
digital  and  contains  as  a  subprogram  the  complete  digital 
flight  control  system  described  above.  The  sas  actuator  com¬ 
mands  are  calculated  by  subtracting  the  pilot's  mechanical 
linkage  control  surface  command  from  the  DFCS  control  surface 
command,  much  as  is  done  in  the  present  ARI  design  for  the 
F-14A  (Ref,  22).  In  the  absence  of  control  signal  limiting 
effects,  the  main  actuator  command  signal  is  identical  to  the 
DFCS  actuator  command  signal.  The  value  of  such  a  control  sys¬ 
tem  is  obvious  in  the  context  of  a  retrofit  program,  and  the 
reliability  and  redundancy  issues  are  more  akin  to  those  in¬ 
volved  with  conventional  control  systems.  The  basic  question 
to  be  answered  about  a  digitally-augmented  system  is  whether 
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Figure  25  Digital Iy~Augmen ted  Conventional  Control  System 


the  limited-authority  sas  channels  are  sufficiently  powerful 
to  improve  the  aircraft  command  response  to  a  level  approach¬ 
ing  the  full -authority  DFCS. 

Testing  of  the  conventional  control  system  and  the 
DFCS  against  each  other  is  described  in  the  following  two  sec¬ 
tions,  and  the  performance  of  the  digitally-augmented  conven¬ 
tional  control  system  is  examined  in  Section  3.6. 


3.4  TRAJECTORY  TRACKING  IN  THE  WIND-UP  TURN 

The  wind-up  turn  maneuver  implemented  here  involves  a 
rapid  roll  into  a  steady  right  turn,  followed  10  sec  later  by 
a  roll  to  a  steady  left  turn.  Although  maximum  thrust  is  ap¬ 
plied  ,  the  aircraft  slows,  which  requires  higher  a  to  hold 
the  large  normal  acceleration.  The  emphasis  in  this  test  is 
on  tracking  accuracy,  so  a  closed-loop  pilot  model  is  used. 

As  discussed  in  Section  3.3.1,  this  pilot  model  attempts  to 
null  trajectory  following  errors  by  moving  the  stick  and  pedals 
in  response  to  these  errors.  The  pilot's  roll  rate-to-rudder 
gain  is  zero  for  these  wind-up  turn  experiments. 

Two  control  systems  are  tested:  the  conventional 
mechanical  control  system  (sas  off)  and  the  full -authority 
DFCS.  The  command  responses  of  the  DFCS  are  illustrated  by 
Fig.  26,  which  shows  that  the  aircraft  response  to  the  step 
commands  of  an  and  pw  is  rapid  and  does  not  exhibit  excessive 
overshoot.  The  response  settling  is  very  good.  Although  the 
sideslip  command  is  zero,  there  is  a  significant  response  dur¬ 
ing  the  initial  roll  and  pull-up.  The  sideslip  settles  quickly 
to  the  steady  value  after  the  initial  transient,  and  exhibits 
very  little  oscillation.  Note  that  the  structure  of  the  DFCS 
is  such  (Fig.  20)  that  the  command  trajectories  are  proportional 
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to  the  pilot's  control  inputs,  The  normal  acceleration  and 
wind-axis  roll  rate  command  inputs  come  from  longitudinal  and 
lateral  stick  respectively,  and  sideslip  command  comes  from 
rudder  pedals.  The  angle-of-attack  response  is  also  included 
in  Fig.  26  to  show  its  rapid  yet  well  controlled  character. 

Note  that  as  speed  drops  during  the  maneuver,  angle  of  attack 
is  increased  by  the  DFCS  to  hold  the  commanded  normal  accelera¬ 
tion. 


TIME  IimI  TIME  (ucl 

Figure  26  DFCS  Command  Response  in  the  Wind-up  Turn 
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The  aircraft  control  surface  motions  that  cause  this 
response  are  shown  in  Fig.  27.  The  rapid  stabilator  and  man¬ 
euver  flap  movement  initially  cause  the  rapid  response  noted 
above,  yet  both  control  surfaces  are  quickly  reduced  to  lower 
values  representative  of  the  steady-state  turn.  Note  that 
full  maneuver  flap  is  used  to  produce  maneuver  quickening. 

Later  in  this  maneuver,  the  stabilator  deflection  is  increased 
to  achieve  higher  angle  of  attack.  The  controller  also  de¬ 
flects  maneuver  flap  to  assist  holding  normal  acceleration  as 
speed  is  reduced.  Hence,  the  DFCS  includes  an  automatic  man¬ 
euver  flap  feature.  Differential  stabilator  is  commanded  to 
its  limit  as  the  aircraft  rolls  into  the  turn,  and  then  quick¬ 
ly  is  reduced  to  a  low  value.  Some  motion  is  necessary  to 
damp  the  Dutch  roll  oscillation,  which  in  this  aircraft  in¬ 
volves  significant  roll  motion  (Ref.  23).  Rudder  motion  serves 
to  damp  sideslip  motion,  and  spoiler  deflection  is  relatively 
small  and  hence  not  shown  here. 

The  same  ACM  pilot  model  flying  the  same  airframe  in 
the  same  maneuver  with  a  conventional  mechanical  control  system 
(sas  off)  achieves  much  different  results.  Figure  28  compares 
roll  angle  responses.  The  poor  Dutch  roll  damping  of  the  con¬ 
ventionally-controlled  aircraft  is  apparent.  The  sideslip 
response  (Fig.  29)  confirms  this.  The  Initial  sideslip  pulse 
is  similar  to  that  observed  for  the  DFCS  aircraft  (Fig.  26), 
but  the  Dutch  roll  that  the  sideslip  pulse  excites  lasts  much 
longer  than  for  the  DFCS  aircraft.  The  angle-of-attack  response 
shown  in  Fig,  29  should  be  compared  to  that  shown  in  Fig.  26; 
this  comparison  shows  the  low-overshoot  and  well-dampled  nature 
of  the  DFCS  response.  Note  that  the  constant  a  characteristics 
of  the  conventional  control  system  is  not  necessarily  a  negative 
attribute;  the  conventional  control  system  is  approximately  an 
angle-of-attack  command  system,  and  the  pilot,  is  expected  to 
compensate  for  the  loss  of  normal  acceleration  at  a  given  a  as 
the  speed  decreases. 
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Figure  27 
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DFCS  Control  Surface  Trajectories  in 
the  Wind-up  Turn 


The  stabilator  trajectories  that  result  from  the  con 
V'-ntional  control  system  are  shown  in  Fig,  30,  and  these  are 
proportional  to  the  pilot  model's  control  stick  inputs.  Es¬ 
pecially  apparent  is  the  pilot  model's  attempt  to  reduce  the 
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DFCS  Control  Surface  Trajectories  in 
the  Wind-up  Turn 


The  stabilator  trajectories  that  result  from  the  con 
ventional  control  system  are  shown  in  Fig.  30,  and  these  are 
proportional  to  the  pilot  model's  control  stick  inputs.  Es¬ 
pecially  apparent  is  the  pilot  model's  attempt  to  reduce  the 
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J.  poorly-damped  Dutch  roll  response  by  using  lateral  stick  (dif¬ 

ferential  stabilator) .  If  this  short-term  compensation  task 
]  can  be  performed  automatically,  the  pilot  could  devote  more 

attention  to  target  tracking. 

[ 

All  of  the  previous  plots  deal  with  the  turn  entry 
I'  and  steady  right  turn  during  the  first  10  sec  of  the  trajec- 

*■  tory.  At  10.1  sec,  the  pilot  model  attempts  to  roll  the  air- 

*  craft  to  a  left  bank,  and  this  roll  input  occurs  at  high  angle 

1  of  attack  for  both  control  systems.  Speed  has  dropped  from 

the  maneuver  entry  speed  of  0.77  Mach  to  0.62  for  the  conven- 
j  tionally-controlled  aircraft  and  speed  has  dropped  to  0.47 

Mach  for  the  DFCS  aircraft,  due  to  its  higher  and  more  ac- 

1  ■ 

j  curate  normal  acceleration  response,  which  has  unavoidably 

produced  higher  drag.  Indeed,  the  DFCS  aircraft  reaches  the 
|  aircraft  pitch  performance  boundary  soon  after  the  left  roll 

is  initiated  because  stabilator  goes  to  its  aircraft-nose-up 
\  limit.  The  aircraft  remains  well  controlled  about  the  other 

‘  axes,  however,  and  even  angle  of  attack  exhibits  only  a  mild 

1  oscillatory  response  around  31  deg  a.  The  roll  to  the  left 

1  bank  (Fig.  31)  is  sluggish  due  to  the  low  velocity  but  occurs 

accurately.  The  mild  oscillations  in  roll  are  probably  due  to 
|  the  small  uncontrolled  a  oscillations. 

j  The  same  maneuver  (entered  at  somewhat  lower  angle  of 

attack)  in  the  conventionally-controlled  aircraft  results  in  a 
j  control-induced  departure  and  incipient  spin.  Angle  of  attack 

is  oscillatory  and  approaches  50  deg.  Among  the  best  indica- 
|  tions  of  the  departure  are  the  plots  shown  in  Fig.  32,  where 

*  body  yaw  rate  is  rapidly  increasing.  The  accompanying  roll 

1  angle  trajectory  shows  that  the  aircraft  initially  rolls  left, 

I  as  desired,  but  rapidly  snap-rolls  to  the  right  through  more 

than  360  deg.  This  undesired  motion  is  caused  by  the  pilot 
!  lateral  stick  input  which  drives  differentia]  stabilator. 
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Figure  30 


Stabilator  Time  Histories  For  a  Conventionally 
Controlled  Aircraft  in  a  Wind-up  Turn 
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Figure  31  DFCS  Aircraft  in  Bank-to-Bank  At 
High  Angle  of  Attack 

At  high  a ,  this  produces  roll  and  yaw  moments  of  'he  opposite 
sign  (in  this  case  left  roll  moment  and  right  yaw  moment), 
whereas  roll  and  yaw  moments  of  the  same  sign  are  necessary  to 
successfully  roll  about  the  velocity  vector.  The  adverse  yaw 
response  causes  the  undesired  snap  roll.  There  is  no  indica¬ 
tion  at  20  sec  that  the  pilot  has  recovered  control  of  the 
aircraft,  and  altitude  is  being  lost  at  16,000  feet  per  minute 
(81  m/s ) . 

The  DFCS  has  enabled  the  pilot  model  to  fly  the  air¬ 
craft  right  up  to  its  performance  boundary  without  experienc¬ 
ing  a  departure  from  controlled  flight.  Besides  expanding  the 
useful  performance  region,  this  has  a  significant  impact  on 
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Figure  32  Control  Induced  Departure  for 
Conventional  Control  System 
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pilot  confidence  and  pilot  training.  With  the  conventional 
control  system,  the  departure  shown  here  can  be  avoided  if  the 
pilot  learns  to  center  lateral  stick  at  high  angles  of  attack. 
Even  if  the  pilot  could  remember  to  do  this  along  with  his 
other  ACM  tasks,  it  is  not  clear  he  has  sufficient  information 
on  which  to  change  his  control  strategy.  The  DFCS  assumes 
this  task,  allowing  the  pilot  to  prosecute  the  ACM  task  with 
improved  concentration. 


3.5  ACCURACY  IN  THE  ROLLING  REVERSAL 

The  rolling  reversal  is  an  air  combat  maneuver  de¬ 
signed  to  convert  ("reverse")  a  defensive  tactical  situation 
to  an  offensive  one.  The  maneuver  is  a  combination  of  a  high-g 
pull-up  accompanied  by  a  360  deg  roll.  Since  it  is  a  "pre¬ 
programmed"  maneuver,  these  tests  are  performed  without  any 
error  feedback  in  the  pilot  model.  The  maneuver  is  so  rapid 
that  velocity  changes  are  not  as  severe  as  in  the  wind-up  turn. 

As  discussed  in  Section  3.3.3,  Initial  tests  indicat¬ 
ed  that  the  conventionally-controlled  aircraft  cannot  execute 
this  maneuver  with  only  lateral  stick  input  due  to  the  low 
rolling  moment  at  high  a.  For  this  reason,  the  roll  rate  com¬ 
mand-  to-pedals  gain  is  inserted  in  the  ACM  pilot  model.  The 
DFCS  does  not  need  pedal  input  to  follow  this  maneuver,  and 
testing  indicates  that  the  only  effect  of  the  extra  pedal  input 
is  a  significant  increase  in  sideslip  (since  pedal  is  the  DFCS 
sideslip  command  channel).  Hence,  in  the  comparison  presented 
here,  the  DFCS  pilot  model  is  simpler  than  the  pilot  model 
flying  the  conventional  aircraft.  In  the  following  compari¬ 
sons,  the  conventional  aircraft  has  all  three  sas  channels  on. 
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The  roll  rate  accuracy  through  the  360  deg  roll  is 
shown  in  Fig.  33.  The  conventional  sas  results  in  only  a  ±  20 
deg/sec  accuracy  to  the  60  deg/sec  command,  with  a  significant 
attitude  effect  apparent,  i.e.,  the  roll  rate  drops  as  the 
aircraft  is  inverted.  The  DFCS  roll  rate  lies  within  ±4  deg/ 
sec  of  the  command,  even  though  rapid  changes  in  the  flight 
condition  (roll  rate  and  normal  acceleration)  are  taking  place. 
The  end  of  the  maneuver  is  much  more  rapid  and  accurate  for 
the  DFCS  aircraft,  whereas  the  conventionally-controlled  air¬ 
craft  roll  response  tails  off  slowly.  This  will  produce  a 
roll  angle  error  which  the  open-loop  pilot  model  is  not  set  up 
to  correct. 


The  longitudinal  responses  are  compared  in  Fig.  34, 
which  shows  that  the  conventional  sas  response  and  the  DFCS 
response  are  quite  similar.  The  DFCS  produces  a  somewhat  more 
even  angle-of-attack  trajectory  and,  due  to  this,  a  much  better 
normal  acclereration  response.  The  increase  in  angle-of-attack 
at  the  end  of  the  roll  compensates  for  the  loss  in  velocity 
that  has  occured  during  the  maneuver. 

Yaw  rate  magnitudes  produced  by  the  two  control  sys¬ 
tems  are  comparable,  although  the  DFCS  produces  much  better 
sideslip  control,  as  illustrated  by  Fig.  35.  Sideslip  trans¬ 
ients  are  reduced  to  about  a  quarter  of  the  value  caused  by 
the  conventionally  controlled  aircraft.  There  is  nc  signif¬ 
icant  steady-state  value  for  the  DFCS  response,  whereas  the 
conventional  sas  allows  about  a  2  deg  average  sideslip  during 
the  60  deg/sec  roll. 

The  rudder  deflections  which  produce  these  responses 
are  shown  in  Fig.  36,  where  the  pilot  rudder  input  (due  to  the 
pilot  model's  internal  roll  rate- to-pedal  gain)  is  apparent  in 
the  conventional  sas  trajectory.  The  stability  augmentation 
system  countermands  the  pilot  input  early  in  the  roll,  which 
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Figure  35  Sideslip  in  The  Rolling  Reversal 
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causes  the  total  rudder  deflection  to  be  fairly  small.  During 
the  latter  half  of  the  roll,  the  augmentation  signal  is  reduced 
and  assists  the  pilot  in  completing  the  maneuver.  The  DFCS 
rudder  trace  is  also  shown,  and  its  similarity  to  the  actual 
rudder  deflection  due  to  the  pilot  and  conventional  sas  is 
notable,  This  ruder  deflection  produces  much  better  sideslip 
control,  as  noted  above,  and  is  generated  by  the  DFCS  from 
lateral  stick  input  only.  No  pilot  rudder  pedal  input  is  nec¬ 
essary  to  produce  proper  roll  response  and  accurate  sideslip 
control  in  the  DFCS  aircraft. 


3.6  DIGITALLY-AUGMENTED  CONVENTIONAL  CONTROL  SYSTEM 

During  the  test  program,  it  became  apparent  that  a 
promising  control  system  modification  approach  could  involve 
the  retention  of  the  conventional  mechanical  linkages,  but 
with  digital  control  of  the  existing,  limited-authority  sas 
actuators.  Ideally,  this  control  system  would  posess  all  of 
the  response  advantages  of  the  full -authority  DFCS  but  would 
be  irnplementable  in  existing  aircraft  by  replacing  the  analog 
sas  computer  with  a  digital,  one.  As  a  p  liminary  test  of  the 
response  of  such  a  system,  the  control  system  implementation 
shown  in  Fig.  22  is  chosen.  The  DFCS  control  laws  are  identi¬ 
cal  to  those  tested  in  the  previous  sections,  and  the  sas  ac¬ 
tuator  commands  are  formed  by  subtracting  the  pilot  control 
command  from  the  DFCS  control  command.  Hence,  for  stabilator, 
differential  stabilator,  and  rudder,  the  actual  control  sur¬ 
face  deflection  approximates  the  DFCS  control  command.  Of 
course,  if  the  difference  between  the  pilot  control  signal  and 
the  DFCS  control  signal  is  greater  than  the  sas  actuator  author¬ 
ity,  this  property  in  lost.  Additionally,  the  DFCS  maneuver 
flap  and  spoiler  signals  are  not  used;  maneuver  flap  is  held 
at  aero  and  spoiler  is  tied  to  lateral  stick  with  no  las  input. 


Finally,  it  should  be  noted  that  the  DFCS  control  law  does  not 
include  control  position  resets  calibrated  to  the  limited  au¬ 
thority  sas  actuators. 

The  digitally-augmented  conventional  control  system 
is  tested  in  the  rolling  reversal,  and  hence  the  results  are 
to  be  compared  to  those  presented  in  Section  3.5.  Roll  rate 
response  accuracy,  shown  in  Fig.  37  significantly  exceeds  that 
of  the  aircraft  with  conventional  sas  during  most  of  the  maneu¬ 
ver,  although  it  is  not  as  good  as  the  full-authority  DFCS. 

During  the  last  1.5  sec  of  the  roll  (5,5  sec  to  7.0 
sec)  the  roll  rate  drops  dramatically  from  its  commanded  value 
of  60  deg/sec.  Noting  that  this  part  of  the  roll  occurs  at 
high  a  and  that  rudder  is  especially  effective  in  rolling  the 
airplane  at  high  a,  it  is  concluded  that  the  9.5  deg  rudder 
sas  limit  used  in  this  simulation  is  not  large  enough  to  com¬ 
plete  this  maneuver.  This  is  confirmed  by  noting  that  the 
control  systems  tested  in  Section  3.5  used  13-18  deg  of  rudder 
at  this  point  in  the  maneuver  (see  Fig.  36).  The  roll  rate 
"bounce"  from  7  to  8.5  sec  is  due  to  the  lack  of  proper  DFCS 
control  reset.  It  is  expected  that  tailoring  the  control  re¬ 
sets  to  the  sas  actuator  limits  can  eliminate  this  "bounce" 
without  affecting  the  rest  of  the  response.  Angle  of  attack 
response  is  similar,  with  accuracy  better  than  a  conventional 
sas  but  not  as  good  as  the  full-authority  DFCS.  An  undesired 
transient  at  the  end  of  the  trajectory  is  apparent  here  also. 

It  is  expected  that  the  use  of  maneuver  flap  for  longitudinal 
response  quickening  can  improve  the  speed  and  accuracy  of  a 
digitally-augmented  conventionally-controlled  aircraft. 

Sideslip  response  is  shown  in  Fig,  38,  where  it  is 
aDparent  that  the  initial  sideslip  transient  is  almost  as  well 
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controlled  as  in  the  DFCS  aircraft.  The  large  spike  as  the 
aircraft  stops  rolling  at  7  sec  can  possibly  be  reduced  by 
proper  calibration  of  the  digital  augmentation  control  resets. 
The  rudder  response,  also  shown  in  Fig.  38,  illustrates  that 
the  9.5  deg  rudder  sas  limit  greatly  restricts  the  yaw  channel 
of  the  digital  augmentation  system.  The  general  rudder  re¬ 
sponse,  however,  is  quite  similar  to  that  achieved  by  the  full 
authority  DFCS,  as  can  be  seen  by  comparing  Fig.  36  and  38. 
Interestingly,  the  yaw  sas  limit  on  the  actual  aircraft  has 
been  expanded  to  19.0  deg  (Ref.  24)  so  a  digitally-augmented 
aircraft  with  this  larger  yaw  sas  limit  could  probably  match 
the  accuracy  of  the  DFCS  quite  well. 

The  rolling  reversal  trajectory  is  flown  by  the  pilot 
model  without  roll  rate  command- to-pedal  interconnect;  hence, 
the  pilot  does  not  provide  any  pedal  input.  The  rudder  deflec 
tion  shown  in  Fig.  38  is  entirely  due  to  the  digital  augmenta¬ 
tion,  which  essentially  operates  as  an  aileron-rudder  intercon 
nect  in  this  application. 


3 . 7  CHAPTER  SUMMARY 

The  improvements  in  handling  qualities  in  the  ACM 
regime  that  a  DFCS  offers  are  demonstrated  in  this  chapter. 
They  are  accomplished  by  simulating  both  the  ACM-oriented  DFCS 
design  and  a  conventional  aircraft  control  system  with  the 
nonlinear  6-DOF  aircraft  model  and  the  ACM  pilot  model.  The 
test  results  specifically  show  that  the  DFCS  offers  a  more 
accurate,  better-controlled  aircraft  response  along  the  two 
typical  ACM  trajectories  simulated.  In  one  case,  the  DFCS  . 
accurately  controlled  the  aircraft  right  up  to  its  performance 
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boundary,  while  the  conventionally-controlled  aircraft  departed 
from  controlled  flight.  Preliminary  testing  of  a  digital  aug¬ 
mentation  system  indicates  that  it  offers  many  of  the  same 
stability  and  control  improvements  as  the  full  DFCS ,  while 
operating  through  existing  sas  channels.  The  DFCS  design 
methodology  demonstrated  in  this  chapter  is  directly  appli¬ 
cable  to  the  design  of  command  augmentation  systems  for  ad¬ 
vanced  vehicle  configurations  such  as  highly-maneuvering  and 
CCV  aircraft. 
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4. 


PILOT  MODEL  HYPOTHESIS  TESTING 


4 . 1  INTRODUCTION 

During  rapid  maneuvering,  an  aircraft's  dynamic  char¬ 
acteristics  can  change  markedly  in  a  matter  of  seconds.  Hence, 
as  the  aircraft  maneuvers,  the  pilot  may  be  called  upon  to 
change  his  control  strategy  dramatically  in  order  to  maintain 
aircraft  stability.  Usually,  the  pilot  has  mastered  the  neces¬ 
sary  procedural  adaption  needed  to  perform  the  maneuver  and 
executes  it  with  precision.  On  some  occasions,  even  the  skil¬ 
led  pilot  may  experience  difficulty,  perhaps  adapting  his  con¬ 
trol  strategy  to  suit  poorly  chosen  criteria,  or  possibly  not 
adapting  at  all.  In  high-performance  aircraft,  this  apparent 
lapse  can  cause  a  pilot-induced  "departure",  i.e.,  a  loss  of 
control  which,  if  not  corrected  immediately,  can  lead  to  a 
spin  and  possible  loss  of  aircraft.  An  example  of  a  spin  de¬ 
parture  is  shown  in  Fig.  32  for  the  case  where  a  simple  closed- 
loop  system  is  chosen  to  resemble  a  pilot's  outer  loop  control 
behavior;  this  simple  pilot  model  does  not  adapt  as  the  aircraft 
executes  a  wind-up  turn.  The  problem  of  interest  becomes  one 
of  determining  if  more  complex  pilot  models  can  be  used  to 
model  and  understand  the  pilot's  discretionary  behavior  in 
departure-prone  maneuvering  tasks. 

In  two  previous  reports,  Refs.  3  and  23,  the  optimal 
control  pilot  model  has  been  employed  to  address  this  problem 
using  the  F-14A  fighter  aircraft.  The  optimal  control  pilot 
model  is  a  complex  optimal  control  system  which  has  been  shown 
to  capture  the  fundamental  aspects  of  a  human  operator's  con¬ 
trol  behavior  (Refs.  26  and  27).  The  optimal  control  pilot 


model  is  determined  by  minimizing  a  quadratic  cost  function  in 
the  presence  of  noise  subject  to  a  human's  physical  limitations 
and  perceptions.  Section  4.2  reviews  the  comp  ■/  ents  of  the 
optimal  control  pilot  model  used  in  this  study. 

In  Ref.  23,  investigations  are  conducted  using  the 
optimal  control  pilot  model  and  F-14A  along  the  demanding, 
wind-up  turn  maneuver  previously  mentioned.  It  is  determined 
in  Ref.  23  that  if  a  pilot  attempts  to  control  an  unstable 
system  optimally,  there  are  levels  of  system  instability  be¬ 
yond  which  the  pilot  looses  control.  In  many  cases,  a  human 
is  physically  incapable  of  controlling  systems  with  substantial 
instabilities.  Reference  23  also  shows  that  the  optimal  con¬ 
trol  pilot  model  exhibits  similar  difficulties  and  can  predict 
when  a  pilot  will  loose  control  due  to  excessively  severe  air¬ 
craft  instabilities. 

A  sequence  of  flight  conditions  are  specified  along  a 
wind-up  turn  maneuver  for  the  F-14A  in  Ref,  23  and  optimal 
control  pilot  models  are  determiend  at  each  flight  condition. 

The  pilot  models  have  no  difficulty  in  stabilizing  the  aircraft, 
indicating  that  controlling  the  F-14A  along  a  wind-up  turn  is 
within  the  physical  capabilities  of  a  pilo*  if  the  pilot  adapts 
correctly  to  the  flight  condition  of  the  aircraft.  As  the  man¬ 
euver  progresses,  however,  Ref.  23  demonstrates  that  there  is 
a  dramatic  variation  in  the  optimal  piloting  strategy,  including 
in  some  cases,  a  change  in  sign  of  the  pilot's  stabilizing 
commands  to  the  aircraft. 

To  determine  if  the  pilot  model  control  strategy  re¬ 
quired  to  complete  a  maneuver  is  somehow  related  to  known  pilot- 
aircraft  control  difficulties,  the  pilot  model's  control  strategy 
is  fixed  in  Ref.  23  at  a  low  angle  of  attack  flight  condition 
(10  deg)  and  the  aircraft  dynamics  are  varied  in  sideslip,  PQ, 
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and  angle  of  attack,  ao>  along  a  wind-up  turn  maneuver.  Vary¬ 
ing  sideslip  is  shown  to  have  no  effect  on  stability  for  a  non- 
adapted  pilot  model;  but,  when  the  aircraft's  angle  of  attack 
is  increased  beyond  about  16  deg,  closed-loop  pilot-aircraft 
instabilities  occur  in  the  form  of  a  lateral  divergence.  The 
onset  of  the  instability  is  shown  in  Ref.  23  to  coincide  with 
the  change  in  sign  of  the  optimal  adapted  pilot  model  control 
gains.  The  pilot  model  results  give  the  indication  that  if 
the  pilot  uses  a  low  angle  of  attack  control  strategy  in  high 
angle  of  attack  flight  conditions,  pilot  control  difficulties 
will  result. 


Continuing  with  the  concept  of  the  pilot  incorrectly 
adapting  to  the  aircraft's  flight  condition,  Ref.  3  varied  the 
pilot  model's  adaptation  point  over  a  number  of  flight  con¬ 
ditions  in  a  wind-up  turn  maneuver.  Five  important  results 
about  pilot-aircraft  stability  were  obtained:  (1)  If  the  pi- 
lot  is  controlling  the  bare  airframe  using  only  lateral  stick, 
the  pilot-aircraft  system  will  depart  from  controlled  flight 
if  the  pilot  remains  adapted  to  any  low  aQ  (below  12  deg  aQ) 
flight  condition  while  the  aircraft  flight  condition  exceeds 
16  aQ;  (2)  If  the  pilot  attempts  to  compensate  for  control 
difficulties  at  16  deg  a  and  adapts  by  using  a  high  aQ  control 
strategy  in  a  low  aQ  flight  condition,  the  aircraft  again  can 
go  unstable  because  of  improper  pilot  control  behavior;  (3) 

The  pilot  model  analysis  did  show,  however  that  the  pilot  can 
compensate  for  control  difficulties  beyond  16  deg  ofQ  without 
adapting  by  using  lateral  stick  and  pedals;  (4)  Even  if  the 
pilot  uses  stick  and  pedals,  pilot-aircraft  instabilities  still 
occur  beyond  25  deg  aQ  if  the  pilot  model  keeps  a  low  aQ  control 
strategy;  and  (5)  If  the  pilot  uses  only  pedal  to  control  the 


*The  F-14A  stability  augmentation  system  (Sas)  is  off  and  an 
aileron-rudder  interconnent  system  (ARI)  is  not  operating. 
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lateral-directional  dynamics  of  the  aircraft,  then  the  pilot 
model  analysis  in  Ref.  3  shows  the  pilot-aircraft  system  re¬ 
mains  stable  for  any  mismatch  between  the  pilot  model  adapta¬ 
tion  point  and  the  aircraft's  flight  condition  up  to  the  limit 
of  angle  of  attack  variation  tested  (30  deg  <*0)-  Of  course, 
this  does  not  indicate  that  the  aircraft  can  be  flown  in  ACM 
with  only  pedals,  but  only  that  stability  would  not  be  in  ques¬ 
tion  if  only  pedals  are  used. 

All  of  these  results  concerning  pilot  model-aircraft 
stability  agree  with  observed  behavior  of  a  human  pilot  flying 
the  bare  airframe  F-14A  aircraft  based  upon  comparisons  with 
piloted-simulation  results  and  discussions  with  NASA  Langley 
personnel  (Ref.  28). 

The  conclusions  just  described  indicate  that  a  proper¬ 
ly  designed  aileron-rudder  interconnect  system  would  probably 
cure  the  high  a  departure  problem,  and  allow  the  pilot  to  use 
a  low  a  piloting  strategy  at  all  angles  of  attack.  An  initial 
F-14A  ARI  design,  which  essentially  eliminates  pilot  differen¬ 
tial  stabilator  input  at  high  a,  while  phasing  in  a  lateral 
stick- to-rudder  crossfeed,  is  tested  in  Ref.  3  and  23  using 
the  optimal  control  pilot  model.  The  pilot  model  analysis 
shows  that  the  ARI  eliminates  the  lateral  divergence  instabil¬ 
ity  but  introduces  a  new  and  different  instability  character¬ 
ized  as  a  growing  oscillation  in  the  lateral  modes,  particu¬ 
larly  in  the  sideslip  angle.  The  piloted  F-14A  at  high  angle 
of  attack  with  the  ARI  on  has  a  pilot  induced  oscillation  prob¬ 
lem  (Ref.  29)  very  similar  to  the  instability  predicted  by  the 
optimal  control  pilot  model  analysis . 

The  remarkable  ability  of  the  optimal  control  pilot 
model  to  predict  instabilities  of  the  actual  pilot-aircraf t 
system  is  qualitative  rather  than  quantitative.  Time-domain 
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comparisons  of  the  predicted  output  of  the  pilot  model  and  the 
measured  output  of  the  actual  pilot  in  a  tracking  situation 
are  not  performed  in  Refs.  3  and  23.  Agreement  between  the 
pilot  model  and  actual  pilot  time  histories  is  judged  by  examin¬ 
ing  time  histories  of  pilot  control  difficulties  and  determining 
if  the  observed  difficulties  agree  with  the  stability  results 
computed  in  Refs.  3  and  23.  It  is  the  purpose  of  this  chapter 
to  perform  a  rigorous  analysis  of  the  aircraft  control  time  his¬ 
tory  predicted  by  the  optimal  control  pilot  model  and  the  meas¬ 
ured  control  time  history  of  a  pilot  flying  the  F-14A  in  a  wind- 
up  turn  tracking  task. 

The  method  used  to  perform  the  analysis  is  a  hypothe¬ 
sis  testing  scheme  developed  in  Refs.  30  and  31  and  briefly 
outlined  in  Section  A. 3  and  Appendix  C.  Given  a  set  of  hypo¬ 
thesized  mathematical  models  used  to  represent  a  physical  sys¬ 
tem,  the  hypothesis  testing  scheme  can  determine  which  model 
best  matches  the  measured  output  of  the  physical  system.  To 
use  the  hypothesis  testing  approach  with  the  optimal  control 
pilot  model,  a  sequence  of  pilot  models  is  constructed  at  flight 
conditions  along  the  wind-up  turn  maneuver.  The  pilot  models 
are  thereby  adapted  to  the  angle  of  attack  of  the  flight  condi¬ 
tions  used  for  their  construction.  Thus,  there  can  be  a  low-a 

o 

control  strategy  pilot  model  and  a  high-aQ  control  strategy 
pilot  model.  A  pilot  model's  control  strategy  is  determined 
by  the  feedback  gain  the  optimal  control  pilot  model  produces 
at  the  flight  condition. 

Recently,  other  researchers  have  been  investigating 
the  concept  that  the  human's  model  of  the  system  under  control 
and  the  actual  system  model  do  not  necessarily  agree.  The 
human's  model  of  the  system  under  control  is  termed  the  "in¬ 
ternal  model",  (Refs.  32  to  35).  Henceforth,  instead  of  say¬ 
ing  that  the  pilot  model  is  nonadapted,  the  pilot  model  will 
be  said  to  have  an  incorrect  internal  model  of  the  aircraft 
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dynamics.  A  low-<yQ  pilot  model  means  the  internal  model  of 
the  aircraft  used  to  construct  the  pilot  model  is  obtained  at 
a  low-ao  flight  condition. 

The  first  step  in  the  hypothesis  testing  procedure  is 
the  constructon  of  a  finite  set  of  pilot  models,  each  with  a 
different  internal  model  representation  of  the  aircraft.  Then 
the  output  of  each  pilot  model  is  compared  with  the  human  data 
to  determine  which  model  best  predicts  the  actual  pilot  response. 
The  output  of  the  hypothesis  testing  scheme  is  the  probability, 

P.j  ,  that  pilot  model  i  best  matches  the  data.  The  hypothesis 
testing  scheme  will  determine  if  the  pilot  does  not  adapt  to 
varying  aircraft  flight  conditions  by  indicating  only  one  pilot 
model  has  a  high  probability  of  matching  the  data. 

The  hypothesis  testing  scheme  is  very  general  and 
does  not  have  to  be  restricted  to  investigating  only  situa¬ 
tions  where  the  pilot  models  internal  model  of  aircraft  var¬ 
ies.  Section  4.4,  for  example,  validates  the  hypothesis  test¬ 
ing  scheme  using  synthetic  data  by  constructing  a  set  of  pilot 
models  which  differ  from  each  other  in  key  pilot  model  param¬ 
eter  values.  One  of  the  hypothesized  pilot  models  is  used  to 
generate  the  data  in  a  monte  carlo  simulation  and  the  hypo¬ 
thesis  testing  scheme  correctly  identifies  this  pilot  model  as 
the  one  which  best  matches  the  data, 

The  actual  piloted  data  used  in  the  hypothesis  test¬ 
ing  scheme  for  the  F-14A  aircraft  is  obtained  from  NASA  Langley 
aircraft  simulation  tests  using  their  Differential  Maneuvering 
Simulator  (DMS).  The  DMS  consists  of  an  enclosed  hemispherical 
screen  within  which  a  pilot  sits  in  a  cockpit  mock-up.  He 
views  a  computer  generated  display  of  a  tracking  situation  on 
the  screen,  The  pilot  generates  control  commands  which  drive 
the  computer  model  of  the  aircraft.  The  computer  portrays  the 
resulting  aircraft  motions  on  the  screen  as  viewed  by  the  pilot. 
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Section  A. 5  presents  time  histories  of  the  aircraft 
states  and  controls  which  result  during  the  DMS  tracking  en¬ 
gagement  under  investigation.  Section  4.6  presents  the  in¬ 
ternal  model  representation  of  the  tracking  dynamics  used  in 
constructing  the  optimal  control  pilot  models.  Appendix  A 
gives  a  detailed  derivation  of  these  equations  of  motion. 
Section  4.7  discusses  the  variables  which  are  perceived  by  the 
pilot  and  modeled  in  the  pilot  model  as  the  observations. 

Also  given  are  the  typical  human  prameters  (perceptual  time 
delay,  neuromusclar  time  constant,  etc.  discussed  in  Section 
4.3)  used  to  construct  the  optimal  control  pilot  model.  The 
hypothesis  testing  results  using  the  actual  data  are  given  in 
Section  4,8.  The  chapter  is  summarized  in  Section  4.9. 


4.2  OPTIMAL  CONTROL  PILOT  MODEL  WITH  THE  PADE  APPROXIMATION 

This  section  briefly  reviews  the  elements  of  the  opti¬ 
mal  control  pilot  model  to  be  used  in  hypothesis  testing.  The 
optimal  control  pilot  model  is  based  on  the  premise  that  a 
motivated,  well-trained  human  controls  a  system  optimally.  To 
construct  the  pilot  model,  the  assumptions  needed  to  specify 
an  optimal  controller  are  formulated,  then  modified  to  reflect 
basic  human  characteristics  and  limitations.  The  pilot  model 
assumptions  are  shown  in  Table  12.  At  chosen  points  along  the 
wind-up  turn  maneuver  under  investigation,  the  pilot's  internal 
model  of  the  aircraft  and  target  dynamics  is  represented  as  a 
linear,  time-invariant  system.  The  n-vector,  Ax(t)  ,  represents 
the  perturbation  aircraft  dynamics,  aircraft  sas,  tracking 
error,  and  target  states.  The  aircraft's  stick  inputs  are 
represented  by  the  m-vector,  Au(t);  FAw(t)  is  the  pilot's  inter¬ 
nal  model  of  white  gaussian  noise  disturbance  inputs. 


TABLE  12 

PILOT  MODEL  ASSUMPTIONS 


T-1063 


AIRCRAFT  AND  TARGET 
PERTURBATION  STATE 
DYNAMICS 

PILOT  OBSERVATIONS 


PILOT  COST  FUNCTION 


PILOT  NEUROMUSCULAR 
DYNAMICS 


PILOT  NEUROMUSCULAR 
LAG 


Ax(t)  ■  F  AX( t )  T  G  Au(t)  +  r.  t) 


4£(t)  -  fH  D] 


'Ax(t-x)' 

_AU(t-T)_ 


+  AV..(  t-T  ' 


iim 

T**-« 


i  h:3 


+  Au  RcAu]dt 


Au(t)  ■  -RLAu(t)  +  AUg ( t )  +  Av^(t) 
0  ...  0 


‘  1 


The  pilot  is  assumed  to  manipulate  the  aircraft  con¬ 
trols  so  as  to  minimize  the  quadratic  cost  function  shown  in 
Table  12.  The  pilot  minimizes  the  cost  function  using  his  per¬ 
ceptions  of  displayed  information  represented  as  the  jK-obser- 
vation  vector,  Ay(t).  As  a  physical  limitation,  the  pilot  has 
a  perceptual  time  delay,  x,  which  delays  his  observations. 

The  pilot  cannot  physically  observe  the  states  perfectly,  hence, 
a  white,  gaussian  observation  noise,  Av^(t)  is  added  to  the 
observations.  The  pilot  manipulates  the  aircraft  controls 
with  his  hand,  hence,  a  neuromuscular  dynamics  model  must  be 
included.  Weighting  the  control  rate  in  the  quadratic  cost 
function  causes  the  control  solution  to  take  the  form  required 
to  model  neuromuscular  dynamics.  The  (mxm)  neuromuscular 
dynamics  matrix,  R^,  in  Table  12  is  diagonal,  with  individual 
elements  representing  the  inverse  of  human  limb  neuromotor 
time  constants.  The  neuromotor  dynamics  are  driven  by  the 
pilot's  internal  control  commands,  Au  (t),  and  by  neuromotor 
noise  represented  as  the  white  gaussian  m-vector,  Av  (t). 
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As  discussed  in  Ref.  27,  the  pilot's  observation  noise 

covariance  scales  with  the  covariance  of  the  states  and  controls, 

while  the  neuromotor  noise  scales  with  the  covariance  of  the 

pilot's  internal  control  commands.  The  covariances  for  Avy(t) 

and  Avu(t)  can  be  determined  if  the  scaling  factors  Pyi  and 

P  .  are  specified,  where  P„ ,  is  the  pilot's  noise- to-signal 
^  ^  t  h  y 

ratio  for  the  i  observation  and  P  . ,  is  the  pilot's  noise- 
f  h  m 

ratio  for  the  i  control. 

To  determine  an  optimal  control  pilot  model,  the  fol¬ 
lowing  parameters  discussed  in  the  previous  paragraphs  and 
shewn  in  Table  12  must  be  known  for  a  given  aircraft  system: 

Pyi ,  Pu£ ,  x,  R^,  H,  D,  and  Qc .  All  of  these  variables  except 
Qc ,  H,  and  D  have  been  measured  experimentally,  and  typical 
ranges  of  their  values  can  be  found  in  the  pilot  model  litera¬ 
ture  (P^  =  -20  db ,  Pu^  =  -30  db,  x=0.2  sec,  xn  =  0.1  sec). 

When  the  pilot  views  motions  in  a  display,  as  is  the  case  for 
the  DMS  simulations,  then  the  elements  in  H  and  D  can  be  de¬ 
duced  from  the  states  portrayed  on  the  display  screen.  The 
values  used  for  the  state  and  control  quadratic  weights  in  Qc 
remain  the  only  unknown  parameters  needed  to  construct  an  op¬ 
timal  control  pilot  model.  These  weights  model  the  pilot's 
trade-off  of  control  effort  against  tracking  accuracy.  Three 
previous  studies,  Ref.  23,  36,  and  37  aid  in  choosing  Qc .  In 
Ref.  23,  the  pilot  model's  control  strategy  (i.e.,  pilot  model- 
aircraft  closed-loop  eigenvalues)  is  shown  not  to  be  especially 
sensitive  to  the  choice  of  Qc .  As  long  as  reasonable  values 
are  chosen  for  the  weighting  elements,  conclusions  about  pilot 
control  strategy  can  be  made  with  confidence.  In  Refs.  36  and 
37,  a  tracking  situation  similar  to  the  one  developed  in  Section 
4.5  is  analyzed  with  the  optimal  control  pilot  model.  Choices 
for  the  weighting  elements  in  Qc  are  given  in  Refs.  36  and  37 
and  are  used  as  an  aid  in  constructing  the  pilot  model  here. 
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It  is  also  possible  to  use  the  hypothesis  testing  scheme  to 
determine  which  quadratic  cost  function  weighting  matrix  (track¬ 
ing  accuracy-control  effort  trade  off)  among  a  given  set  of 
weighting  matrices  best  matches  the  actual  pilot's  control 
strategy.  Section  A. 5  shows  hypothesis  testing  results  for 
different  cost  function  weighting  matrices  and  serves  to  veri¬ 
fy  that  the  tracking  accuracy-control  effort  trade-off  assumed 
in  this  study  closely  approximates  that  actually  practiced  by 
the  pilot. 

The  hypothesis  testing  scheme  requires  an  analytic 
pilot  model,  hence,  the  pure  time  delay  in  the  model  is  re¬ 
placed  by  its  Pade  approximation.  The  Pade  approximation  al¬ 
lows  the  pilot  model  to  be  constructed  as  a  linear,  time- 
invariant  model  with  the  pure  delay  approximated  by  additional 
pilot  model  states.  The  solution  to  the  optimal  control  prob¬ 
lem  posed  in  Table  12  with  the  pure  time  delay  replaced  by  its 
Pade  approximation  is  given  in  Ref.  23  and  summarized  in  Table 
13.  A  block  diagram  is  shown  in  Fig.  39.  The  pilot'  n  observa¬ 
tions  (2  in  Table  13)  are  degraded  by  noise,  then  passed  through 
a  lead-lag  network  representing  the  Pade  approximation.  The 
resulting  signal  is  processed  by  a  Kalman  filter  (3  in  Table 
13)  which  generates  a  best  estimate  of  the  states,  controls, 
and  lagged  observation  states.  The  state  estimates  are  multi¬ 
plied  by  the  feedback  matrix,  C,  (4  in  Table  13)  to  form  the 
pilot's  internal  control  command.  The  gain  matrix,  C,  in  the 
pure  time  delay  and  Pade  approximation  pilot  models  are  the 
same.  The  pilot's  internal  control  command,  Auc(t)  is  the 
input  to  the  neuromuscular  dynamics  model  (5  in  Table  13) 
from  which  the  pilot's  aircraft  control,  Au(t)  is  generated 
and  sent  to  the  aircraft  model's  control  actuators. 

With  the  Pade  approximation,  the  pilot  model  described 
in  Table  13  and  Fig.  39  can  be  written  as 
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TABLE  13 

ELEMENTS  OF  THE  PILOT-AIRCRAFT  SYSTEM  MODEL 


EQUATION 

EQUATION  PARAMETERS 

RELATION  TO  PILOT  PERFORMANCE 

1.  Aircraft  Dynamic  Model  In¬ 
cluding  Internal  Target  Model 
of  Pilot  and  Tracking  error 
dynamics: 

AX 

Aircraft,  target,  and 
tracking  error  atate 
variables 

Pilot  Bust  obaerve  thi»  well 
enough  to  command  aircraft  and 
and  to  provide  stability. 

dx  =  Fax  +  CAu 

Au 

Aircraft  control  vari¬ 
ables 

Pilot  must  use  this  to  coaaand 
aircraft  and  to  provide  stabil¬ 
ity. 

F 

Aircraft,  target,  and 
tracking  error  dynamics 
(stability  derivatives 
and  inertial  coupling) 

Aircraft  must  be  stable  enough 
for  pilot  to  control,  subject 
to  normal  human  capabilities 

G 

Aircraft  control  effects 
(sensitivity  to  control 
deflections) 

Aircraft  must  respond  to  exter¬ 
nal  commands  in  a  way  which  the 
pilot  can  understand 

2.  Pilot's  Cues: 

Ay  *  HAx  ♦  DAu 

H 

State  variable  display 
selection  and  transforma¬ 
tion 

Cues  must  be  sufficient  for 
command  and  stabilization. 

Pilot's  Cues  Delayed  by,  r, 
using  a  Pade  approximation 

t 

Pilot  time  delay 

A*  =  -  |  ♦  AJ  +  AVy 

D 

Control  variable  display 
selection  and  transforma¬ 
tion 

A 

AyE  *  -  Ay  -  Ayy  ♦  -  Az 

Pilot  induced  zero  mean 
gaussian  noise  in  observa¬ 
tions 

Noise  in  observation  has  di¬ 
rect  effect  on  estimation  per¬ 
formance  of  the  pilot 

3.  Pilot  Estimation  Model  with 
Pad^  Approximation: 

*-E 

State.  Ax,  control,  Av 
and  perception-delay 
state.  Az,  variables 
estimated  by  pilot 

Estimates  of  motion  variables 
must  be  accurate  enough  to  pro¬ 
vide  effective  closed-loop  con¬ 
trol 

4e  1  F^E  *  K(A£e-HeA«> 

F  G  0 

fe  =  c  -rl  0 

H  D  [ -|l ) 

fe 

K 

Dynamic  model  assumed 
(i.e.,  "learned")  by  the 
pilot  including  neuromus¬ 
cular  lags  and  perception 
time  delays 

Estimation  gains  which 
weight  the  difference  be¬ 
tween  the  pilot’s  obser¬ 
vations  and  his  predic¬ 
tion  of  pilot-aircraft 
response . 

The  better  the  pilot's  know¬ 
ledge  of  the  aircraft  and  his 
own  capabilities,  the  better  he 
can  cope  with  noisy  measure¬ 
ments 

Less  noise  in  the  pilot's  ob¬ 
servation  of  cues  leads  to 
high  K  and  more  reliance  on 
observed  motions. 

^E 

Pilot’s  delayed  observa¬ 
tion  of  motion,  control, 
delay  states 

Noise  in  observations  has  di¬ 
rect  effect  on  estimation  per¬ 
formance  of  the  pilot. 

H,  =  [-H  -D  |*  T 1 ] 

he 

Pilot’s  t ransformat ion  of 
estimated  variables  to 
agree  dimensionally  with 
observations 

Pilot  disorientation  would  de¬ 
grade  performance. 

U.  pilot  Control  Model: 

*HC  =  c*x  4  Ayu 

A& 

Subset  of  ANp  correspond¬ 
ing  to  aircraft  target 
and  tracking  error  mo¬ 
tions 

Ax  must  be  close  to  Ax  for  pre¬ 
cise,  stable  control." 

c 

Control  gains  which  trans¬ 
form  pilot's  estimates  of 
aircraft  motions  to  con¬ 
trol  actions 

Pilot  attempts  to  tradeoff  air¬ 
craft  motions  and  available 
control  "power".  Improper  con¬ 
trol  strategy  could  degrade 
command  response  and  destabil¬ 
ize  the  system. 

“'-U 

Pilot  induced  zero  mean 
neuromotor  gaussian 
noise . 

Reflects  the  fact  that  the 
pilot  cannot  estimate  some  of 
the  stales  perfectly. 

5.  Pilot  Neuromuscular  Model: 

Au 

Pilot  mode]  control 
output  variables 

Aircraft  control  input 

Au  -  -R.  Au  *  Au„ 

L  -  -c 

Neuromuscular  lags 

Neuromuscular  system  smooths 
pilot  outputs  and  could  prevent 
pilot  from  stabilizing  a  fast 
instability. 

*ot 

Pilot  internal  control 
commands 

Result  of  conscious  effort  to 
provide  "best"  control, 
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*PM  =  e 


FPMAt 


At 


F„ut 


PM  =  X  6  PM  dtKPM 


e  *JT 

AX  =  [AUT 


??ht 


^PM 


f™1  at 


=  W 


PM 


A^£ ,  aJ]T 


COO  0 

EE-mE  7K 

o  -|k 


Vu  0  0 

0  KVyKT  -KVy 

0  -VyKT  Vy 

[0  -,T  ,»]' 


Vu,  Vy  are  the  covariance  matrices  of  the  pilot's  neuromotor 
noise  Avy  and  observation  noise,  Avy,  respectively.  The  mat¬ 
rix,  K,  is  the  pilot  model's  Kalman  filter  gain.  The  matrices 
K  and  C  are  determined  by  solving  Riccati  equations  which  de¬ 
pend  on  the  pilot's  internal  model  description  of  the  aircraft 
(1  in  Table  13).  The  pilot  model  is  unadapted  when  the  internal 
model  description  of  the  aircraft  used  to  determine  the  pilot, 
model  is  not  the  same  as  the  mathematical  model  of  the  aircraft 
at  the  flight  condition  under  investigation. 


.103 


HYPOTHESIS  TESTING 


4.3 


The  hypothesis  testing  philosophy  for  investigating 
pilot  control  strategy  is  to  construct  a  number  of  pilot  mod¬ 
els,  each  using  an  internal  model  of  the  aircraft  at  different 
points  along  the  wind-up  turn  maneuver;  then,  determine  which 
pilot  model  best  represents  the  actual  pilot  control  behavior. 

If  the  pilot  does  not  adapt  as  the  aircraft  flight  condition 
changes,  then  only  one  of  the  pilot  models  will  best  represent 
the  data.  If  the  pilot  adapts  then  the  hypothesis  testing  scheme 
will  switch  among  the  pilot  models  as  the  pilot's  control  stra¬ 
tegy  changes. 

Mathematically,  the  hypothesis  testing  scheme  is  im¬ 
plemented  by  viewing  each  pilot  model  represented  in  Eq.  23  as 
a  dynamic  system  whose  input  is  the  pilot  observation  vector, 
with  process  noise,  Awk,  and  whose  output  is  the  pilot 
control  vector  Au^.  The  measured  output  vector  is  the  actual 
measured  pilot  control, 

AVk  =  11 0  01  4*k  +  *k 

is  a  zero-mean  gaussian  measurement  noise  with  covariance 
V.  In  the  data  from  NASA  Langley,  the  DMS  measurement  noise 
is  essentially  zero,  but  V  is  included  here  to  keep  the  deri¬ 
vation  general,  A  Kalman  filter  using  the  pilot  model  as  the 
plant  and  the  control  measurements  as  the  observation  vector 
can  be  constructed  as 

A$Jk(  +  )  =  Akk(-)  +  Kh  [Au^  -  Auk(-)]  (24) 

4Sk*i(-)  *  *PM  4Sk  *  rpM 


(25) 


"I" 

'  I ' 

kh  •  P(-) 

0 

{[I  0  0]P(-) 

0 

.0. 

.0. 

+  V} 


P(-)  =  4>PMP(  +  J4>PM  *  WPM 
P(  +  )  =  1 1  -  KH  1 1  0  0]]  P(-) 

In  this  application,  the  above  Kalman  filter  is  assumed  to  be 
in  steady  state,  If  the  pilot  model  is  correct,  the  filter 
residuals 


are  a  zero-mean  white  gaussian  noise  sequence  with  covariance 


E{ArkAr£}  = 


I  0  0]  P(-) 


+  V  s  S 


(26) 


If  the  pilot  model  is  not  correct,  then  Ark  may  not  be  a  white 
gaussian  noise  sequence  and  does  not  satisfy  Eq.  26. 

In  the  hypothesis  testing  scheme,  the  different  resi¬ 
duals,  Ar^ ,  for  each  hypothesized  pilot  model  are  compared 
against  each  other  to  determine  which  prediction,  A&k(-),  best 
causes  Ark  to  be  a  white  gaussian  noise  sequence.  The  mathe¬ 
matical  details  of  this  comparison  process  are  presented  in 
Appendix  C  and  are  based  on  work  in  Refs.  30  and  31.  The  final 
result  is  that  a  rt-cursive  equation  can  be  constructed  for 
computing  the  probability,  ^(k),  that  pilot  model  i  best 
fits  the  measured  data.  Given  N  pilot  models  there  are  N  pro¬ 
babilities,  P„  . (k),  and  the  sum  of  the  probabilities  must  be 

n  )  1 

one.  A  modification  is  made  to  the  hypothesis  testing  scheme 
to  insure  that  no  pilot  model  probability  can  go  below  e  ,  which 
in  this  case  is  taken  to  be  0.05.  This  insures  that  if  the 
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pilot  switches  control  strategy,  the  modified  scheme  can  change 
its  probabilities  to  reflect  the  switch. 


Care  must  be  taken  in  applying  the  Kalman  filter  shown 
in  Eqs.  24  and  25  to  a  situation  where  measured  pilot  data 
comes  from  an  actual  nonlinear  pilot  flight  test  or  simulation 
as  is  the  case  here.  The  filter  requires  perturbation  varia¬ 
bles,  Aum  and,  4^,  but  only  total  values,  u^  ^ ,  and, 
are  available  from  the  flight  test.  A  common  procedure  in 
classical  control  designs,  Ref.  38,  is  used  to  rectify  the 
situation.  If  necessary,  the  measured  pilot  controls  and  ob¬ 
servations  will  first  be  passed  through  a  high-pass  filter 
(sometimes  known  as  a  wash-out  filter)  to  separate  the  low 
frequency  nominal  states  from  the  high  frequency  perturbation 
states.  The  high-pass  filter  is  converted  to  discrete- time 
using  the  bilinear  transform,  Ref.  39,  in  order  to  preserve 
its  frequency  domain  characteristics.  A  value  of  4  sec  is 
used  for  the  high-pass  filter  time  constant. 


4.4  VERIFICATION  WITH  SYNTHETIC  DATA 

This  section  validates  the  pilot  model  hypothesis  test¬ 
ing  scheme  described  in  the  previous  section  with  synthetic  data. 
The  synthetic  data  is  created  by  a  monte-carlo  simulation  of  a 
pilot  model  for  a  low-order  plant.  The  precise  pilot  model 
parameters  used  to  create  the  data  are  known;  hence  it  can  be 
determined  if  the  scheme  performs  correctly.  Each  test  has 
three  hypothesized  pilot  models.  One  pilot  model  is  the  known 
correct  model.  The  other  two  pilot  models  differ  from  the 
known  pilot  model  by  varying  a  key  pilot  model  parameter. 

The  low-order  plant,  Ref.  23,  used  to  construct  the 
pilot  models  has  dynamics 
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The  closed-loop  pilot  model  feedback  gain  for  this  low-order 
plant  is  independent  of  in  the  cost  function.  The  pilot 
model  parameters  are  a  time  delay,  t,  of  0.1  sec,  a  neuromotor 
time  constant,  r  ,  of  0.1  sec,  observation  noise-to-signal 
ratios,  P  of  "20  db  and  a  neurometor  noise-to-signal  ratio, 
Pu^  ,  of  -30  db.  The  state  Ax^  and  its  rate  Ax^  are  assumed  to 
be  observed  by  the  pilot  model. 

Four  tests  are  performed.  In  all  four  tests,  the 
initial  condition  on  the  state,  Ax-^O),  is  one  and  all  noise 
sources  are  operational.  In  the  first  test,  two  of  the  hypo¬ 
thesized  pilot  models  differ  from  the  known  pilot  model  in  the 
choice  of  the  feedback  gain,  C.  The  probability  time  histor¬ 
ies  are  shown  in  Fig.  AO.  The  pilot  model  which  has  the  in¬ 
correct  low  feedback  gain  is  initially  chosen,  then  the  pilot 
model  with  the  incorrect  high  feedback  gain  is  chosen,  then 
the  true  pilot  model  gradually  increases  in  probability.  Fig¬ 
ure  AO  is  an  example  of  what  happens  in  hypothesis  testing 
when  the  different  models  are  very  similar.  The  incorrect 
initial  choices  are  believed  to  be  caused  by  the  state  initial 
condition  which  asymptotically  damps  out  after  3  seconds. 

The  second  test  uses  three  pilot  models  where  two  of 
the  pilot  models  have  incorrect  values  for  the  observation 
time  delay.  The  pilot  model  which  uses  the  correct  time  delay 
eventually  has  the  highest  probability  as  shown  in  Fig.  Al. 
Note  that  visual  observation  of  the  predicted,  Uj<(-),  and  ac¬ 
tual  pilot  model  controls,  u  .  ,  in  Figs.  AO  b,  c,  and  d  can- 

in  i  k 

not  distinguish  which  hypothesized  pilot  model  performs  best. 
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Figure  40  Pilot  Model  Hypothesis  Testing  Results  For 
Different  Feedback  Gains 

The  third  and  fourth  tests  are  shown  in  Figs.  42  and 
43.  In  both  cases,  three  neuromotor  noise  covariances  are 
used  to  form  three  hypothesized  pilot  models.  In  both  cases 
the  pilot  model  with  the  correct  neuromotor  noise  covariance 
reaches  the  highest  probability.  In  Fig.  42,  the  pilot  model's 
control  and  observation  synthetic  data  is  high-pass  filtered 
as  discussed  in  the  previous  section  before  the  probability 
sequences  are  computed.  Figures  41a  and  42a  are  virtually 
Identical  indicating  that  the  high-pass  filtering,  if  neces¬ 
sary,  should  not  greatly  affect  the  results. 

This  section  shows  that  the  pilot  model  hypothesis 
testing  scheme  discussed  in  Section  4.3  can  choose  the  pilot- 
model  which  gives  the  best  fit  to  the  data.  The  next  section 
presents  the  data  from  the  flight  test  that  is  used  for  the 
analysis  in  Secton  4.8 
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Figure  42  Pilot  Model  Hypothesis  Testing  Results 
For  Different  Neuromotor  Noise 
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A. 5 


FLIGHT  TEST  SELECTION 


The  simulated  flight  test  selected  for  study  is  a 
wind  up  turn  tracking  task  using  the  NASA-Langley  Differential 
Maneuvering  Simulator  (DMS).  The  pilot  is  instructed  to  track 
a  target  aircraft  which  is  displayed  on  a  large  hemi-spherical 
screen,  within  which  the  cockpit  is  mounted.  A  sky-earth  re¬ 
presentation  is  also  displayed.  The  target  aircraft's  trajec¬ 
tory  is  prestored.  There  are  no  intentional  disturbances  af¬ 
fecting  either  the  target  aircraft  or  the  piloted  aircraft, 
although  the  prestored  target  trajectory  exhibits  rapid  and 
continuous  maneuvering. 

The  pilot  is  instructed  not  to  use  rudder  pedals. 
Lateral-directional  control  by  the  pilot  is  accomplished  using 
only  lateral  stick.  The  simulation  is  performed  with  the  pitch 
and  yaw  sas  operational,  although  the  roll  cas  is  off.  The 
yaw  sas  improves  the  lateral-directional  stability  of  the  air¬ 
craft  at  high  angles  of  attack.  This  means  that  the  stability 
boundaries  of  the  bare-aircraf t  determined  in  Refs.  1  and  23 
cannot  be  used  directly  to  evaluate  pilot  performance  and  the 
performance  of  the  hypothesis  testing  scheme.  The  sas  states 
also  have  to  be  included  in  the  aircraft  model  increasing  the 
dimension  of  the  problem. 

The  view  through  the  cockpit  wnich  is  useful  in  con¬ 
structing  the  tracking  error  equations  is  shown  in  Fig.  44 . 

The  vertical  and  lateral  components  of  tracking  error  are  en 
and  £y,  respectively,  which  the  pilot  can  directly  nerceive. 

The  difference  between  the  pilot's  fixed  pipper  and  the  air¬ 
craft’s  velocity  vector  are  the  angle  of  attack,  ? ,  and  side¬ 
slip  angle,  p,  which  the  pilot  cannot  directly  perceive.  The 
difference  between  the  aircraft's  stability-axis  roll  angle, 

4>v,  and  the  stability-axis  roll  angle  of  the  target,  is 
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b)  View  Through  Cockpit 
Figure  44  Pursuit-Tracking  Task 
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6<pv  and  is  directly  perceived  by  the  pilot.  The  pilot  does  not 
have  any  motion  cues.  Any  computational  time  delay  in  the  DMS 
between  the  pilot's  stick  movement  and  changes  on  the  screen 
is  assumed  to  be  small  and  is  not  modeled  in  the  analysis. 

A  complete  record  of  the  aircraft  state  and  control 
histories  is  recorded  at  32  samples  per  second.  Additionally, 
the  tracking  angles  and  range  to  the  target  are  recorded,  and 
target  motion  is  inferred  from  these  measurements.  This  ap¬ 
proach  dictated  the  target  aircraft  identification  procedure 
discussed  in  Section  4.6  and  required  the  tracking  equation 
simplifications  shown  in  Appendix  B. 

The  wind-up  turn  tracking  task  lasted  90  sec.  At  50 
sec  into  the  run,  the  pilot  is  forced  to  saturate  lateral  stick 
almost  continuously  in  order  to  continue  tracking.  Because  of 
this,  only  the  first  40  sec  of  the  run  are  analyzed.  Original¬ 
ly  it  was  intended  to  perform  hypothesis  testing  on  both  the 
lateral-directonal  and  longitudinal  dynamics,  coupled  and  un¬ 
coupled.  The  '.omplexity  of  the  sas  states  increases  the  dimen¬ 
sion  of  the  coupled  Kalman  filter  in  Eq .  24  from  26  to  38  states 
placing  the  computation  requirements  of  a  coupled  analysis  out¬ 
side  the  scope  of  this  analysis,  Only  the  uncoupled  lateral- 
directional  dynamics  are  to  be  investigated  with  the  pilot 
model . 

Variations  of  angle  of  attack  and  lateral  tracking 
error,  e^,  are  shown  in  Figs.  45a  and  45b  for  the  40  sec  in¬ 
terval.  The  aircraft  remains  at  a  low  angle  of  attack  flight 
condition  for  the  first  10  seconds,  then  the  aircraft  angle  of 
attack  gradually  increases  passing  through  the  n  =  16  deg  sta¬ 
bility  boundary  discussed  in  the  introduction  of  this  chapter 
without  incident.  As  previously  mentioned,  the  lack  of  in¬ 
stability  at  high  angles  of  attack  is  due  to  the  yaw  sas  since 
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Figure  45  Aircraft  Time  Histories  (Continued) 
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former  DMS  simulations  with  the  bare  aircraft,  at  high  of,  using 
only  lateral  stick,  Ref.  23,  exhibit  piloting  difficulties  of 
beyond  16  deg  a . 

The  trajectory  of  the  target  aircraft  as  seen  through 
the  windscreen  is  shown  Fig.  45c,  The  target  aircraft  starts 
at  the  bottom  of  the  screen  and  is  eventually  brought  to  the 
pipper  area  by  the  end  of  40  sec,  Figure  45d  shows  the  pilot's 
stabilator  commands  from  the  longitudinal  stick.  Although 
there  are  no  disturbances  in  the  system,  the  pilot's  control 
movements  exhibit  considerable  activity.  The  intended  opera¬ 
tion  of  the  high-pass  filter  is  effectively  shown  in  Fig.  45d. 
The  high-pass  filter  moves  the  pilots  stabilizing  perturbation 
actions  back  to  the  zero  axis  from  the  a  nominal  6  =  -2  deg. 

S 

The  aileron  control  command  of  the  pilot  lateral  stick 

is  shown  in  Fig.  45e.  In  this  case,  the  pilot's  command  is 

centered  near  the  zero  axis  and  there  is  no  need  to  use  the 

high-pass  filter,  The  pilot  similarly  manages  to  keep  near 

the  zero  axis  in  Fig.  45b  and,  as  shown  in  Fig.  46,  also  keeps 

6pv  near  the  zero  axis.  Since  only  the  lateral-directional 

axis  is  to  be  investigated  it  is  decided  that  the  high-pass 

filter  is  not  necessary  and  6a,  e  .  and  6<)>  can  be  used  di- 

ay  v 

rectly  in  Eqs.  24  and  25,  If  the  longitudinal  axis  would  have 
been  investigated,  the  high-pass  filter  would  be  required. 

The  pilot  manages  to  keep  sideslip  angle  close  to 
zero  as  shown  in  Fig.  45f.  A  small  sideslip  angle  is  an  in¬ 
dication  the  pilot  is  purposefully  flying  a  coordinated  turn 
and  regulating  sideslip.  Some  of  the  instabilities  at  high 
angle  of  attack  simulated  in  Ref,  3  are  characterized  as  a 
high  frequency  growing  oscillation  in  p  which  is  not  evident 
in  Fig.  45f,  NASA  Langley  provided  a  second  wind-up  turn 
tracking  task  which  had  the  unstable  sideslip  oscillation  time 
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history  but  unfortunately  could  not  be  analyzed  within  avail¬ 
able  resources. 


Using  the  aircraft  time  histories  shown  in  Fig.  45, 
the  pilot  model  hypothesis  test  can  be  well  posed;  when  the 
angle  of  attack  of  the  aircraft  increases  beyond  5  deg  (see 
Fig.  45a),  significantly  changing  the  aircraft's  dynamics, 
does  the  pilot's  control  behavior  (Fig.  45e)  remain  fixed  to  a 
5  deg  a  strategy  or  does  the  pilot  change  his  strategy  to  match 
the  aircraft's  changing  dynamics?  To  begin  answering  this  ques¬ 
tion,  the  hypothesized  dynamical  model  of  the  aircraft-target 
dynamics  must  first  be  developed. 


4.6  TRACKING  TASK  MODEL 


The  tracking  task  model  ((1)  in  Table  13)  is  composed 
of  four  components  at  the  adaptation  point;  the  linearized 
dynamics  o  the  subject  aircraft,  the  sas  dynamics,  the  lin¬ 
earized  tracking  error  dynamics,  and  the  linearized  target 
model.  A  derivation  of  the  first  two  models  is  presented  in 
Appendix  A  while  the  latter  two  models  are  derived  in  Appendix 
B.  A  summary  of  all  the  linearized  models  is  given  in  Table  14. 


Combining  all  the  models  in  Table  14  results  in  the 
tracking  task  mode]  used  here  and  is  shown  in  Eq .  27.  The 
states  Av,  Ar,  Ap,  and  A<£v  in  the  model  are  the  perturbation 
body  y-axis  velocity,  body-axis  yaw  rate,  body-axis  roll  rate, 
and  stability  axis  roll  angle  of  the  aircraft,  respectively. 

The  state,  Ae^,  is  the  perturbation  lateral  tracking  error.  The 
stains  Apw-j<  and  A0V^,  are  the  target's  perturbation  stability-axis 


The  variable  AXgAS  is 


roll  rate  and  roll  angle,  respectively, 
the  state  associated  with  the  wash  out  filter  in  the  yaw  sas 
channel  of  the  aircraft. 
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The  model  is  typical  of  air-to-air  combat  investiga¬ 
tion  models  first  performed  in  Ref.  36  and  later  in  Ref.  37 
using  the  optimal  control  pilot  model.  In  both  references  the 
aircraft  models  are  simplified  linear  time-invariant  longi¬ 
tudinal  dynamics  implemented  on  simulators.  Simple  target 
dynamics  are  used  both  in  the  simulation  and  the  model.  In 
our  case  use  of  the  complete  target  aircraft  model,  as  used  in 
the  DMS,  in  the  tracking  task  model  is  infeasible.  An  ap¬ 
proximation  whereby  only  the  target  aircraft  states  Ap^  and 
A((iv^  are  modeled  in  the  tracking  task  model  must  be  made. 

The  approximation  is  made  by  determining  what  target 
dynamics  are  essential  to  the  pilot  in  formulating  his  inter¬ 
nal  model  of  the  target.  Since  the  stability  axis  roll  angle 
is  the  important  target  state,  only  that  state  and  its  deriva¬ 
tive  are  included  in  the  target  model.  In  the  pilot's  internal 
target  model,  unmodeled  target  model  dynamics  are  lumped  into 
a  zero  mean  gaussian  noise  term  Aw,p  which  drives  ApwT-  The 
implication  of  this  assumption  is  that  the  human  has  no  expli¬ 
cit  information  regarding  Aw^,  and  must  model  its  effect  as 
white  noise.  A  detailed  discussion  of  this  type  of  white  noise 
assumption  is  given  in  Ref.  37.  The  parameters  in  the  simplified 
target  model  in  Table  14  can  be  determined  by  identification 
based  on  measured  values  of  tt>v-p*  The  approach  taken  here  is 
to  perform  an  exhaustive  search  over  ^  anc*  t*ie  covariance 
of  AWj  in  Table  14  using  the  hypothesis  testing  scheme  in  Ap¬ 
pendix  C.  <t>vT  is  not  available  as  a  direct  measurement  from 
the  NASA  DMS  data.  Using  the  available  measurements  and  the 
en,  Cy  dynamics  shown  in  Eqs.  B-12  and  B-13,  is  recon¬ 
structed  from  the  data  as 


50, 


tan 


-1  -V(Ey»P  >■*■«; 
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(28) 
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Back  differencing  is  used  to  calculate  the  derivatives  in  Eq. 

28.  Figure  46  shows  plots  of  the  resulting  calculations.  The 
movements  in  the  aircraft  stability-axis  roll  angle,  $v,  fol¬ 
lows  the  movements  of  6vt  as  the  pilot  tracks  the  target  air¬ 
craft.  The  target's  normal  acceleration  is  also  calculated 

A 

and  shown  in  Fig.  46.  Using  4>v^  in  Fig.  46,  a  series  of  tar¬ 
get  model  parameters  are  postulated,  and  hypothesis  testing  is 

used  to  choose  the  set  which  best  matches  the  actual  target 

2 

data.  This  search  yielded  a^  =  -  4.4,  &2  ~  "  and  E{w,p}  = 

4.0  (deg^/sec^ ) for  the  pilot's  internal  target  model  parameters. 

4.7  PURSUIT  TRACKING  PILOT  MODEL  PARAMETERS 

The  parameters  left  to  be  specified  for  the  pilot  model 
are  the  pilot  cues  (2  in  Table  13),  a  preliminary  choice  of 
the  pilot's  quadratic  state  weighting  matrix,  Q. ,  in  Table  12, 

w 

the  standard  human  parameters  discussed  in  Section  4.3,  and 
the  flight  conditions.  The  flight  conditions  are  a  priori 
specifications  of  hypothesized  aircraft  dynamics  at  which  the 
pilot  may  formulate  a  control  strategy,  The  perturbation  dy¬ 
namics  at  the  chosen  flight  conditions  form  the  pilot's  inter¬ 
nal  model  of  the  aircraft. 

From  Fig.  44,  the  pilot  cues  are  the  lateral  error,  Ae^, 
and  stabili ty-axis  roll  angle  error,  AS<J>v.  As  is  customary  in 
optimal  control  control  pilot  models,  the  rates  of  these  states 
are  also  assumed  observed  by  the  pilot  resulting  in 

AyT  =  iAcy,  A cy>  A6<t»v,  A&4>v ] 

The  derivatives  of  the  states  in  Ay  are  computed  by  back  dif¬ 
ferencing  the  measurements  Ae^T  and  A54>v.  From  Fig.  45a,  three 
flight  conditions  which  primarily  differ  in  angle  of  attack 
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are  chosen  for  the  pilot's  internal  model;  low,  aQ  -  5.2  deg, 
medium,  aQ  -  9.3  deg,  and  high,  ofQ  =  19.2  deg. 

The  quadratic  cost  function  weights  quantify  the  pilot's 
tradeoff  of  tracking  error  and  control  effort.  High  weights 
imply  that  the  pilot  expends  significant  control  effort  in 
order  to  minimize  tracking  error.  He  can  be  said  to  be  a  "high- 
gain"  or  "tight"  pilot.  On  the  other  hand,  low  quadratic  weights 
will  typify  a  "low-gain"  or  "loose"  piloting  technique.  The 
specific  values  for  the  quadratic  weights  used  in  the  tests 
reported  here  are  given  in  Table  15.  The  tracking  error  rate 
weight  in  Qc  is  chosen  to  be  one-fourth  of  the  tracking  error 
weight  as  suggested  in  Ref.  36.  Based  on  Fig.  45 f ,  lateral 
velocity  is  weighted  and  because  the  pilot  seemed  to  place 
more  emphasis  on  reducing  A<j>v  variations  rather  than  A6<t>  vari¬ 
ations  in  Fig.  Aba,  Abv  is  weighted  in  Qc .  At  the  low  aQ  and 
medium  ao  flight  conditions  in  Table  1A,  a  high  weight  (2)  and 
a  low  weight  (1)  Qc  are  specified.  The  rest  of  the  pilot  model 
parameters  are  standard  (t=0.2  sec,  tn  =  0.1  sec,  Py  *  -  20 
db,  Pu  =  -30  db) . 

The  resulting  pilot  model  closed-loop  eigenvalues  and 
feedback  gain  elements  are  shown  in  Tables  16  and  17  respec¬ 
tively.  The  eigenvalues  and  gains  exhibit  little  change  when 
Qc  changes  between  cost  functions  1  and  2  in  Table  15.  The 
closed-loop  eigenvalues  are  not  very  satisfactory  from  a  con¬ 
trol  point  of  view.  At  =  9.3  deg  the  sas/tracking  error 
mode  combination  has  poor  damping  whil.e  the  roll  angle  eigen¬ 
value  at  aQ  =  5.2  deg  is  almost  neutrally  stable.  The  pilot 
model  control  gains  at  ao  =  5.2  deg  and  aQ  =  9.3  deg  are  simi¬ 
lar,  while  the  high  pilot  model  gains  have  significant  dif¬ 
ferences.  These  differences  indicate  the  hypothesis  testing 
scheme  should  be  able  to  make  a  clear  decision  between  the 
high  and  low  aQ  strategies. 


1.2A 


•-j|r.iir-v,i',duaeaq 


TABLE  15 

SQUARE  ROOT  OF  PILOT  MODEL  COST  FUNCTION  WEIGHTS 
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AIRCRAFT 
ANGLE  OF 
ATTACK 

COST 

FUNCTION 

LATERAL 

VELOCITY 

ROLL 

ANGLE 

TRACKING 

ERROR 

TRACKING 
ERROR  RATE 

5.2 

1 

0.25 

0.50 

5.0 

1.25 

5.2 

2 

0.50 

1.00 

10.0 

2.5 

9.3 

1 

0.25 

0.25 

5.0 

1.25 

9.3 

2 

0.25 

0.25 

10.0 

2.5 

19.5 

1 

0.25 

0.25 

5.0 

1.25 

A. 8  HYPOTHESIS  TESTING  RESULTS 

This  section  presents  the  results  from  the  hypothesis 
testing  scheme  in  two  formats.  In  the  first  format,  the  one- 
step  predicted  pilot  model  control  time  history  is  overlaid  on 
the  actual  pilot's  lateral  stick  control  time  history.  The 
closer  these  two  trajectories  match,  the  better  the  indication 
is  that  the  assumed  pilot  model  is  correct.  The  second  format 
is  composed  of  plots  of  the  probabilities  that  a  pilot  model 
is  the  best  of  those  tested.  The  first  format  gives  an  abso¬ 
lute  indication  of  pilot  model  performance  while  the  second 
format  gives  a  relative  indication  of  performance  between  the 
pilot  models. 

Initial  trials  of  the  hypothesis  testing  calculation 
procedure  with  the  piloied  simulation  data  indicated  a  severe 
mismatch  occurs  between  the  Kalman  filter's  estimate  of  its 
output  error  covariance,  S  in  Eq .  26,  and  the  actual  output 
error  covariance,  which  is  calculated  as 
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For  a  large  number  of  data  points  (k  is  large),  the  matrix 

should  approach  S  for  a  linear,  time-invariant  system.  The 

mismatch  occurs  because  the  pilot's  internal  model  of  the  tar- 

2 

get  driving  covariance,  ,  appears  to  differ  from  the  actual 

value  calculated  in  Section  4.5,  and  also  due  to  the  existence 
of  unmodeled  nonlinearities  and  the  time-varying  nature  of  the 
actual  piloted  simulation.  The  pilot  model  probabilities, 
which  strongly  depend  on  S,  cannot  be  calculated  in  the  pres¬ 
ence  of  this  mismatch. 


One  approach  to  resolve  this  mismatch,  suggested  in 
Ref.  35,  is  to  adjust  the  pilot's  internal  model  of  the  target 
driving  covariance  until  the  mismatch  is  minimi;, ’.ed .  Since 
this  requires  a  large  number  of  iterations  of  the  pilot  model 
Vcval  culations  and  the  hypothesis  testing  procedure,  this  ap¬ 
proach  is  not  pursued  here.  Rather,  a  simpler  approach  is 
used  whereby  an  estimate  of  S,  based  on  Eq .  29  applied  to  the 
entire  1280  samples  of  data  is  used  in  the  hypothesis  testing 
scheme.  This  procedure  yields  well-behaved  pilot  model  proba¬ 
bilities  from  the  hypothesis  testing  scheme.  The  usefulness 
of  this  modification  is  confirmed  by  the  agreement  observed  in 
the  next  paragraphs  between  the  hypothesis  testing  results  and 
the  pilot  control  trajectory  comparisons. 

Comparisons  of  the  actual  pilot  control  time  history 
to  the  pilot  model  predictions  are  illustrated  in  Fig,  47.  The 
low  a  pilot  model  with  normal  piloting  technique  shown  in  Fig. 
47a  matches  the  actual  pilot's  responses  quite  well,  except 
perhaps  in  the  regions  circled.  The  ao  pilot  model  with  es¬ 
pecially  "tight"  piloting  technique  (Fig.  47b)  exhibits  no 
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a)  Low  a  Pilot  Model  With  Nominal  Quadratic  Cost  Weights 
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b)  Low  a  Pilot  Model  With  High  Quadratic  Cost  Weights 


Figure  47 


Comparison  of  Actual  and  Predicted  Pilot  Controls 
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c)  Medium  a  Pilot  Model  With  Nominal  Quadratic  Weights 


TIME  (i.c)  TIME  (mc) 

d)  High  a  Pilot  Model 

Figure  47  Comparison  of  Actual  and  Predicted  Pilot 
Controls  (Continued) 


29 


I 


discernable  difference  from  the  low  aQ  nominal  pilot  model. 

The  medium  a  pilot  model  {Fig.  47c)  does  not  match  the  pilot 
control  inputs  as  well  as  the  low  a  pilot  model  except  perhaps 
within  the  circled  regions.  The  high  aQ  pilot  model  does  a 
poor  job  of  predicting  pilot  behavior  even  in  high  oQ  flight. 
This  pilot  does  not  signlfii  v.-tly  adapt  his  control  strategy 
to  match  the  changes  in  aircraft  dynamics.  Note  that  no  con¬ 
clusion  regarding  quadratic  weight  ("tight"  or  "normal"  con¬ 
trol  ef fort-tracking  error  trade-off)  can  be  made  based  on 
visual  inspection  of  these  figures. 

The  hypothesis  testing  results  are  shown  in  Fig.  48. 
Recall  that  these  results  take  the  form  of  a  set  of  probabili¬ 
ties  for  each  of  the  hypothesized  pilot  models,  with  the  sum 
of  the  probabilities  of  all  hypotheses  in  a  given  test  equal 
to  one.  Allowing  the  hypothesis  testing  algorithm  to  choose 

between  the  low  a  and  medium  of  pilot  models  results  in  the 

o  o 

probability  histories  shown  in  Fig.  48a.  The  low  ofQ  model 
best  matches  the  pilot's  responses  except  in  the  circled 
regions.  The  hypothesis  testing  result  confirms  the  qualita¬ 
tive  analysis  of  the  control  time  histories.  The  addition  of 
the  high  aQ  pilot  model  (48b)  to  the  set  of  hypotheses  pro¬ 
duces  n  change;  at  no  time  does  the  high  o<0  model  predict  the 
actual  pilot's  behavior,  even  when  the  aircraft  really  is  in  a 
high  ft  flight  condition, 

Choice  of  pilot  tracking  error-control  effort  trade¬ 
off  (quadratic  weight  level)  could  not  be  made  from  a  visual 
examination  of  the  control  responses,  but  the  hypothesis  test¬ 
ing  algorithm  clearly  indicates  that  the  "normal"  weights  give 
a  much  better  fit  to  the  data  than  the  "high"  weights.  This 
indicates  that  the  "high-gain"  or  "tight"  pilot  model  is  not 
as  good  a  representation  of  what  the  pilot  actually  did  as  is 
the  normal  or  "low-gain"  pilot  model  values. 
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c)  Pilot  Model  at  a  =  5  deg  with 
Different  Quadratic  Weights 


Figure  48  Pilot  Model  Probabilities  (Continued) 

These  results  apply  to  the  pilot  trajectory  analyzed 
here,  It  is  apparent  that  analysis  of  many  more  pilot  trajec¬ 
tories  is  necessary  before  general  conclusions  can  be  made. 
Recall  from  Chapter  3  that  the  use  of  pedals  at  high  aQ  to 
roll  the  F-14A  is  desirable.  The  pilot  in  this  test  was  spec¬ 
ifically  requested  not  to  use  rudder  pedals,  and  so  it  is  not 
clear  what  this  pilot's  high-ao  responses  would  have  been  if 
pedals  had  been  allowed, 


4 . 9  CHAPTER  SUMMARY 

Applications  of  the  optimal  control  pilot  model  in  Refs. 
3  and  23  show  great  promise  for  expanding  and  understanding  sta¬ 
bility  and  performance  characteristics  of  piloted  aircraft  in 
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maneuvering  flight.  The  results  in  this  chapter  have  substan¬ 
tiated  these  results  by  postprocessing  simulated  flight  test 
data  to  determine  actual  control  strategies  used  by  a  human 
pilot.  The  primary  result  is  that  the  pilot,  throughout  most 
of  the  simulation  studies,  did  not  change  his  internal  model 
of  the  aircraft's  dynamics  from  a  low  angle  of  attack  flight 
condition  to  others  that  the  aircraft  traversed ;  i . e . ,  the 
pilot  does  not  appear  to  adapt  his  control  strategy  to  match 
changes  in  the  aircraft's  behavior.  It  is  suggested  that  the 
reason  the  pilot  did  not  adapt  is  the  same  as  the  reason  given 
in  Ref.  3  for  the  minimum  control  effort  strategy.  To  adapt, 
the  pilot  would  have  to  increase  his  control  effort  by  subjec¬ 
tively  increasing  his  feedback  gains  as  shown  in  Table  17; 
increasing  control  effort  is  something  pilots  apparently  try 
to  avoid.  The  optimal  control  pilot  model  predicts  control 
signals  which  closely  match  those  of  the  actual  pilot  when  the 
optimal  control  pilot  model  has  a  fixed  low  aQ  control  strategy 
The  implications  of  these  results  are  clear, 


•  If  the  pilot  model's  low  angle  of  attack 
control  strategy  destabilizes  the  air¬ 
craft  at  high  angles  of  attack  then  a 
human  pilot  will  also  have  stability 
problems,  since  it  cannot  be  assumed 
that  he  will  adapt  to  the  actual  air¬ 
craft  flight  condition. 

•  An  automatic  control  law  (sas,  cas,  AKI , 
DFCS  etc.)  should  be  tested  with  a  non¬ 
adapting  pilot  model  to  assess  its  capa¬ 
bilities.  A  very  simple  nonadapting 
pilot  model  is  used  in  Chapter  3  to 
successfully  test  the  DFCS  designed  in 
Chapter  2,  for  example. 


A  number  of  other  results  have  been  obtained  in  this  study 


•  A  modified  hypothesis  testing  scheme 
combined  with  a  linear  pilot  model  is 
developed  that  can  be  used  to  analyze 
nonlinear  piloted  aircraft  data. 

•  The  hypothesis  testing  scheme  and  the 
simple  procedure  used  to  handle  total 
value  aircraft  states  and  controls 
(i.e.,  high-pass  filtering)  are  suc¬ 
cessfully  validated  against  computer 
generated  synthetic  data  for  which  the 
true  hypothesis  is  known 

•  The  pursuit-tracking  target/aircraft  mod- 
el  developed  in  this  report  can  be  usecl 
Tn  any  future  pilot-aircraft  studies 
performed  with  (as  in  this  study)  or 
without  (Refs.  3  and  23)  pilot  data  to 
predict  pilot  control  behavior.  Spec¬ 
ifically,  lateral-directional  and  longi¬ 
tudinal  tracking  can  be  analyzed  together 
to  insure  that  all  appropriate  couplings 
are  examined. 


The  hypothesis  testing  scheme  with  the  optimal  control  pilot 
model  has  considerable  room  for  improvement  and  continued  a- 
nalysis.  Recommended  extensions  are  outlined  in  Chapter  5. 


5. 


CONC LUSIONS  AND  RECOMMENDATIONS 


5.1  SUMMARY  AND  CONCLUSIONS 

The  main  results  of  this  report  are  that  a  digital 
flight  control  system  for  air  combat  can  be  constructed  which 
demonstrates  better  handling  qualities  than  a  conventional 
control  system  along  typical  acm  trajectories,  and  that  it  is 
possible  to  choose  (from  a  set  of  pilot  models)  the  pilot  model 
that  best  matches  actual  pilot  control  motions. 

Chapter  2  of  this  report  is  devoted  to  the  construc¬ 
tion  of  an  Air  Combat  Digital  Flight  Control  System  (DFCS), 

This  was  demonstrated  for  the  Grumman  F-14A  as  an  example. 

This  DFCS  exhibits  the  two  major  properties  necessary  for  a 
command  augmentation  system: 

•  Modification  of  vehicle  stability  charac¬ 
teristics  to  provide  frequency  and  damping 
that  match  MIL-F-8785B  Level  1,  Category 

A  requirements  at  all  flight  conditions 
in  the  design  regime. 

•  Modification  of  vehicle  control  response 
characteristics  to  provide  fast  and  ac¬ 
curate  response  to  pilot  commands.  Addi¬ 
tionally,  this  response  is  uniform  over 
the  design  regime. 

It  should  b*  emphasised  that  the  design  regime  spans 
the  entire  range  of:  velocity,  normal  acceleration  and  roll 
rate  expected  in  air  combat  maneuvering  (ACM)  flight.  Sample 
rates  and  control  calculation  requirements  are  similar  to  thoae 
of  a  flight  control  algorithm  designed  by  TASC  that  has  been 
test  flown  as  part  of  the  NASA  VALT  program  (Ref.  7). 
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Chapter  3  details  the  testing  of  the  Air  Combat  DFCS 
against  the  conventional  control  system,  A  summary  of  the 
results  and  conclusions  from  these  tests  is  as  follows: 


•  The  DFCS  provides  precise  stability  aug¬ 
mentation  in  steady  turns  and  rapid  roll¬ 
ing'.  The  Dutch  roll  mode  is  well  damped, 
and  longitudinal  maneuvers  are  quickened 
by  the  automatic  maneuver  flap  feature 
built  into  the  DFCS. 

•  The  DFCS  provides  precise  command  response, 
with  more  accurate  normal  acceleration 

and  much  more  accurate  roll  rate  response. 
Automatic  turn  coordination  operates  by 
nulling  sideslip  excusi.ons  even  in  rapid 
maneuvering  flight. 

•  The  DFCS  handles  control  saturation  well, 
providing  good  stability  and  accurate 
command  response  right,  up  to  the  perfor¬ 
mance  bounds  of  the"airframe I  Control 
saturation  is  used  by  tfiebFcs  (as  is 
common  in  ACM  flight)  to  produce  rapid 
and  accurate  command  response. 

•  Conventional  piloting  technique  of  the 
aircraft  with  a  conventional  mechanical 
control  system  produced  a  control-induced 
departure  and  incipient  spin  in  one  ACM 
test . 

•  Loss  of  lateral  stick  roll  control  power 
for  the  conventional  control  system  in 
high-a  requires  the  pilot  to  compensate 
by  using  rudder  pedals.  The  DFCS  does 
this  automatically,  enabling  the  same 
piloting  technique  to  be  used  in  all 
flight  conditions, 

*'  The  DFCS  can  be  implemented  with  conven¬ 
tional  sensor  inputs,  Except  for  pilot 
control  inputs,  the  present  sensor  suite 
ohould  be  sufficient . 

•  Tests  of  digital  control  of  the  present 
limited-authority  sas  actuators  with  the 
conventional  pilot's  mechanical  control 
system  gives  preliminary  indication  that 
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most  of  the  DFCS  advantages  can  be  gained 
by  digital  control  of  the  present  limited- 
authority  sas  actuators, 


Chapter  A  of  this  report  presents  a  method  for  eval¬ 
uating  actual  pilot  control  actions  (with  data  from  the  NASA 
Differential  Maneuvering  Simulator)  with  hypothesized  pilot 
model  predictions.  The  different  pilot  strategies  tested  here 
are  adapted  to  low-a ,  medium-a ,  and  high-a  flight  conditions. 
Also,  a  pilot  model  based  on  tight  tracking  error-high  control 
effort  is  compared  to  a  pilot  model  based  on  loose  tracking 
error  -  lower  control  effort.  The  processing  of  piloted- 
simulation  trajectory  data  led  to  the  following  results  and 
conclusions : 


•  The  hypothesis  testing  method  can  choose 
the  candidate  model  even  when  differences 
are  not  apparent  during  visual  inspection 
of  the  trajectories. 

•  The  low  angle  of  attack  optimal  control 
pilot  model  predicted. the  pilot  input 
well,  whereas  the  medium  angle  of  attack 
model  only  appeared  to  be  better  in  cer¬ 
tain  very  short  trajectory  segments. 

The  pilot’s  control  strategy  did  not 
correspond  to  the  high  angle-of-attack 
model  at  all,  even  during  high  angle- 
of-attack  flight. 

•  During  the  piloted-simulation  evaluated 
here,  the  pilot  did  not  adapt  his  control 
strategy  to  the  angle-of-attack  at  which 
the  aircraft  was  flying.  Hence a  pilot/ 
aircraft  stabiTity  study  such  as  that 
performed  by  TASC  in  Ref.  3  should  pro¬ 
vide  useful  information. 


5 . 2  RECOMMENDATIONS 


In  the  area  of  control  system  design,  investigation 
of  accelerometer  feedbacks  rather  than  aerodynamic  angle  feed¬ 
back  is  desirable.  This  can  be  achieved  by  a  digital  estimator 
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which  uses  accelerometer  output  to  estimate  the  aerodynamic 
angle  feedbacks  necessary  in  the  present  control  algorithms, 
or  the  control  algorithms  can  be  designed  for  accelerometer 
feedbacks  instead  of  aerodynamic  angle  feedbacks.  Theoretical 
solutions  to  both  approaches  are  available,  with  the  second 
being  desirable  from  a  control  system  simplicity  viewpoint. 

In  any  event,  aerodynamic  angles  will  still  be  necessary  for 
gain  scheduling.  The  choice  between  the  two  depends  on  sensor 
suite  and  sensor  accuracy  and  noise  properties.  The  output  of 
this  study  would  consist  of  control  system  designs  better  adapted 
to  the  present-day  fighter  sensor  suites, 

In  the  area  of  control  system  gain  scheduling,  a  working 
ACM  control  system  must  be  designed  for  additional  types  of 
flight  conditions,  The  three  parameters  used  in  this  study 
remain  the  most  important  (speed,  normal  acceleration,  and 
roll  rate),  although  the  ranges  may  require  some  extension. 

The  addition  of  flight  conditions  with  different  altitude, 

Mach  number,  eg  location,  and  weight  allows  a  determination  of 
the  importance  of  these  effects  and  enables  an  examination  of 
Air  Combat  DFCS  performance  changes  due  to  these  flight  condi¬ 
tion  variations.  The  output  of  this  task  would  be  an  evalua¬ 
tion  of  the  importance  of  the  various  gain-scheduling  parameters. 

In  the  area  of  digital  augmentation  of  a  conventional 
aircraft,  design  of  a  version  of  the  Air  Combat  sas  optimized 
for  operation  through  the  limited-authority  sas  actuators  should 
be  aggressively  pursued.  This  approach  may  provide  a  signifi¬ 
cant  improvement  in  air  combat  handling  qualities  without  major 
modification  to  the  conventional  mechanical  control  system. 
Retrofitting  a  digital  sas  computer  during  the  F-14A  CILOF 
program  can  possibly  be  justified  by  these  improvements . 

Previous  aircraft  stability  and  control  studies  and 
pilot-aircraft  stability  studies,  along  with  the  air  combat 
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evaluation  technique  presented  here  form  the  basis  of  an  analy¬ 
sis  capability  upon  which  an  improved  analog  sas  with  aileron- 
rudder  interconnect  for  good  high  angle-of-attack  handling  qual 
ities  can  probably  be  designed.  Although  the  handling  qualities 
improvements  are  not  likely  to  be  as  great  for  this  approach 
as  for  the  digital  sas  approach,  the  analysis  tools  do  make 
this  approach  feasible. 

Even  though  the  pilot  model  testing  performed  in  this 
study  indicated  that  the  pilot  flying  the  aircraft  simulation 
analyzed  here  did  not  adapt  to  the  actual  flight  condition  of 
the  aircraft,  it  cannot  be  discerned  whether  this  is  typical 
of  ACM  pilots  in  general.  Analysis  of  pilot  data  over  many 
trajectories  would  be  necessary.  An  additional  data  set  gen¬ 
erated  by  NASA  Langley  is  available  and  can  be  examined  using 
the  pilot  model  hypothesis  testing  programs  developed.  There 
are  a  series  of  improvements  to  the  approach  used  here  which 
should  be  pursued  in  the  event  of  additional  pilot  data  analy¬ 
sis,  the  most  important  of  which  is  improving  the  method  used 
to  select  the  pilot's  internal  model  of  the  target  driving 
covariance.  This  would  enable  the  use  of  the  Kalman  filter's 
estimate  of  it's  output  error  covariance  in  the  hypothesis 
testing  procedure.  Other  questions  that  can  be  answered  in 
further  data  analyses  include: 

•  Do  other  pilot  models  predict  pilot 
response  as  accuratey  as  the  optimal 
control  model? 

•  Does  the  pilot's  control  effort-tracking 
error  trade-off  vary  along  an  ACM  trajec¬ 
tory? 

•  Do  different  pilots  exhibit  different 
adaptation  strategies  or  control  effort¬ 
tracking  error  tradeoffs? 

•  Do  additions  to  the  pilot's  observation 
set  improve  the  pilot  model's  match  of 
pilot  response? 
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•  In  the  multi-axis  case,  can  attention 
allocation  between  axes  be  identified? 

Some  of  these  questions  can  be  answered  by  further  analysis  of 
available  data,  while  some  of  them  require  additional  data, 

More  realistic  target  and  target  dynamics  models  can 
be  incorporated  in  the  piloted  simulation  analysis.  While  the 
analysis  proceeded  reasonably  well,  additional  target  data 
would  reduce  the  need  for  simplifications  to  the  target  and 
target  dynamics  models  made  in  Appendix  B.  It  is  not  known  if 
these  simplifications  reduced  the  accuracy  of  the  pilot  model 
analysis  performed  here;  an  investigation  of  more  accurate 
models  would  answer  this  question. 
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APPENDIX  A 

AIRCRAFT  EQUATIONS  OF  MOTION 

mu  11 "  -  -*»  1 " v  1  ’  . .  IIJ'  1 


A . 1  OVERVIEW 

This  appendix  reviews  the  equations  of  motion  for  a 
maneuvering  aircraft.  These  derivations  make  extensive  use  of 
vector-matrix  differential  equation  ("state-space")  notation, 
and  the  resulting  equations  are  written  in  a  form  suitable  for 
analysis  using  concepts  of  modern  control  theory. 


A. 2  NONLINEAR  DYNAMIC  EQUATIONS 

The  derivation  of  nonlinear  rigid-body  equations  is 
reviewed  in  this  section.  The  equations  are  developed  using 
"flat-earth"  assumptions,  i.e.,  the  effects  of  earth  curvature 
and  rotation  are  assumed  negligible.  This  means  that  earth- 
fixed  and  inertial  reference  frames  are  equivalent. 

For  the  moderate  velocities  of  interest  in  this  report 
(typically  below  305  m/s  (1000  fps)),  the  equivalence  of  earth- 
fixed  and  inertial  reference  frames  is  a  good  assumption.  The 
origin  of  the  inertial  reference  frame  used  here  is  located  on 
the  surface  of  the  earth,  with  the  x- ,  y- ,  and  z-axes  in  a 
north-east-down  orientation.  Since  the  simplest  statement  of 
Newton's  Second  Law  is  given  in  an  inertial  reference  frame, 
this  frame  plays  an  important  part  in  the  derivation  of  the 
dynamic  equations. 


For  point-mass  analysis,  the  velocity  reference  frame 
is  convenient.  The  origin  lies  at  the  eg  of  the  vehicle,  and 
the  x-axis  points  along  the  velocity  vector.  The  y-axis  is 
horizontal  and  the  z-axis  lies  in  the  vertical  plane  which 
includes  the  x-axis.  The  velocity  heading  angle,  £,  describes 
the  heading  of  the  projection  of  the  velocity  x-axis  on  the 
horzontal  plane  relative  to  north,  and  the  flight  path  angle, 

Y,  describes  the  inclination  of  the  velocity  vector  relative 
to  the  horizontal. 

The  wind  axes  also  include  an  x-axis  along  the  velo¬ 
city  vector,  but  the  z-axis  lies  in  the  plane  of  symmetry  of 
the  vehicle  veing  described.  Hence,  as  t.he  vehicle  rolls  about 
the  velocity  vector,  the  wind  axes  roll  also,  The  roll  angle 
relative  to  the  velocity  axes  is  <J> v . 

Detailed  vehicle  state  equations  can  best  be  expressed 
in  body- fixed  rather  than  velocity-fixed  axes.  These  are  the 
axes  in  which  the  pilot,  the  sensors,  and  the  control  surface 
locations  are  defined.  Body  axes  are  the  only  axes  in  which 
the  moment-of-inertia  matrix  is  constant,  and  also  dynamic 
data  collected  from  sting-mounted  wind  tunnel  models  or  from 
flight  tests  usually  are  expressed  in  body  axes. 

The  various  body- fixed  axis  systems  have  a  common 
origin,  located  at  the  body  center  of  mass,  and  are  fixed  in 
orientation  with  respect  to  the  vehicle.  Generally,  the  body 
x-axis  extends  forward  out  the  vehicle's  nose,  the  y-axis  ex¬ 
tends  out  the  right  wing,  and  the  z-axis  extends  out  the  bottom 
of  the  vehicle.  The  x-z  plane  is  usually  a  plane  of  geometric 
symmetry,  if  the  vehicle  has  one.  There  are  a  number  of  possi¬ 
ble  body-fixed  reference  frames,  and  the  one  fixed  by  the  builde 
is  simply  referred  to  in  this  report  as  the  body-axis  system. 

For  any  nominal  flight  condition,  body- fixed  axes  can  be  chosen 


so  that  the  x-axis  is  aligned  with  the  velocity  vector,  and 
the  z-axis  is  in  the  body-axis  plane.  This  set  of  body-fixed 
axes  is  referred  to  as  the  stability-axis  system. 

A  transformation  from  inertial  to  body  axes  is  com¬ 
posed  of  a  right-handed  yaw  through  an  angle  i|j,  then  a  right- 
handed  pitch  through  an  angle  6,  and  then  a  right-handed  roll 
through  an  angle  <j> .  The  body  axes  and  wind  axes  are  related  to 
each  other  through  the  aerodynamic  angles,  angle  of  attack,  or , 
and  angle  of  sideslip,  p.  Angle  of  attack  represents  the  air¬ 
craft  body  pitch  angle  above  the  velocity  vector,  and  sideslip 
is  the  angle  that  the  aircraft  nose  is  yawed  left  of  the  velo¬ 
city  vector. 

In  summary,  the  inertial  axes,  body  axes,  wind  axes 
and  velocity  axes  can  be  related  to  each  other  by  the  trans¬ 
formations  Illustrated  in  Fig.  A.l.  The  inertial -to-body -axis 
transformation  procedes  as  follows: 

H  j  ( 4> ,  6  ,<|0  =  h|<*)  hJ<0) 


10  0 

-  — 

cos0  0  -sinB 

cosi|j  sim|i  0 

s 

0  cos((i  sinij) 

0  10 

-  s  imp  cost})  0 

0  -  s in4>  cosiji 

sin0  0  cosB 

0  0  1 

(  A- 1 ) 

For  orthonomal  matrices  such  as  these,  the  matrix  inverse, 

(  )  ,  is  equal  to  the  transpose,  (  )  .  Equation  A-2  il¬ 

lustrates  the  inertial-to-velocity  axis  transformation  matrix. 

hx(y.O  =  h^h(y)  HjH(£> 
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Reference  Axis  Transformations  (Arrows 
Indicate  Right-Hand  Rotation) 
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The  velocity- to-wind  axis  transformation  is  defined  by 


c  o  &  (Ji  y  s 


and  the  wind-body  transformation  matrix  is  given  by 


H*<a,-{3)  =  H^<«)  hJ(-P) 
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(A-4) 


In  the  remainder  of  this  section,  the  vehicle's  equa¬ 
tion  of  motion  is  derived  as  a  single  state-vector  equation  of 
the  form 


x  =  f(x,u)  (A-5) 

where  x  is  the  state  vector,  u  is  the  control  vector,  f  is  the 
vector  system  dynamics  equation,  and  disturbances  are  neglect¬ 
ed,  The  state  vector  is  a  12-element  vector,  and  the  nonlinear 
state  equations  are  readily  derived  as  four  sets  of  three  equa¬ 
tions  representing 

•  Translational  Kinematics 

9  Rotational  Kinematics 

•  Translational  Dynamics 

•  Rotational  Dynamics, 

The  kinematic  equations  relate  the  vehicle's  translational  and 
rotational  velocities  to  its  position  in  inertial  space,  and 
thus  involve  body-axis/inertial-axis  relationships.  The  dy¬ 
namic  equations  describe  the  changes  of  the  vehicle  velocity 
components  caused  by  the  applied  forces  and  moments;  they  are 
best  derived  in  a  body-fixed  frame  of  reference. 
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A  second  nonlinear  vector  equation  is  also  derived  to 
relate  the  output  of  the  aircraft  to  its  state  and  control. 

The  output  equation, 


y  =  h(x,u)  (A-6) 

is  necessary  in  any  control  design  problem  because  the  goal  of 
control  system  design  is  to  cause  the  output  to  follow  a  de¬ 
sired  trajectory.  The  specific  components  of  the  output  vector 
depend  on  the  vehicle  and  the  tactical  situation. 

A. 2.1  Kinematics 


The  translational  and  angular  position  of  the  vehicle 
are  given  relative  to  inertial  space  by  the  inertial  position 
vector,  x-j.  and  by  the  inertial-body  Euler  angle  vector,  vB: 


(A-7) 


(A-8) 


It  is  important  to  note  that  the  Euler  angle  "vector"  is  not  a 
true  vector  in  physical  space;  it  is  an  ordered  triple  of  right- 
handed  rotations  which  occur  about  different  axes  of  different 
reference  frames. 


The  translational  and  angular  rate  vectors  often  are 
expressed  in  body  axes,  as  in  the  following: 
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(A-9) 


(A-10 ) 


The  body-axis  translational  rate  vector,  vfi,  is  ari  expression, 
in  body  axes,  of  the  derivative  of  the  inertial  position  vector. 
This  relationship  supplies  the  first  part  of  the  nonlinear 
state  equations  of  motion 

Xj  =  Hg  Vg  (A- 1.1) 

where  Hg  is  the  inverse  of  the  inertial-body  transformation 
derived  in  Eq.  (A-l ) . 


The  body  angular  rate  vector  also  can  be  related  to 
the  derivative  of  the  Euler  angle  vector  by  noting  that  the 
Euler  angle  derivatives  occur  in  three  different  reference 
frames.  The  resulting  transformation  is  constructed  in  Eq, 
(A-12),  where  the  individual  transformations  are  the  same  as 
those  of  Eq.  (A-l) : 
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o  ^ 

-  — B  (A-12) 

In  this  equation, 
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The  ordering  of  the  transformations  in  Eq.  (A-12)  arises  from 
the  ordering  of  the  Euler  angles.  As  can  be  seen  from  Fig. 
A-l,  the  angular  rate,  q,  occurs  about  the  Xg  axis;  the  rate 
0,  occurs  about  the  axis;  and  tji  occurs  about  the  z^  axis. 
The  inverse  of  Eq,  (A-ll)  supplies  the  rotational  kinematic 
part  of  the  vehicle  nonlinear  state  equations,  and  it  is  given 
by: 


vB  =  L"1  wj  (A-13) 


A. 2 . 2  Dynamics 


The  dynamics  of  the  vehicle  involve  Newton's  Second 
Law,  which  equates  the  applied  force  to  the  time  derivative  of 
inertial  translational  momentum  of  a  body,  For  rotational 
motion,  this  equivalence  becomes  one  between  torque  and  the 
derivative  of  angular  momentum,  measured  in  an  inertial  ref¬ 
erence  frame. 


An  expression  for  the  inertial  translational  accel¬ 
eration,  expressed  in  body-axis  variables,  can  be  derived  from 
Eq.  (A-10)  by  taking  the  derivative  of  both  sides: 

-1  =  HB^B  +  ^B-B  (A- 14) 

Here  , 

Hg  =  H*  uig  (A-15) 

and  uig  is  the  cross-product  equivalent  matrix  for  Wg  given  by 
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This  leads  to  the  body  axis  equation 

VB  =  H®*!  -  OiBvB 


(A- 16) 


(A-17) 


The  applied  specific  forces  consist  of  gravitational 
forces  and  contact  forces.  In  inertial  axes,  the  gravity  force 
is 

0  “ 

0  (A-18) 

e 

The  specific  contact  force  can  be  broken  into  two  components, 
one  of  which  is  due  to  aerodynamic  forces,  Fg ,  and  one  of  which 
is  due  to  thrust,  Tg : 

'x/m 

Fg/m  =  Y/m 
Z/m 

b/” 

B/" 

B  /m 
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(Capital  letters  are  conventionally  used  in  aerodynamics  to 
denote  the  force  components.)  The  specific  equations  for  the 
aerodynamic  and  thrust  foioes  are  given  in  Ref.  23. 

The  translational  dynamic  equation  is  formed  by  equa 
ting  the  sum  of  the  aerodynamic  and  gravitational  specific 
forces  to  the  inertial  translational  acceleration  of  the  ve¬ 
hicle.  Then  Eq.  (A-17)  becomes 

VB  *  <E„  +  TB)/m  H*g,  -  SBvB  (A-21) 


To  construct  the  rotational  dynamic  equation,  an  ex¬ 
pression  for  the  time  derivative  of  angular  momentum  measured 
in  inertial  axes  is  necessary.  The  angular  momentum,  hg ,  is 
most  easily  expressed  in  body  axes;  neglecting  rotating  mach¬ 
inery,  it  is  the  product  of  the  moment-of-inertia  matrix  (con¬ 
stant  in  body  axes)  and  the  angular  rate  vector 

hB  =  h  ~B  (A-22 ) 

where  the  inertia  matrix  contains  all  products  and  moments  of 
inertia : 
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(A-23 ) 


The  time  derivative  of  the  angular  momentum,  expressed  in  iner¬ 
tial  axes,  is  derived  as 


HBIBujB 


1 B-B 


BBB  B  B  B-B 


( A-24 ) 
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The 

components . 


contact  moments  consist  of  aerodynamic  and  thrust 
These  are  defined  as 


(A-25 ) 


(A-26) 


(Capital  letters  are  conventionally  used  for  the  moment  compo¬ 
nents.)  The  specific  equations  for  the  aerodynamic  and  thrust 
moments  are  given  in  Ref.  23.  The  rotational  dynamic  equation 
is  formed  by  equating  the  applied  torques  to  the  derivative  of 
the  angular  momentum: 

"s  =  hVb  +  G„)  -  ijVkSb  (a'27> 


A . 2 . 3  S ummary  of  State  Equations 


The  nonlinear  state  equations  are 


*1  =  Yg 

D 


-B  =  LB  -B 


Yb  “  (-B  +  -B^/m  +  HI^I  "  WB-B 


(A-28) 

(A-29 ) 

(A-30 ) 
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(A-31) 


These  equations  fall  into  the  general  state  equation  form 


X  =  f(x,u) 


by  defining  the  state  vector  as 


T  T  T  T 
*1  vb  vb  wB 


(A-32) 


(A-33) 


and  noting  that  the  aerodynamic  forces  and  moments  are  functions 
of  the  states,  controls,  disturbances  and,  to  some  extent,  the 
state  history. 

A . 2 . 4  Output  Equations 

The  output  vector,  consists  in  general  of  nonlinear 
functions  of  the  aircraft  states,  x,  and  controls,  u,  and  is 
given  by 


1  ~  h(x.u)  (A-34) 

For  a  conventional  fighter  aircraft,  the  four  basic 
commanded  motions  are  longitudinal,  lateral,  normal,  and  di¬ 
rectional  motions.  Longitudinal  motion  results  in  a  velocity 
magnitude  change  and  can  be  commanded  by  V  or  V.  Lateral  (roll¬ 
ing)  motion  is  used  to  orient  the  maneuver  plane  and  can  be 
commanded  by  p,  p  ,  or  <t> .  Normal  and  directional  plane  mo¬ 
tions  are  two  degree-of-freedom  motions,  and,  in  general,  re¬ 
quire  two  commands.  In  the  normal  plane,  acceleration  (an  or 
q)  and/or  attitudes  (0,  a,  or  y)  can  be  commanded,  with  the 
two-element  directional  command  vector  chosen  in  an  analogous 
wav.  All  of  these  commands  are  desirable  in  one  situation  or 
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another.  In  ground  attack,  both  flight  path  control  (y)  and 
independent  fuselage  pointing  (a)  might  be  desirable.  In  air 
combat  maneuvering,  normal  acceleration  (an)  is  certainly  a 
useful  command,  as  is  stability-axis  rollrate  (Pw)*  Output 
equations  for  each  of  these  quantities  are  needed  so  that  they 
can  be  used  in  control  system  design. 

The  aircraft  velocity  and  the  aerodynamic  angles  are 
given  in  terms  of  the  body  axis  velocities  as  follows 


V 

Ju2  +  v^  +  w^ 

p 

* 

tan"1  (v/Ju2  +  w^) 

a 

tan-1  (w/u) 

(A-35 ) 


The  three  components  of  earth-relative  acceleration  expressed 
in  wind  axes  are  given  in  terms  of  body-axis  quantities  as 
follows 


fv 


where  it  should  be  noted  that  these  quantities  are  different 
than  those  usually  named  a^  and  an ,  as  the  acceleration  of 
gravity  is  included  in  these  terms.  Hence,  a  command  system 
set  up  on  this  basis  automatically  compensates  for  gravity. 

The  earth-relative  acceleration  in  wind  axes  can  also 
be  specified  in  terms  of  the  velocity-axis  angular  rates  as 
follows 


•W, 


=  Hj  «x  =  H 


\}b  +  “B 


(A-36) 
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HW  JtT  =  HW  HV  X 
1  V  I  "I 


1  0 


0  v 


=  0  cos«j)  sin<}) v  0  Vcosy  0  i  (A-37) 


0  -  sin<|>  cosijiv  0  0  -V  y 


The  second  two  equations  result  in  the  following  equations  de¬ 
scribing  the  velocity  angle  dynamics 


V  cosy  1  =  cos$v  +  an  sin<t>^ 

V  y  =  -  ay  sin<f>v  +  an  cos0v 


(A-38) 


The  wind-axis  roll  rate  is  the  first  component  of  the 
body  angular  velocity  expressed  in  wind  axes,  which  is 


=  Hg(a ,p )  wB 


(A-39) 


Equations  relating  the  body  angular  rate  in  wind  axes  to  the 
aerodynamic  and  velocity  angular  rates  can  be  derived  with 
reference  to  Fig.  A-l. 


%  =  “v<'*v>hIh(^  0  +  H>v>  i 
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+  H  ^  ( -  P  ) 
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which  results  in 


Pw 

-  (siny  +  asinp 

qw 

(cosy  sin<t>v  +  YC°S(t>v  +  acosp 

Jw. 

(cosy  cos4)v  y s inif) v  -  p 

(A-40) 


(A-41) 


These  last  equations  are  useful  in  deriving  simple  ACM  pilot 
models  and  setting  up  point-mass  aircraft  trajectory  genera¬ 
tors.  The  earth-relative  velocity  and  the  inertial-velocity 
angles  are  also  candidates  for  inclusion  in  command  systems. 
The  equations  are  given  below. 


V 

h 

*  Hg (4»  »e  ,0 ) 

i 

3  > 

i _ 

nr 

i _ 

LwJ 
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+  y?  +  *i 
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= 

arc  cos  +  y^  j 

.  /  •  /•  2  ,  •  2  .  *2  \ 

Y 

arc  sm  [~z,  -Jx,  +  y  +  z  T  1 

L  J 

L  \  1  1  1  ij 

(A-42 ) 


(A-43) 


For  a  fighter  pilot  requiring  rapid  sustained  orieri 
tation  changes,  an  acceleration-oriented  maneuvering  set, 


=  fv.  a  ,  a,  P,  a  ,  p  1  T 


(A-44) 


could  be  useful.  The  maneuvering  set  gives  the  pilot  direct 
control  over  normal  accelration,  an,  and  roll  rate  about  the 
velocity  vector,  p^.  Independent  fuselage  pointing  is  pro¬ 
vided  about  the  velocity  vector  using  angle  of  attack,  a,  and 
sideslip,  p,  commands.  The  air-relative  velocity  magnitude, 
V,  is  commanded,  and  the  aircraft  can  be  directed  to  make  a 
flat  turn  (no  bank  angle)  with  the  lateral  acceleration,  ay, 
command . 


A. 3  LINEAR  DYNAMIC  EQUATIONS 

For  aircraft  stability  and  control  response  analysis, 
linearized  versions  of  the  nonlinear  state  and  output  equa¬ 
tions  are  necessary.  The  linearized  aircraft  state  equations 
are  given  in  Ref.  23,  along  with  a  description  of  the  linear¬ 
ized  aerodynamic  force  and  moment  relationships.  Lineariza¬ 
tion  of  the  output  equations  described  in  Section  A. 2. A  is 
given  here. 

The  command  augmentation  system  design  methods  of 
Chapter  4  require  linearized  versions  of  the  output  vector, 
which  is  a  function  of  both  aircraft  states  and  controls,  so 
the  following  perturbation  output  vector  equation  results: 

Ayd  =  +  D^0'iio'1  (A-45) 

The  individual  rows  of  H  and  D  depend  on  the  chosen 
elements  of  the  command  vector,  and  the  following  equations 
describe  the  available  constituents.  The  perturbation  wind- 
axis  velocity  vector  is  related  to  the  perturbation  body-axis 
velocity  vector  as 
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AV 

AP 

Ad 


<V0.P0)  h”<« 


O 


P0)  AXB 


(A-46 ) 


Jw  is  a  diagonal  matrix  which  has  elements  1,  V  ,  and  VQ  cos 


The  following  linearized  output  equation  gives  the 
relationship  between  the  body-axis  state  variables  and  the 
perturbation  wind-axis  accelerations: 


AV 

Aa> 

Aa 


n 


s  HB(ao'^o)  )A¥B  +  A^B  '^B 

o  o  B 


iB  +  °B  — B 
0  0  0 


*v<“o>jw1<vo^>hb(VI,o>as:b 


(A-A7 ) 


where 


Lw<%> 


0  sin  a  0 

0  0  1 

0  -cos  a  0 

o 


(A-48 ) 


These  equations  are  easily  evaluated  using  general  computer 
routines  that  have  been  developed  for  this  type  of  analysis. 
Equation  A-47  requires  both  the  nominal  and  perturbation  body- 
axis  velocity  derivatives,  Vg  and  AVg.  Vg  is  part  of  the 

nominal  flight  condition  specification,  while  AVg  consists  of 
three  rows  of  the  linear  system  differential  equation.  Intro¬ 
ducing  these  three  rows  causes  the  accelerations  to  be  functions 
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of  the  perturbation  Euler  angles,  body-axis  translation¬ 
al  and  angular  rates,  and  the  perturbation  control  deflections. 


The  perturbation  wind-axis  roll  rate  depends  on  both 
the  perturbation  body-axis  angular  rate  and  on  the  perturba¬ 
tion  body-axis  velocity,  which  affects  the  body-to-wind  axis 
transformation  matrix.  The  linearized  equations  for  body  an¬ 
gular  rate  in  wind  axes  are  as  follows 


Ap, 

Aq, 

Ar. 


w 


■w 


0 


(A-49) 


Linearization  of  the  wind  axis  quantities  V,£,y,  and  <i>v  can 
best  be  pursued  by  regarding  them  as  the  state  rates  in  a  sys¬ 
tem  whose  input  consists  of  the  commanded  values  of  V,  an ,  ay 
and  pw.  Assume  =  0.  The  linear  state  equations  are 
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An  additional  portion  of  the  linearised  relationships 
deals  with  the  aircraft  control  system.  Only  the  yaw  sas  is 
modeled  in  the  pilot-model  testing  work,  as  the  roll  cas  is 
off  and  the  longitudinal  dynamics  are  not  analyzed.  As  shown 
in  Fig.  23,  the  yaw  sas  feeds  back  washed-out  yaw  rate  and 
low-pass-filtered  lateral  accelerometer  output.  The  low-pass 
filter  bandwidth  is  very  large  and  hence  Lhe  filter  i9  modeleo 
as  a  simple  gain.  Hence,  the  yaw  sas  adds  one  state  to  the 
aircraft  model.  The  yaw  sas  actuator  deflection  is  given  by 
the  equation 

5r,sas  =  kry  *  ^  r  '  k2  AXsaS 
and  the  yaw  sas  state  equation  is 

'“‘sas  *  <-V»>AxM1  *  <1/x)r 

where  x  is  the  yaw  rate  washout  time  constant,  k-,  is  the  lat¬ 
eral  acceleration  feedback  gain  and  k£  is  the  yaw  rate  feed¬ 
back  gain. 


(A-51) 


(A-52 ) 


163 


APPENDIX  B 


AIR  COMBAT  EQUATIONS  OF  MOTION 


The  full-order,  nonlinear  equations  of  motion  for  the 
tracking  dynamics  model  in  the  air  combat  scenario  are  developed 
in  this  section.  Simplified  linear  equations  of  the  nonlinear 
tracking  dynamics  model,  amenable  to  computer  analysis,  are 
also  derived  for  use  in  Chapter  4.  The  addition  of  pipper 
dynamics,  although  not  needed  in  this  work,  is  a  straightfor¬ 
ward  extension  of  this  approach. 

The  target  racking  angles  are  defined  as  the  pilot's 
view  angles  relative  to  the  aircraft  x-axis.  These  angles,  eR 
and  Cy,  are  shown  in  Fig.  44b.  Derivation  of  the  tracking 
angle  dynamics  best  occurs  in  the  wind  axis  system  because  it 
is  aligned  with  the  velocity  vector.  By  expressing  the  target- 
to-tracking  aircraft  displacement  vector  in  wind  axes  and  divid¬ 
ing  by  the  target-aircraft  range,  the  tracking  angle  equation 
result'’ 


(B-l ) 


where,  for  relatively  small  wind  axes  angles, 


--WT 


R/R 

-«-y  +  e> 

-Un  +  a) 


(B-2 ) 


The  derivative  of  the  tracking  angle  equation  follows  directly 
(assuming  R  =  0) 


3  64 


(B-3 ) 


,T  "  -1^} 


The  aircraft  velocity  in  wind  axes,  by  definition,  is 


( B  -  4 ) 


and  one  can  define  the  target's  velocity  relative  to  the  track¬ 
ing  aircraft  velocity  as  follows 


The  second  term  in  Eq.  B-3  can  be  redefined  in  terms  of  a  cross- 
product  of  the  tracking  angle  vector  (vVT)  and  the  wind-axis 
angular  rates.  The  lateral  and  normal  tracking  angle  dynamic 
equations  finally  result 

ey  =  -  <5Vy/R  +  (£n  +  “><PW  ’  asinp)  +  ay/V  -  p  (B-6) 

and 


V  +  <5V. 

6V. 


x 


(B-5 ) 


En  =  ‘  aVz/R  '  (ey  +  (Pw  ‘  «sinp)  -  an/V  -  a  (B-7) 

Equations  for  the  dynamics  of  6Vy  and  6VZ  can  be  derived  by  a 
similar  process.  The  results  are 


6Vv  s  anT  sin ~  a  +  ,(6VZ)(PW  -  asinp) 


(6Vx)(ay/V) 


(B-8 ) 
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(B-9) 


6VZ  =  -  anT  cos(6<I)v)  +  an  -  (6Vy)(pw  -  a  sin  p) 

+  (<5VX)  <an/V) 

where  64>v  is  defined  as  shown  in  Fig.  44b. 

An  attempt  was  made  to  calculate  the  target  accelera¬ 
tion  (an,p)  and  roll  (<j>v,p)  time  history  from  the  available  data 
(en,  ey,  an,  ay,  a,  p)  but  was  unsuccessful  due  to  the  amount 
of  differentiation  of  the  input  time  histories  necessary. 
Because  of  this  difficulty  the  tracking  equations  are  simpli¬ 
fied  by  assuming  that  <5V  and  p  -  asinp  are  identically  zero. 

A  W 

Further  simplification  occurs  by  noting  that  near  steady  state 
the  perturbation  velocities  can  be  expressed  as 

6Vy  =  RUy/v)  +  6^y  At  (B'10) 

6VZ  =  R(-an/V)  +  6Vz  At  (B-ll) 

Using  R/V  as  the  time  interval  At  in  Eq.  B-10  and  B-ll  allows 
these  short  term  approximations  to  be  inserted  into  Eqs.  B-6 
and  B-7.  The  derivatives  6V  and  6VZ  can  then  be  eliminated 
from  Eqs.  B-6  and  B-7  by  inserting  Eqs.  B-8  and  B-9.  This 
produces  the  final  simplified  nonlinear  equations. 

ey  =  (-1/V)  (anTsin(6<j)y)  -  ay)  -  p  (B-12) 

en  =  (1/V)  (anT  cos(60v)  -  an)  -  a  (B-13) 

These  equations  enable  the  identification  of  the  target  tra¬ 
jectory,  and  also  form  the  basis  for  the  final  linearized  equa 

tions,  which  assume  an>p  =  an  ,  and  641  =  0.  These  linearized 

1  0  vo 

equations  are 

Ac v  =  <-i/V  )(a  -  )(A64»  )  +  (l/VQ)Aa  -  A(3 
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(B-14) 
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l  ‘  -  ian>  ’  L*  <B-15> 

f 

,s  The  lateral  tracking  equation,  Eq.  B-1A,  is  used  in  Chapter  A 

ar  in  the  construction  of  an  air  combat  pilot  model. 
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APPENDIX  C 
HYPOTHESIS  TESTING 


This  appendix  gives  the  detailed  equations  for  the 
hypothesis  testing  procedure  used  in  Chapter  4.  The  procedure 
is  based  on  a  model  structure  identification  algorithm  modi¬ 
fied  to  allow  for  changes  in  the  true  hypothesis. 

The  hypothesis  testing  philosophy  for  investigating 
pilot  control  adaptation  in  flight  is  to  construct  a  number  of 
different  pilot  models  and  then  determine  which  pilot  model 
best  represents  the  pilot  control  behavior.  The  i1"*1  pilot 
model  is  represented  by  the  discrete-time  equation 


-Aul 

■■  mm 

Au 

4XE 

=  *PM 

AXE 

+  rPM^k  + 

Az 

Az 

k+1  k 


(C-l) 


The  pilot  model  is  taken  from  Ref.  23  where  the  pure  time  delay 
in  the  optimal  control  pilot  model  is  replaced  by  a  Pade  approx 
imation.  The  steps  for  arriving  at  Eq.  C-l  are  discussed  in 
Section  4.1.  In  Eq.  C-l,  Au^  is  the  control  predicted  by  the 
pilot  model,  Ax^.  is  the  state  vector  of  the  pilot  model's  inter 
nal  Kalinan  filter,  Az^  are  states  associated  with  the  Pade 
approximation ,  A£  ar?  the  assumed  observations  of  the  pilot, 
and  Aw^  is  the  process  noise  of  the  pilot  model. 

The  method  for  determining  which  pilot  model  provides 
the  best  fit  at.  any  given  time  is  taken  from  Refs.  30  and  31 
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and  illustrated  in  Fig.  C-l.  The  notation  used  in  the  figure 
is  as  follows: 


Au  measured  value  of  pilot  perturbation  con- 

'”,,K  trol  at  time  t^ 

Ay.:  measured  values  of  perturbation  states 
assumed  observed  by  pilot 


Auj^(-):  one-step  predicted  value  of  Aum  ^  based 

upon  |Aum  l,  ....  and  (AZl . 

Ayb  , }  for  the  iz  pilot  model  Kalman 
filter 

Ar^:  one-step  predicted  residual  at  the  kL 

step  for  the  i  pilot  model  Kalman  filter 

P„  probability  that  the  i^  pilot  model  is  best, 

’  based  upon  knowledge  of  {Au  ,  ...»  Au  .  ,} 

and  {Ajrj . A^.j}  -“>1>  "“■k  1 


The  pilot  model  Kalman  filter  is  constructed  using  Eq.  C-l  as 
the  plant,  ^Aw^  as  the  plant  process  noise,  Ay^  as  the  known 
plant  inputs,  and  the  observation  is 


4Vk  =  11  0  01 


"Au  " 

a*e 

Az 

m  m 

k 


+  A^k 


<  C-2 ) 


For  our  purposes,  Av^  is  the  measurement  noise  encountered 
when  attempting  to  measure  the  actual  pilot  controls.  The 
data  analyzed  in  Chapter  A  exhibited  very  low  measurement  noise 
on  the  pilot  controls  and  Av^  is  assumed  to  be  zero.  The  pilot 
model  Kalman  filter  just  described  is  not  to  be  confused  with 
the  internal  Kalman  filter  in  the  pilot  model. 
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Figure  C-l 


PROBABILITY 

CALCULATION 


i  ter  Approach  to  Pilot 
Concho.  rategy  Hypothesis  Testing 


The  output  of  the  i1'  pilot  model  filter  is  a  state 
esti’1  ite  [AuT(-),  AzT(-)]^  of  the  pilot's  internal 

peri^rbation  states  and  controls  based  upon  the  first  k-1 
input-output  pairs, 
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The  update  part  of  the  pilot  model  Kalman  filter  is 
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(C-5 ) 


The  matrix  K1 
filter  gain  for  the  i ^ 


is  the  steady-state  pilot  model  Kalman 
pilot  model. 


For  notational  purposes,  let 

Uk  s  {A^i»,l . (C“6) 

Applying  Bayes'  rule,  one  can  write: 


PH,i(k)  “  p<H,i|uk) 

*  P(H , i  | tk *  Uk.!> 

_  p^!!„.k|H'i.uk.l>P<H»1luk.1) 
P<iVklUk-l> 

,  pUa  m,kH’1’ulk-i):p(H’1)uk-i 

E  PWH11,k|H.i.Uk.1)P(H,i  Uk.1> 
i=l 

L  P<AHm,k|H,i,Uk.1)PHti(k-l) 
i-1 


171 


Using  the  fact  that  Aw^  is  assumed  to  be  a  zero-mean  gaussian 
sequence,  along  with  the  assumption  the  Kalman  filter  state 
vector  |Au^(-),  AXp(-),  Az^-)]^  is  also  gaussian,  it  follows 
that 


P<AHm,U  « > i * Uk_ x )  = 


(2n)m/2  det 


■3s(Arj)T(Si)"1(ArJ) 


where 


S1  =  [I  0  0]  Pi(-) 
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0 
L  o 


(c-S) 


(C-9) 


P*(-):  a  priori  steady  state  erro^  covariance 

t  h 

matrix  computed  using  i  pilot  model 
Kalman  filter 


Note  that  the  PH  ^(k)  are  computed  recursively.  In 
addition,  it  should  be  mentioned  that  the  approximation  shown 
in  Eq.  C-8  is  based  on  the  assumption  that  the  "true"  hypothe¬ 
sis  remains  constant.  If  the  true  hypothesis  changes  with 
time,  the  probability  density  function  in  Eq .  C-8  must  be  con¬ 
ditioned  on  the  hypothesis  values  at  all  previous  times.  This 
leads,  as  time  progresses,  to  an  expanding  "tree"  of  hypotheses. 
It  is  necessary  to  calculate  Pu  .(k)  for  each  of  these  hypothe- 
ses.  Since  this  imposes  an  excessive  computational  burden,  an 
efficient  suboptimal  approach  to  the  problem  is  developed  in 
Ref.  30.  The  suboptimal  approach  is  nearly  optimal  in  situa¬ 
tions  where  the  probability  of  switching  hypotheses  at  any 

given  time  is  low.  The  suboptimal  approach  is  derived  by  not- 

t  h 

ing  that  if  the  i  pilot  model  is  true  then  the  recursive  cal¬ 

culation  in  Eq .  C-7  causes  P,,  .  to  approach  1.  If  Pu  .  should 

ri ,  1  n  1 1 
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equal  one,  it  will  remain  "trapped"  at  one  regardless  of  changes 
in  the  pilot's  control  strategy.  To  prevent  this,  Eqs .  C-2  and 
03  are  modified  in  the  following  fashion.  Define  a  parameter 
e  corresponding  to  the  minimum  value  that  one  of  the  ^(k) 
can  fall  to  and  still  "recover"  if  H,i  suddenly  becomes  the 
true  hypothesis,  (In  the  simulation  studies  reported  in  Ref. 

30,  £=0,05  was  found  to  be  a  good  value  and  is  used  in  Chapter 
4).  Define 


6<k)  4  £  [e  -  PH  (k)]  (010) 

{i:PH>.(k)<£}  H 


A 

m(k)  =  number  of  ..(k)  that  are  greater  than  2e 

(C-ll) 


where  the  symbology  in  the  equation  for  6 ( k >  is  meant  to  de¬ 
note  the  fact  that  the  sum  is  taken  over  all  i  such  that  PH  ^ 

<  e.  First  update  the  P„  .  (k)  using  Eqs.  07  and  08  then  reset 

n  j  1 

these  probabilities  as  follows: 


• 

If  PH , : 

j (k)  <  e , 

set  PH  A(k)  “ 

e 

• 

If  PH,tOO  >  2e 

and  6 ( k )  >  0, 

set 

PH,i<k)  * 

[PH.i<k>] 

6(k) 

m(k) 

(C-12 ) 

■  ^ 

old 

where 

denotes  the 

value 

obtained 

from  Eq.  07.  old 

The  reset  procedure  guarantees  that  all  of  the  P^  ^(k)  are 
greater  than  (or  equal  to)  e  and  that  they  sum  t.o  one.  Re¬ 
turning  to  the  previous  example,  the  problem  of  P^  ^(k)  be¬ 
coming  so  small  that  it  cannot  recover  when  the  true  hypothe¬ 
sis  swi'ches  to  H..  is  now  avoided.  In  addition,  the  procedure 
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can  be  implemented  with  only  N  parallel  filters.  The  simula¬ 
tion  results  in  Chapter  4  indicate  that  this  suboptimal  proce¬ 
dure  is  effective  in  detecting  changes  in  the  true  hypothesis. 

The  above  hypothesis  testing  procedure  can  be  generalized, 
th 

The  i  pilot  model  matrices  in  Eq.  C-l  can  be  time-varying  and 
the  pilot  model  can  even  be  nonlinear.  The  time-varying  Kalman 
filter  can  be  used  in  the  former  assumption  and  the  extended 
Kalman  filter  can  be  used  in  the  latter  assumption. 


APPENDIX  D 
LIST  OP  SYMBOLS 


In  general,  matrices  are  represented  by  capital  letters 
and  vectors  are  underscored;  exceptions  to  these  rules  are 
only  made  when  they  are  contradicted  by  standard  aerodynamic 
notation.  Capital  script  letters  are  used  to  denote  scalars  in 
some  cases. 


VARIABLE  DESCRIPTION 

a  Control  deflection  weighting  scalar 

a  Normal  acceleration  (inertial  acceleration 

a  normal  to  the  velocity  vector  in  the  body 

x-z  plane  -  including  gravitational  accelera¬ 
tion) 

b  Wing  span 

C  Gain  Matrix  (Type  1) 

Pilot  control-strategy  feedback  matrix 

Cio  Partial  derivative  of  the  nondimensional  co- 

A  efficient  of  force  or  moment  1  with  respect 

to  the  nondimensional  variable  2  (scalar) 

c  Mean  aerodynamic  chord 

D  Control  observation  matrix 


F  System  dynamics  matrix 

F  Aerodynamic  contact  force  vector 

f  Vector-valued  nonlinear  function  for 

system  dynamics 

G  Control  input  allocation  matrix 

G  Thrust  moment  vector 

g  Magnitude  of  g  (=  |  £  |) 
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&  Gravitational  acceleration  vector 

Vector-valued  nonlinear  function  for 
output  equations 

H  Aircraft  state  observation  matrix 

Hv  Total  pilot  observation  matrix  for  pilot 

model 

2 

Hh  Euler  angle  transformation  from  Frame  1  axes 

to  Frame  2  axes 

h  Altitude 

h  Angular  momentum  vector 

I  Identity  matrix 

Ig  Rotational  inertia  matrix 

J  Cost  functional 

K  Gain  matrix  (Type  0) 

Pilot  Kalman  filter  gain  matrix 

L  Angular  velocity  transformation  matrix 

Aerodynamic  moment  about  the  x-axis  (scalar) 

1  Number  of  outputs 

M  Aerodynamic  moment  about  the  y-axis  (scalar) 

Weighting  matrix  on  state-control  rate  and 
control-control  rate  products 

M  Aerodynamic  contact  moment  vector 

m  Mass  of  the  vehicle 

Number  of  controls 

Number  of  pilot  model  probabilities  that 
exceed  a  specified  level 

N  Aerodynamic  moment  about  the  z-axis  (scalar) 

n  Number  of  states 

P  Riccati  matrix  in  the  optimal  regulator 

problem 

Fj,  Probability  that  a  pilot  model  among  a 

M  set  of  pilot  models  is  corrent 


176 


P  Pilot  noise-to-signal  ratio  for  neuromotor 

noise 

p  Rotational  rate  about  the  body  x-axis 

Q  State  weighting  matrix 

q  Rotational  rate  about  the  body  y-axis 

1  2 

Free  stream  dynamic  pressure  ) 

R  Control  or  control-rate  weighting  matrix 

R,  Matrix  with  diagonal  consisting  of  the  in¬ 

verse  of  human  neuromuscular  time  constants 

r  Rotational  rate  about  the  body  x-axis 

r  Kalman  filter  innovations  vector 

S  Reference  area  (usually  wing  area  for 

aircraft ) . 

Covariance  matrix  of  Kalman  filter 
innovations  vector 

T  Thrust  force  magni  tude(  |  T  j) 

T  Thrust  force  vector 

t  Time 

U  A  set  of  control  vectors  ordered  in  time 

u  Body  x-axis  velocity  component 

u  Control  vector 

uc  Pilot  model  control  command 

V  Inertial  velocity  magnitude  (|V|) 

Covariance  matrix  for  pilot  control 
measurement  noise 

V  Velocity  vector  of  body  observed  from  iner¬ 
tial  axes 

Vu  Pilot  neuromotor  noise  covariance  matrix 

V  Pilot  observation  noise  covariance  matrix 

v  Body  y-axis  velocity  component 
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V 


^u 


Y 

y 

1 

Z 


z 


z 


VARIABLE 

(GREEK) 

a 

P 


r 


Pilot  control  measurement  noise  vector 

Pilot  neuromotor  noise  vector 

Pilot  observation  noise  vector 

Eody  z-axis  velocity  component 

Total  disturbance  vector  (including  pilot 
noise ) 

Aerodynamic  force  along  the  x-axis  (scalar) 
Position  along  the  x-axis 
State  vector 

Longitudinal  distance  between  actual  c.g, 
location  and  point  used  for  aerodynamic 
moment  measurements  (expressed  in  body 
axes ) 

Aerodynamic  force  along  the  y-axis  (scalar) 
Position  along  the  y-axis 
Pilot  observation  vector 
Lagged  pilot  observation  vector 
Aerodynamic  force  along  the  z-axis  (scalar) 
Position  along  the  z-axis 

* 

Lagged  states  associated  with  the  Paae 
approximation 

DESCRIPTION 

Wind-body  pitch  Euler  angle  (angle  of  attack) 

Negative  of  wind-body  yaw  Euler  angle 
(sideslip  angle) 

Discrete  time  control  effect  matrix 


Y  Inertial-velocity  axis  pitch  Euler  angle 

(flight-path  angle) 
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The  number  of  pilot  model  probabilities 
that  are  less  than  a  specified  level 

Differential  stabilator  deflection 

Maneuver  flap  deflection 

Rudder  deflection 

Symmetric  or  collective  stabilator  deflection 
Spoiler  deflection 
Tracking  error 
Damping  ratio 

Inertial-body  pitch  Euler  angle 
Eigenvalue 

Aircraft  orientation  vector 

Inertial -velocity  axis  yaw  Euler  angle 
(velocity  heading  angle) 

Correlation  coefficient 

Air  density 

Real  part  of  an  eigenvalue  in  radians/sec 
Pilot  time  delay 

Human  neuromuscular  time  constant 
Discrete  time  state  transition  matrix 
Inertial -body  axis  roll  Euler  angle 
Stability-axis  roll  angle 
Inertial-body  axis  yaw  Euler  angle 
System  inverse  matrix 

Frequency  in  radians/sec;  imaginary  part 
of  an  eigenvalue 


—1 

Rotational  rate  vector  of  Reference  Frame 

2  with  respect  to  Reference  Frame  1  and 
expressed  in  Frame  1  coordinates. 

2  11  1 

so  E 2  *s  left-handed.  Thus , 

Frame  1  and  Frame  2  are  not  interchange¬ 
able.  ) 

VARIABLE 
(SUBSCRIPT  OR 
SUPERSCRIPT 

DESCRIPTION 

a 

Aircraft  dynamics 

B 

Body  axes 

c 

Commanded  value 

D 

Discrete  time  matrices  for  optimal  control 

d 

Desired  value 

E 

Earth  axes  (north,  east,  down) 
Conglomerated  pilot  model  matrices 

I 

Inertial  axes 

i 

Pilot  model  index 

k 

Time  index  for  discrete  time  quantity 

SL 

Aerodynamic  moment  about  the  x-axis 

m 

Aerodynamic  moment  about  the  y-axis 

n 

Aerodynamic  moment  about  the  z-axis 

0 

Nominal  value  about  which  linearization 
occurs 

PM 

Pilot  Model 

s 

Psuedoinverse  weighting  matrix 

T 

Target  dynamics 

V 

Velocity  axes 

W 

Wind  fxes  (same  as  stability  axes  for 

P0  s  ff0  =  0 ^ 

180 


X 


Component  along  the  x.-axis 
y  Component  along  the  y-axis 

z  Component  along  the  z-axis 

X  Aerodynamic  force  along  the  x-axis 

Y  Aerodynamic  force  along  the  y-axis 

Z  Aerodynamic  force  along  the  z-axls 


OPERATOR 

C) 

r> 


(-) 


(+) 

r> 


DEFINITION 


Time  derivative 

Matrix  equivalent  to  vector  crossproduct. 
Specifically,  if  x  is  the  three-d imensit nal 
vector 


X  = 

HI  «■ 

X 

y 

,  then  x  = 

"  0 
z 

“Z 

0 

y 

“X 

z 

_*y 

X 

0 

and  the  cross  product  of  x  and  f  is  equal 
to  the  product  of  the  matrix  x  and  the 
vector  f, 

x  x  f  =  xf 

The  one  step  predicted  estimate  of  a 
quantity  immediately  before  meas- 
surement 

The  best  estimate  of  a  quantity  immediately 
after  a  measurement 

Difference  between  a  vector  and  its  expected 
value 


(A)  Estimated  value 

T 

(  )  Transpose  of  a  vector  or  matrix 

(  ) * ^  Inverse  of  a  matrix 
(  )'  Limited  value 

(  )  Perfect  model  following  values 


181 


A<  ) 
6(  ) 


ZLi 

c> 

ACRONYM 
ACM 
ARI 
cas 
CGT 
Cl  LOP 
c.g. 
DFCS 
DOF 
DMS 
NASA 
NL 
PI 
PIO 
sas 
VALT 


Perturbation  about  the  nominal  value  of 
a  variable 

The  variable  is  the  difference  between  the 
target  variable  and  the  aircraft  variable 

Exponential  of  (  ) 

Right  pseudoinverse  of  a  matrix 

Maximum  value,  usually  due  to  displace¬ 
ment  of  an  actuator. 

Clockwise  angular  rotation 

Nondimensional  or  average,  value 

CORRESPONDING  PHRASE 

Air  Combat  Maneuvering 

Aileron-Rudder  Interconnect 

Command  Augmentation  System 

Command  Generator  Tracker 

Conversion  in  lieu  of  procurement 

Center  of  Gravity 

Digital  flight  control  system 

Degrees  of  Freedom 

Differential  Maneuvering  Simulator 

National  Aeronautics  and  Space  Administration 

Nonlinear 

Proportional  Integral 
Pilot  Induced  Oscillation 
Stability  Augmentation  System 
VTOL  Approach  and  Landing  Technology 
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